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ABSTRACT 


An experimental study of gas film cooling was conducted on a 3X size 
model turbine vane. Injection in the leading edge region was from a 
single row of holes angled in a s panwise direction. The suction surface 
injection configurations investigated were: injection angle, 3 - 18°, 35 D i 

location of holes from stagnation, x-/L s - 0.083, 0.150, 0.228. Experi- 
ments were conducted at reduced freestream pressure and temperature to 

5 

simulate a leading edge Reynolds number, Re^ = 1.6 x 10 and a coolant- 
to-mainstream density ratio of Pq/p^ ~ 2.18. Measurements of the local 
heat flux downstream from the row of coolant holes, both with and without 
film coolant flow, were used to determine the film cooling performance 
presented in terms of the Stanton number ratio, STp C /ST Q . Results for a 
range of coolant blowing ratio, M = 0 to « 2.0, indicate a reduction in 
heat flux of up to 15 to 30 percent at a point 10 to 11 hole diameters 
downstream from injection. An optimum coolant blowing ratio corresponds 
to a coolant-to-freestream velocity ratio in the range of 0.5. The shallow 
injection angle resulted in superior cooling performance for injection 
closest to stagnation, while the effect of injection angle was insignificant 
for injection further from stagnation. 


NOMENCLATURE 


area; constant in cooling effectiveness correlation 

film cooling correlating parameter 

exponent for (x '/MS) in effectiveness correlation 

exponent for coolant Reynolds number in effectiveness 
correlation 

drag coefficient for cylinder-in-crossflow 
skin friction coefficient 
geometric modeling factor 
specific heat at constant pressure 
specific heat at constant volume 

correction factor to account for real flov/ effects in slug 
flow analysis 

coolant hole diameter 

turbulence screen wire diameter 

Eckert number, V 2 /C.,{T t - T 1( ) 
w p vr 

force 

Grashof number, p 2 g£ {Ty - TJL^/jj 2 

' c/i 

local acceleration of gravity 

constant relating units of force, mass and length, 
ci - 32.17 lb -ft/1b. r sec^ 

C i-.i i 

test section duct height 

boundary layer shape factor 

local heat transfer coefficient, q'VO- - T, t ) 


NOMENCLATURE (cont'd) 


local heat transfer coefficient with film cooling and 
T C * V q FC^ T AM " T W,FC } 

laterally averaged heat transfer coefficient with T c ; 

2 2 

momentum flux ratio, PqJq/pJJ^ 

angle between surface and gradient of gravitational 
potential field 

o 

acceleration parameter, (o^/V^) dVV /dx 
proportionality constant 
characteristic length? length of coolant slot 
length of film cooling hole 
shear-layer length of separation bubble 
reattachment length of separation bubble 
blowing ratio, p c V c /p 0 V m 
M ach number 

mass flux of entrained fluid into coolant dnt 
turbulence screen mosh spacing 
meter (length or dimension) 
mass flowrate 

film cooling variable, function of B 

number of film cooling holes 

unit normal vector 

pressure 

pressure 

dimensionless pressure, p/p ( ..V^ 

Prandtl number 
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weight flow of coolant per unit width of plate 
heat flux 

characteristic radius 
radius of curvature 
Reynolds number 
gas constant 
coolant slot height 

spacing between coolant hole centerlines 

Stanton number 

temperature 

average boundary layer temperature 
wall thickness 

time averaged streamwise velocity 

fluctuating velocity component in streamwise direction; 
dimensionless velocity In x-direction* u/V w 

velocity 

velocity component in vertical direction 
dimensionless velocity in y-direction, v/V , 
test section duct width 
coolant weight flow 

film cooling parameter according to Spalding 

distance in streamwise direction 

distance downstream of cooling holes, -'x-x..; also x/L 

i 

dimensionless streamwise distance frem cooling holes, x/d 
vertical distance from surface 
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NOMENCLATURE (cant'd) 

dimensionless vertical distance, y/d^; also y/L 
penetration distance cf coolant jet into boundary layer 
spanwise or lateral coordinate 
dimensionless spar, wise distance, z/d^j also z/L 

streamwise injection angle; thermal diffusivlty 

dimensionless coolant jet trajectory parameter 

component of compound injection angle measured from surface 

component of compound injection angle measured clockwise 
from positive x-axis in x-z plane 

spanwise injection angle 

ratio of orifice or venturi diameter to pipe diameter 

coefficient of thermal expansion 

distance from effective origin of boundary layer 

distance from boundary layer origin to point where cooling 
effectiveness is unity 

ratio of specific heats, C. 5 /C y 

boundary layer thickness 

displacement boundary layer thickness 

cooling effectiveness, identical to n 

point source thermal diffusivity 

cooling effectiveness parameter 

dimans ion! ess coordinate normal to surface 

momentum boundary layer thickness 

dimensionless temperature (T-j- - T)/(Y-j » f w ) 


xi 

NOMENCLATURE (ccnt'd) 

U dynamic viscosity 

v kinematic viscosity 

£ dimensionless streamwise coordinate; coordinate along 

coolant jet axis 

p density 

p' dimensionless density 

$ viscous dissipation function 

$' heat flux ratio, qp C /qg 

X angle between x-axis and tangent to coolant jet trajectory 

oi power of temperature ratio for gas viscosity variation 

with temperature 

to 1 frequency 

SUBSCRIPTS 

ADW adiabatic wall condition 

B bul k 

C cool ar.t 

CHAR characteristic 

FC film cooling condition 

i injection 

ISO isothermal condition 

IVJ impervious wall condition 

0 jet 

L lower 

LE leading edge 

0 initial condition 
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nDMEriCLATUflE (cont’d) 


OPT optimum condition 

0PT,NEAR optimum condition for q" measurement location nearest the 
injection holes 

P pressure surface 

RE at reattachment 

REF reference 

RMS root-mean- square 

5 suction surface 

SEP at separation 

SP span 

T total condition 

U upper 

W wait 

X in stream*, vise direction 

0 without film cooling 

frees tr earn 

A based on boundary layer energy thickness 

6 based on boundary layer momentum thickness 


I. INTRODUCTION 


1.1. General Discussion 

In many engineering applications, the successful performance of an 
overall system depends upon maintaining tolerable temperature levels at 
various times and locations within the working cycles of the individual 
components* This is especially true in the gas turbine industry as tur- 
bine efficiency and work output increase as the maximum cycle operating 
temperature is elevated. Consequently, to survive the severe thermal 
environment created by high combustor exit temperatures, the initial 
stages of the turbine must be cooled. Convective cooling, internal to 
turbine vanes and blades, has been used extensively; however, increased 
operating pressures and temperatures require external cooling to augment 
internal cooling configurations (1) . Gas film cooling is one such 
external cooling technique. 

Film cooling is accomplished by injecting, at discrete locations, a 
relatively cool, secondary fluid into the boundary layer of a surface 
exposed to a hot, primary mainstream flov/. The injected coolant film is 
analogous to an external heat sink, providing a blanket protection that 
progressively deteriorates as the coolant film moves farther from the 
site of injection. 


■r 

Numbers in parentheses refer to REFERENCES. 
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The present investigation is concerned with determining the cooling 
effectiveness of a specific gas film cooling injection geometry located 
in the leading edge region of an inlet turbine vane* The following two 
sections will discuss the general film cooling problem, its relevance to 
gas turbine design and the impetus behind the present choice of injection 
configurations and parameters. 

1*2. Statement of the Problem 

Since trends point toward future gas turbine operating pressures 
and temperatures up to 4.05*lG 6 N/m 2 (40 atm) and 2200°K (350CPF), there 
is a need to thoroughly understand and improve upon current film cooling 
techniques in order to adequately augment internal cooling schemes. 

The basic problem associated with film cooling is the determination 
of the degree to which a particular injection configuration can effec- 
tively reduce the local heat flux conditions downstream of injection. 
Figure 1 illustrates the basic gas film cooling configuration with a 
relatively cool fluid injected into the hot gas boundary layer at x - x. . 
With zero film coolant flow, the surface heat flux for prescribed hot gas 
and wall temperatures is qj. With film coolant injection, the coolant 
film acts as an effective heat sink, reducing the surface heat flux down- 
stream from injection to q^. Thus, external film cooling reduces the 
amount of internal heat removal necessary to maintain the prescribed 
wall temperature for x > x^. As the cool fluid elements move farther 
from injection the coolant film loses its "identity" as a protective 
layer due to an increase in the fluid element temperature and mixing with 
hot gas stream. Consequently, the reduced heat flux, qj^* approaches a 



condition characteristic of the normal, uncooled boundary layer heat 
transfer for large distances downstream from injection. The primary 
objective of any film cooling study is to determine how effective a 
given coolant injection geometry is in reducing the local heat flux 
under a range of coolant and freestream flow conditions in order to 
identify configurations and flow conditions for optimum or most effec- 
tive film cooling performance. 

The subject investigation is concerned with a twofold problem in 
relation to gas film cooling for application to the design of high pres- 
sure and high temperature gas turbine engines: 

(1) effective cooling of the high heat flux leading edge region 
of an inlet turbine vane subject to geometric injection con- 
straints which may be imposed by significant leading edge 
surface curvature, and 

(2) determining the heat flux reduction characteristic of the 
configuration implied in (1) at flow conditions that simulate 
the actual turbine environment. 

The manner in which these two problems are confronted in the current 
study will now be discussed. 

1.3. Research Objecti ves 

The purpose of the subject investigation is to provide meaningful 
information on gas film cooling of turbine vanes in the form of corre- 
lations for the local surface heat flux reduction for use in the design 
of high temperature and high pressure gas turbine engines. 
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Film cooling configurations are considered that offer a practical 
solution to reducing the heat load in the leading edge region of an inlet 
turbine vane. Some of the most severe heat loads within the entire work- 
ing cycle of the gas turbine engine are experienced at the leading edge 
of a first-stage, inlet turbine vane. Consequently, a film cooling con- 
figuration that can provide both uniform and effective film coverage 
would be required in this critical area. Conventional injection geomet- 
ries, consisting of a row or multiple rows of holes angled to the surface 
in a streamwise direction, as depicted in Fig. 2(a), are frequently 
employed for turbine vane film cooling. However, the use of film cooling 
in the leading edge region is restricted by the surface which limits the 
minimum injection angle that can be used for coolant injection in the 
streamwise direction. As indicated in Fig. 2(a), for streamwise injec- 
tion through a typical vane wall thickness of 1.27*10 m (0.050 inch) and 
leading edge radius of 5.08x10 m (0.20 inch), the shallowest possible 
injection angle would be near 41° with respect to the surface tangent. 
Steep injection angles enhance mixing between the coolant and mainstream 
and should be avoided. 

An alternative film coolant injection configuration for the leading 
edge region is the use of coolant holes that are normal to the mainflow, 
yet angled to the surface in the lateral or spanwise direction as shown 
in Fig. 2(b). It has been shown (2) that a single hole with this orien- 
tation provides a wider field of film coverage than streamwise injection. 
Thus, in view of the surface curvature problem for leading edge cooling, 
the use of spanwise film coolant injection may provide the designer 
greater freedom to utilize shallow injection angles for improved film 
cooling performance. 


/ 
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In the present study, experiments were conducted to determine the 
film cooling performance for a single row of spanwise angled holes 
located in the leading edge region of a model turbine vane. The heat 
transfer and flow environment of the prototype high temperature, high 
pressure first stage turbine vane were simulated using a 3X size model 
vane positioned in a two-dimensional channel located downstream from a 
hot-gas generator. The inlet Reynolds Number (Re) , Mach number (Ma), 
and freestream turbulence intensity (Tu) were chosen to match the inlet 
conditions of the prototype vane. The model vane was oriented within 
the channel to reproduce the local Mach number distribution (Ma jj ) and 
flow deflection of the prototype vane. The model vane test surfaces 
were water cooled (internal ly) , independent of the gas film cooling 
system, in order to simulate a large freestream gas-to-wall temperature 
ratio (T t /T w ) representative of the boundary layer density variations 
for the high temperature prototype vane. 

The film cooling experiments were conducted using a single row of 
spanwise angled holes (3 « 18°, 35°) , a hole spacing-to-diameter ratio 
(S-]/d 0 ) of 4.0, with the row of holes positioned at three different 
locations in the leading edge region. The injection of cool air (T^ - T^) 
resulted in a hot gas-to-coolant temperature ratio of approximately 2.0, 
Miniature heat flux sensors were used to measure the local heat flux both 
with and without film coolant injection and the film cooling performance 
was presented in terms of the local Stanton number ratio, STp^/ST . 


Symbols introduced in parentheses are defined in the Nomenclature. 
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The objectives of the experimental film cooling study were to: 

(1) measure the reduction in heat flux to the vane surface due to 
coolant injection from a single row of spanwise angled holes, 

(2) determine the influence of the coolant injection angle, 3, 

(3) determine the influence of the location of the row of holes 
in the leading edge region, and 

(4) determine, for the simulated turbine environment {T t /T« ~ 2.0), 
the film coolant flow conditions that result in the optimum 
film cooling performance. 

The nature of the film coolant injection geometry studied, the range of 
cooling parameters investigated, and the presentation of heat transfer 
data obtained under simulated gas turbine conditions are the features of 
this investigation that will contribute to gas turbine design capabili- 
ties. 
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II. LITERATURE SURVEY 
II. 1. Background Information 

Prior to a discussion of the literature concerning previous film 
cooling studies relevant to gas turbine cooling applications, a brief 
orientation will be presented. Familiarization with the various film 
cooling techniques in light of their utility for turbine applications, 
introduction to and definition of terms that characterize the "cooling 
effect" of injection, and a brief description of the physical signifi- 
cance of key cooling parameters are items that will be included to make 
the presentation of the literature review more meaningful. 

A variety of methods can be used to effectively film cool a surface. 
The cool, secondary fluid may enter the external boundary layer through 
slots, holes, or a combination of the two. Figure 3 shows the basic geo- 
metric configurations used in gas turbine cooling. Tangential slots and 
slots with holes are well adapted to combustor liner cooling. Use of 
these techniques is not suitable for smaller components, such as turbine 
blades and vanes, for fabrication, structural , and aerodynamic reasons. 
Rows of angled holes, staggered rows of holes and interrupted, tangential 
slots in the vano and blade trailing edge are generally the accepted 
configurations for turbine applications. Although holes provide s prac- 
tical answer for small component cooling, continuous, tangential slot 


Hot Gas 




Figure 3. Typical Turbine Blade Film Cooling Configuraticix 



injection has been demonstrated to show greater cooling potential in 
light of its minimal interaction with the hot primary mainstream. 

In order to illustrate the comparative capability of cooling with 
a slot and a row of holes, Fig. 4(a) is presented. Sketches of the two 
cooling techniques are followed by representative cooling results. The 
basic features of the cooling techniques common to both configurations 
are the wall temperature downstream of injection (T^), the coolant blow- 
ing ratio, M = (pVWtpV) . the coolant injection angle (a), a charac- 
teristic injection dimension (slot width, $, or hole diameter, d Q )» and 
the downstream distance (x) from the point of injection. 

For an insulated wall, the downstream wall temperature (T adia- 
batic wall temperature) represents the degree to which the hot gas (T T ) 

OD 

has been reduced in temperature at x due to injection of cool air (T^). 
This effective reduction in Ty can be conveniently expressed as 

co 

(Ty - t adw)/(T t ~ T c ) and is commonly given the name, adiabatic wall 
cooling effectiveness 

The results presented in Fig. 4(a) show the variation in n^j with 
the dimensionless distance from injection (x-x-/d 0 ) for the two cooling 
techniques. A comparison is made between a single row of holes, angled 
in the freestream direction at 35° to the surface, and a two-dimensional 
slot at 30° to the wall. The holes are at 3 diameter spacings. The data 
shown represents both the effectiveness along the centerline of a hole 
and that characteristic of a lateral average across the span of the 
entire row. The slot effectiveness is represented with an equation by 
Hartnett, Birkebak and Eckert (3) that correlated their data for 
(x ‘/MS) > 60, where x 1 = x - x,. Since the data for the row of holes 


Adiabatic Effectiveness 


MSS 


Slot Cooling 


f 

Cool Film 


' r 'C 



Slot Coolxng 
Correlation 
(3) a = 30° 


S l /d 0 “ 3 *° 

A z/& Q = 0.0 (2) 

O lateral ave 
u a 0.61 rad {35°} 


Dimensionless Distance From Injection, (x-x^J/d 


Figure 4 (a) 


.Comparison of Slot Cooling and a 
Single How of Coolant Holes. 
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represents an optimum condition for that geometric configuration with 
M = 0.5, the slot effectiveness was calculated using an H of unity (an 
optimum coolant mass flux condition for slot cooling). The slot height 
S was then adjusted in order to inject the same amount of coolant mass 
as in the row of holes. For the given hole spacing (S^/d Q = 3.0), S 
was calculated to be 0.131dg. 

The slot is shown to be 3 times more effective near injection 
(x7d 0 - 15) and is still 40% better than the row of holes at 40 diam- 
eters downstream of injection. This comparison is based on the lateral 
average of the effectiveness. Though the slot configuration appears to 
be more effective than a row of holes, the decay rate of the cooling 
effectiveness for each technique is too severe for use in cooling exten- 
sive surface areas. A superposition of the effects of single row, hole 
cooling, by using multiple rows, would seem to be a logical extension 
in the development of practical cooling configurations. 

As the number of rows of holes are increased over a specified sur- 
face area and the spacing between holes decreases, the closely-packed 
array approaches, in the limit, transpiration cooling (a cooling technique 
in which the coolant is injected uniformly over the entire surface through 
a porous media). Figure 4(b) shows a comparison both pictcrially and 
graphically between transpiration and multiple row (also termed for many 
rows, "full -coverage") film cooling. The cooling parameters are essen- 
tially identical to those previously discussed with the added feature 
that the coolant flow is uniformly distributed over a transpired surface 
while, for multiple row cooling, injection occurs at discrete locations. 

Figure 4(b) presents data from Choe, Kays and Moffat (4) for an 11 
row array of holes at 5 diameter spacings and angled perpendicular to the 


Local Stanton Number, ST, 


Full-’Coverage Film Cooling 



Transpiration Cooling 



1.5 2 3 4 5 6 7 B 9 30 6 1.5 

Local Reynolds Number, Re x 


Figure 4(b), Comparison of Transpiration and Full 
Coverage Film Cooling. 
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main flow. An exit coolant temperature equal to the external wall tem- 
perature was prescribed (6^ = T}» The local Stanton number (ST^) is 
shown as a function of the local Reynolds number (Re^) for a blowing 
ratio, M, of Q.53 (represented by symbols}. The local Stanton number 
variation with Re^, is also shown as solid lines for an uncooled boundary 
layer and a transpired boundary layer at two different values of the 
blowing ratio, H (M ** 0.0027 and 0.0055). 

The results from the first 4 rows of holes (represented by the 
initial 4 symbols above the M - 0 line) show the initial penalty that 
must be paid in the form of increased heat transfer just due to the 
magnitude of the turbulent mixing generated by injection. Between rows 
5 and 11 the heat transfer to the surface is reduced by nearly 40?£ over 
the uncooled wall (M = 0). Following the final injection row, the heat 
flux continues to decrease due to the termination of injection-related 
turbulence and the cool mass addition accumulated from the previous 11 
rows. This so called "full-coverage 11 film cooling is compared to the 
heat transfer conditions that would exist over the same test area if the 
boundary layer was fully transpired. Blowing ratios (synonotnous with 
mass flux ratio or M) of 0. 0027 and 0.0055, representing one tenth and 
one fifth, respectively, of the total 11 row coolant flow rate are shown 
in relation to the full -coverage case. For transpiration mass flow 
rates approaching the total, hole coolant flow, there is at least an 
order of magnitude difference between the film cooled Stanton number, 
that can presently be achieved with full -coverage configurations, and 
the fully transpired boundary layer. 

With film cooling from multiple hole configurations shown to be one 
practical means of achieving effective overall cooling, information 
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related to single row cooling under a variety of conditions is still 
necessary as a building block for predicting the performance of full- 
coverage configurations. Additional research to understand the physics 
of injection and to apply this knowledge to creative configuration con- 
cepts is necessary if the gap between practical cooling and transpiration 
cooling is to be closed. 

Having placed the application of film cooling techniques in perspec- 
tive, it is now appropriate to discuss in detail the manner in which the 
potential of various cooling configurations can be characterized. In 
addition, the relationship between the cooling potential and pertinent 
injection parameters will be noted. 

11-1*1, Concept of Cooling Effectiveness 
II. 1.1.1. The Adiabatic Wall Temperature Effectiveness 

The primary objective of film cooling research to develop the 
capability of predicting wall temperature distributions on a film cooled 
surface, Eckert {5) postulated that the heat flux to a surface under- 
going film cooling can be approximated by 

q FC = h 0 ^ T ADW,FC ~ T W,Fc) W 

where 

qpQ " wall heat flux with film cooling 
hg - heat transfer coefficient without film cooling 
t ADW FC “ Makati c wall temperature with film cooling 
T w - actual wall temperature with film cooling 
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The film cooling adiabatic wall temperature reflects the 

degree to which the external hot gas has been reduced in temperature 
near the wall due to addition of a cooler fluid. This reduction in the 
local hot gas temperature is expressed as the following ratio 

\ " T ADW 

^ADW = T t - T c ■ ^ 

CO 

and given the name adiabatic wall cooling effectiveness, Film 

cooling studies conducted over an adiabatic surface yield adiabatic wall 
temperature distributions that are functions of such parameters as the 
amount of mass injection, the angle of injection, freestream turbulence 
characteristics and the local nature of the boundary layer. The measured 
adiabatic wall temperature represents a characteristic temperature of 
the boundary layer at a particular location downstream from injection 
which the wall sees as the driving force for heat transfer. 

Therefore, if one chooses a particular cooling scheme, the adiabatic 
effectiveness can be expressed as a known function of the distance from 
injection. Knowing the adiabatic wall temperature distribution from the 
effectiveness data then allows Eqn. (1) to be used to calculate the 
amount of heat removal needed in a design to maintain a specified wall 
temperature. A note of caution is in order, however, concerning the use 
of Eqn. (1). It was assumed that the local heat transfer coefficient 
was unaffected by the addition of mass into the boundary layer. Use of 
h Q as the local heat transfer coefficient under film cooling conditions 
is only valid far from the injection site. As film coolant is introduced 
into the hot gas boundary layer, disturbances are created by the injection 
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process (i.e. highly energetic turbulent conditions are established) 
which enhance the transport of mass, momentum and energy. The region 
affected by these hydrodynamic "blowing effects' 1 usually extends to 
fifty coolant hole diameters downstream of injection, but may reach as 
far as SO to 100 diameters as shown in Fig. 5(a) and (b). Before Eqn. 

(1) can be applied near injection, the variation in the local heat 
transfer coefficient must be known. 

For slot-type cooling, Fig. 5(a) shows the variation in the ratio 
of the local heat transfer coefficient with injection (h 1 ) to that with- 
out injection (h 0 ) for up to 90 slot widths downstream (6). The slot 
was angled at 30° to the surface and the blowing ratio, M, varied from 
0.34 to 0.87. At M = 0.87, as much as a 40?i increase in the local con- 
vective heat transfer coefficient was reported. This increase depends 
on M and can persist up to 40 slot widths downstream. The increase in 
(h'/hg) is probably related to the turbulence generated by the inter- 
action of the coolant stream and the main flow at injection. 

Figure 5(b) shows a similar trend in the heat transfer coefficient 
for injection from a row of holes angled at 35° to the surface and at 3 
hole diameter spacings (7). The ratio of the laterally averaged (across 
the span of the holes) heat transfer coefficient with injection, IP , to 
that without injection, hg, is shown as a function of M for downstream 
distances of 6.7 to 71.8 hole diameters. On the average, Tf'/hg decreases 
from unity at M = 0 to a minimum of 0.92 at M = 0.5. For M > 0.5, the 
injection-related turbulence may be acting to increase IP/hg in a 
monotonic fashion out to M - 2.0. At this condition lp/hg increases to 
nearly 1.25 for regions near injection (xVdg < 10). For the lower 
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values' of blowing ratio, h '/h^ changes only slightly downstream; however, 
for M > 1.0 the difference between the heat transfer coefficient near 
injection and far downstream begins to become more significant. 

In summary, the adiabatic wall temperature approach to film cooling 
requires some basic qualifications. The use of Eq. (1) is valid only 
far from injection where hydrodynamic "blowing" effects have diminished. 
To use Eq. (1) near injection knowledge of the heat transfer coefficient, 
h'» is essential. An alternate approach used to characterize film cool- 
ing effectiveness is outlined in the following section. 

11,1.1.2. Cooling Effectiveness From Heat Transfer Data 

The cooling effectiveness of an injection configuration can be 
defined from direct heat transfer measurements. Under this approach, 


where 


^FC = h FC (T T “ T W,Fc) 

CO * 


hp C - local heat transfer coefficient with film cooling 


( 3 ) 


A cooling effectiveness can be defined from Eq. (3) to reflect the 
reduction in heat flux due to film cooling. Thus, 



h re CT T„ " T W.FC ] 

h 0 CTf. - 'W 


( 4 ) 


whore 

qJJ c heat flux to wall without film cooling 
T w q - wall temperature without film cooling 
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When p C - Q , Eq. (4) is given the name, isothermal cooling 
effectiveness, Hjjqj where 


^ISO 


’o " q FC _ , h FC 

— 5o h o 


(5) 


This method of characterizing the cooling effectiveness incorporates 
into hp C all of the phenomena of coolant injection; namely, (a) the heat- 
sinking effect responsible for lowering the hot gas temperature near the 
wall, and (b) the hydrodynamic effects of injection. Mare simply, 


hpQ s hp C (Tq, hydrodynamic injection effects) (6) 

In comparing the adiabatic wall effectiveness concept and the iso- 
thermal effectiveness concept one finds that both methods require an 
investigation of two unknowns; namely, (a) the determination of the 
coolant temperature effect on either or (q^/qjj) and (b) the iden- 
tification of the hydrodynamic ‘'blowing effects" as manifest in either 
F'/hg or (ppc/qg) Tq = Tf . The isothermal effectiveness approach 
may prove mare convenient in that heat transfer measurements are made to 
identify both items (a) and (b) above. The adiabatic wall effectiveness 
method first requires insulated wall temperature measurements and then a 
separate series of heat transfer measurements to identify F'/hg. 

The film cooling effectiveness, whether adiabatic or isothermal, will 
be a function of numerous cooling parameters. The following discussion 
is intended to introduce these pertinent parameters in such a way as to 
build a foundation for a basic understanding of the phenomenon of gas 
film cooling. 
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II. 1.2. Introduction to Parameters Affecting Cooling Effectiveness 
The purpose of introducing the present section is to give the reader 
a general picture of the physical significance of some of the more impor- 
tant parameters that affect film cooling effectiveness. Documentation of 
statements made concerning the importance of each item will appear in the 
next section. The following discussion is intended only as a means to 
introduce the reader to the items which later will be presented in detail. 
Application of film cooling techniques to protect critical components 
requires a detailed understanding of numerous hydrodynamic, geometric and 
thermodynamic variables that are characteristic of discrete, secondary, 
coolant injection. The dominant factors affecting the degree of inter- 
action between a jet of cool injected fluid and an external, hot primary 
stream are: the coolant-freestream momentum flux ratio (or coolant-free- 
stream mass flux ratio), the injection angle of the coolant jet relative 
to the surface and the coolant-freestream temperature ratio. 

The trajectory of a jet injected into a hot gas mainstream depends 
to a great degree upon the momentum flux ratio, I, where 



P r V? v r 
-£4= (vA) M 

pji v “ 


( 7 ) 


The greater the coolant jet momentum, the deeper the penetration of the 
jet into the mainstream. Furthermore, as noted in Eq. (7), increases in 
I will mean subsequent increases in the amount of coolant mass added to 
the boundary layer (i.e. mass flow of coolant is proportional to M). 
However a condition may be established at sufficiently large values of I 
in which the coolant mass that is added is no longer in the proximity of 
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the wall due to the extreme penetration. The resulting jet trajectory 
is conducive to enhanced mixing with the hGi mainstream* realizing neg- 
ligible cooling effects at the wall. 

The trajectory of the coolant fluid is also related to both the 
degree of mainstream mass entrained into the injected jet and the drag 
forces acting on the exposed frontal area of the issuing jet. In the 
near regions of injection, the development of the coolant jet trajectory 
is dominated by any initial component of momentum parallel to the main- 
stream that may be imparted by angled injection. The combination of 
coolant momentum flux, mass entrainment and drag forces determines the 
extent to which the injected film covers the surface. In addition to 
keeping the injected jet as near the wall as practical, the coolant jet 
must function in a thermodynamic sense. The coolant temperature relative 
to the freestream level must provide the necessary heat sink capacity to 
aid in reducing the hot gas temperature near the wall. Variations in 
the coolant temperature must, however, be considered in light of its 
effect on the momentum flux ratio. 

Covering the surface to be protected with a continuous, uniform, 
cool film is the ultimate objective with discrete (hole) injection. The 
two most significant parameters associated with uniform film coverage 
are the center-to-center, coolant-hole spacing and the number of rows of 
coolant holes employed. The larger the spacing ratio, the greater the 
opportunity for the mainstream, hot fluid to freely move between and 
under the discrete coolant jets. This promotes greater interaction 
between the hot primary and cold secondary fluids. Reducing the distance 
between adjacent coolant jets creates a greater barrier to the oncoming 
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mainstream, resulting in jet trajectories that are less penetrating and, 
at the same time, reducing the effective interaction zone between coolant 
jets. Adding additional rows of holes with staggered spacing configure" 
tions allows consecutive downstream cooling holes to carry on when up- 
stream coolant injection has lost its cooling capacity. Multiple, 
staggered rov/s of holes approach the highly effective transpiration 
cooled configuration. Limitations on the amount of coolant fluid avail- 
able would weigh heavily upon the decision to incorporate multiple-row 
or full -coverage film cooling. 

The details of the coolant hole itself are coupled with the overall 
cooling performance of the injection jet. In order that the injected 
coolant mass leave the surface in some predetermined manner (i.e. injec- 
tion angle relative to the surface) care must be exercised in choosing 
the hole size in relation to the surface material thickness. The hole 
length-to-diameter ratio (L'/dg) must be of the order of 2 or greater to 
insure that the issuing jet does not leave the surface in a manner char- 
acteristic of thin-plate orifice flow. Small L'/dg ratios promote rapid 
mixing of the primary and secondary fluids, diffusing the coolant momen- 
tum equally in all directions and reducing the importance of angled in- 
jection to a second order effect. 

When modeling a cooling configuration, the hole diameter in relation 
to the boundary layer characteristics is significant. As mentioned ear- 
lier, the coolant jet trajectory is influenced by drag forces created by 
the flow of the deflecting mainstream over and around the exposed frontal 
area of the cylinder-shaped core of the issuing jet. Equivalence of the 
momentum boundav’y layer thickness-to-coolant hole diameter ratio (o./d-j) 
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in both a simulated and actual case insures that pressure forces or form 
drag terms will produce similar effects on the coolant jet trajectory. 

The nature of the external flow field may have a significant effect 
on the manner in which a coolant film performs following injection. The 
existence of large scale mainstream turbulence of sufficient energy may 
effectively act as an additional mixing agent between the coolant and 
mainstream fluids. Freestream turbulence intensities near 20 % (repre- 
sentative values at the exit of typical combustors) and the presence of 
turbulent "eddies' 1 with scales of the order of the film-cooled boundary 
layer thickness can contribute significantly to the deterioration of any 
injected film. In addition to freestream turbulence effects, injecting 
into boundary layers of varying thicknesses will result in a greater 
deflection of the issuing jet toward the surface for a thinner boundary 
layer than for a thick one at the same freestream velocity. Finally, 
wall curvature may cause the film cooled boundary layer to be more stable 
or unstable, depending upon the nature of the curvature. Instability of 
the film coolant layer will lead to separation of the film from the sur- 
face, thereby promoting mixing with the primary mainstream. 

The previous few pages have dealt with establishing a basis for 
understanding the physics of film injection. Those parameters considered 
important in determining the cooling potential of discrete secondary 
coolant injection have been presented and discussed in light of their 
physical significance. Documentation of the importance of such param- 
eters can be found in a vast library of film cooling studies carried out 
over the past 30 years. A sunwary of the most significant investigations 
dealing with the effects of the major film injection parameters will now 
be presented. 
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II. 2. Literature Review of Important Cooling Parameters 

II. 2.1. Effect of Coolant Injection Angle 
For a given cooling scheme, it is the penetration of the coolant jet 
into the mainstream and consequent mixing that primarily causes the cool- 
ant to lose its "identity" as a protective, heat-sinking mechanism. This 
penetration depends on the momentum of the coolant relative to that of 
the mainflow as well as the manner in which this momentum is directed (as 
accomplished by the injection angle). The coolant is injected through 
the surface in a direction that can be described by two basic angles, a-j, 
and ct 2 as shown by the sketch of the coolant coordinate system in Fig. 
6(a). Film coolant may be injected through the surface along the stream- 
wise direction with some streamwise angle a-j, maintaining ctg - 0°. This 
is represented by Fig. 6(b). For s tr oamwi s e injection (a-j, o:g - 0°), 
ctg is generally assumed to be zero and the streamwise configuration is 
described by the single angle, a (subscript dropped). Alternately, the 
film coolant may be injected in a spanwise or lateral direction (i.e. 
perpendicular to the mainstream flow) with at some value relative to 
the surface and a 2 » 90°. This is represented by Fig. 6(c). For lateral 
injection (ou, a? = 90°), cig is assumed to be 90° and the lateral con- 
figuration is described by a sinctle angle, 6 . In general there could be 
some compound angle of injection (a-j, a g) with the film coolant having 
components in both the streamwise and spanwise directions. 
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Figure 6. Definitions of Coolant Injection Angles. 


The great majority of film cooling data in the literature deals 
only with streamwise (a) injection, since that configuration imparts a 
component of coolant momentum parallel to the mainstream, thus improving 
the cooling effect downstream from injection. Interest in spanwise injec- 
tion (3) has stemmed from practical limitations, of fabricating streamwise 
angled holes on curved surfaces, and from the possible enhancement of 
film coolant uniformity in the lateral (i.e. spanwise) direction. Rela- 
tively little film cooling data are available for spanwise injection, 
especially from a row of angled holes under conditions of high values of 
the local acceleration. 

The effect of changing the injection angle has been well documented 
for 2-D slot injection and for streamwise injection from a single row of 
holes. The effects of changing the spanwise injection angle have been 
investigated for a single hole; however, no information is available for 
a row of spanwise angled holes. Papell (8) and Metzger et al. (9,10,11) 
have demonstrated the coolant may penetrate farther into the boundary 
layer as the angle of injection is increased from a = 0° (tangential slot) 
to a = 90°. Figure 7(a) is from Papell and shows the variation in 
with increasing injection angle for M of 0.5 to 4.5. The data is repre- 
sentative of one location downstream from injection corresponding to 26 
slot heights (x l /S - 26). For tangential slot injection (a = 0°) 

Is the highest for all M. As a increases n ADW decreases. For M = 4.5, 
the decrease in effectiveness is more severe with increasing a, possibly 
due to excessive coolant penetration into the boundary layer at such 
large values of M. 

Figure 7(b) describes the benefit of decreasing a from 90° to 35° 
for a single hole (2, 12). The cooling effectiveness is shown as a 
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function of the distance from injection for M = 0.5, The data is repre- 
sentative of that along the coolant hole centerline extended in the 
streamwise direction. The frees tream Reynolds number at injection based 
on the coolant hole diameter was 0.44 x 10 . The shallower angle of 35° 
shows a cooling effectiveness nearly twice that of a - 90° for all the 
downstream distances investigated. It is expected that the same trend 
of improved cooling effectiveness with decreasing injection angles would 
apply in a comparison of a row of holes. 

For spanwise injection, the effects of the coolant angle for a single 
hole are presented in Fig. 7(c) from Goldstein et al. (2), The figure 
shows contours of r^gy as a function of lateral distance (z/dg) and 
streamwise distance (x/dg) from the hole. Three injection angles are 
compared at M = 0.5: a = 35°, 3 = 35° and 3 = 15°. Compared to a = 35° 
the results for 3 - 35° show a somewhat wider field of film coverage; 
however, the maximum cooling effectiveness decreased from 0.4 to 0.3, 
Further changes in 3 from 35° to 15° show that the film coverage is 
extended even more. The highest measured is still lower than for 
a = 35° at the same M of 0.5. This evidence of extended film coverage 
at shallow lateral angles would seem to imply that a row of spanwise- 
angled holes might be used to achieve more uniform film cooling downstream 
from injection than streamwise angled holes (at same spacing ratio). The 
cooling effectiveness for a row of spanwise-angled holes is expected to 
show even less spanwise variation than shown for the single spanwise- 
angled hole, with the cooling effects of neighboring coolant jets super- 
imposing between holes. 


DIMENSIONLESS DISTANCE DOWNSTREAM 


Figure 7 (c) . Lines of Constant Cooling Effectiveness For 
Streamv/iso and Lateral Injection Angles (2) 
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11,2.2. Effect of Blowing Ratio 

In general, the coolant jet trajectory (i.e, penetration into the 
external boundary layer) is a function of the momentum flux ratio, I. 
However, for most cooling studies conducted with negligible density varia- 
tions between the coolant and mainstream, the blowing ratio (mass flux 
ratio) M has served as the major correlating parameter. M is defined 


M = 


p c v c 


( 8 ) 


It should be noted that for p c /p ot ~ 1 , M = /T and, thus, not only serves 
as an indication of the relative mass addition but also reflects the 
relative momentum fluxes. 

For 2-D, tangential slot cooling with the coolant gas temperature 
near the freestream level (pQ - p^), values of M near 1.0 have been shown 
to yield the maximum cooling effectiveness (13, 14). Figure 3(a) shows 
the variation of the mass transfer equivalent of namely n^» with 

the blowing ratio for slot injection (n^ is the impervious wall effec- 
tiveness with concentration levels analogous to the temperature levels 
of h^ DW )- Results for two downstream locations of 50 and 100 slot 
heights are presented. Near M = 1.0 the effectiveness was a maximum for 
both locations. As K increases from unity the velocity ratio must 
increase (p c « p w ) , yielding lower values. For M > 1.0 the relative 
differences in velocity between the coolant and freestream approach a 
condition in which mixing at the interface between the separate films 
increases. 

It is possible, however, to increase M well past unity and have the 
effectiveness improve steadily. This is accomplished for a tangential 
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slot by holding the ratio of velocities near unity and allowing the cool- 
ant temperature (density) to change appropriately. This permits minimum 
interface mixing* due to the velocity ratio fixed at unity, yet allows 
the momentum of the coolant along the surface to increase. Burns and 
Stollery (15), using a mass transfer analogy for measuring a concentra- 
tion cooling effectiveness, n^, as shown in Fig. 8(b), have demonstrated 
the gain in effectiveness by using a denser or cooler film with a veloc- 
ity ratio, Vq/V^ - 1.0. The effectiveness, rjjys increases at a fixed 
distance from injection (x'/S) for all increases in Pq/p^,* It appears 
from the evidence of Fig. 8(a) and Fig. 8(b) that the velocity ratio, 

V c /V ro , may be a significant parameter in determining the overall cooling 
effect of increasing M (i.e, as M increases ri^py may decrease if 
Vq/V^ >1). For 2-D slot inj’ection, V C /V M near unity results in the 
maximum effectiveness achievable for all P c /p to . 

Extension of the results pertaining to the blowing ratio effect on 
cooling from a 2-D slot to 3-D hole cooling is not possible. Coolant 
injected from a row of angled holes into a deflecting mainstream exhibits 
distinctively different behavior than that observed for the two-dimensional 
slot. Due to the ability of the mainstream gas to move between adjacent 
coolant jets and continually deform them downstream, an optimum blowing 
ratio for three-dimensional hole cooling can be expected to be quite dif- 
ferent from the two-dimensional case. Liess and Carnel (16), Eckert et 
al . (17., 18), Goldstein et al. (2, 19), Eriksen (7, 20, 21), Amana et al. 
(22), and Louis et al. (23), have shown that an optimum cooling effective- 
ness from a row of holes is achieved at M = 0.5 for streamwise injection 
angles below 35°. 


Impervious Wall 
Effectiveness, n 



Dimensionless Distance From Slot, x/Q 


Figure 8 (b) . Effect of Density Ratio on Tangential 
Slot, Cooling Effectiveness (15) . 
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Figure 8(c) shows the film cooling effectiveness as a function of M 
for a single hole and a row of holes both angled streamwise at 35° (2). 

The data is representative of 5 to 81 hole diameters downstream along a 
centerline through a coolant hole. Up to 11 hole diameters from injec- 
tion the optimum M for highest n ^ occurs at M = 0*5. For M > 0.5 in 
this region, n ADW decreases due to greater jet penetration into the main- 
flow. However for distances from injection greater than 11 hole diameters, 

reaches another maximum following M = 0.5. At M near 2,0, the 
adiabatic effectiveness at x'/dp = 80.67 is nearly the same as at M - 0.5 
for x'/dp = 31,47. This improvement attributed to the 

spreading of the coolant jets, resulting in a cooling reinforcement as 
they begin to coalesce. 

It is interesting to note that the optimum M of 0.5 actually repre- 
sents a velocity ratio, V c /V ra = 0.5 (Pq/p m * 1 for the aforementioned 
study). The significance of this value was pointed out by Pederson (24). 
Figure 9 represents a portion of the results presented by Pederson in 
which the coolant density was changed from Pq/p^ of 0.75 to 4.17. Figure 
9 shows a spanwise averaged n T ^ (based on concentration measurements using 
a mass transfer analogy) as a function of M for a row of 35° streamwise 
angled holes at 3 diameter spacings. The data is for x'/dp = 10,29. 
Examination of Fig. 9 will reveal that the maximum effectiveness for 
each Pq/Pco occurs near a velocity ratio between 0.4 and 0.5. Apparently, 
an equilibrium condition is established at V c /V ro = 0.5 between disturbances 
created by injection (i.e. local turbulence due to mixing brought on by 
increased coolant jet velocities) and additional heat sinking capacity 
due to increased cool mass addition. 
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The effect of M an for spanwise injection is not well docu- 
mented. Goldstein et al. (2) have found that for lateral coolant injec- 
tion from a single hole, the field of film coverage is slightly reduced 
in going from M = 0,5 to 1.0, Figures 10(a) and (b) are effectiveness 
contours as a function of lateral distance (z/d Q ) and streamwlse dis- 
tance (x Vdg) for a - 35°, 3 - 35° and 3 = 15°. There does not appear 
to be a significant change in r^j levels in going from M = 0.5 to M - 
1.0 for 3 = 35° or 3 = 15°. At M = 1.0 the lateral coverage is reduced 
from M = 0.5, No definite conclusions can be drawn from this single 
hole data in order to predict the overall effect of M on the cooling 
effectiveness for a row of spanwise angled holes. 


II. 2.3, Effect of Injected Coolant Temperature 
The coolant temperature is the driving mechanism responsible for 
lowering the effective boundary- layer gas temperature. Not only does 
lowering the coolant temperature provide an additional heat-sinking capa- 
city to the injected film, but it also yields a lower trajectory of the 
injection jet for constant M (i.e, constant coolant mass addition). 
Referring again to Fig. 9, dealing with the work of Pederson, the 
effect on n IW of increasing the coolant density (i.e. decreasing the 
coolant temperature) from Pq/p^ " 0.75 to 4.17 is presented. At 
x'/dg = 10.29, n IW increases at any M for increasing Pq/p^* It is not, 
however, the increase in p r (i.e. decrease in T^.) which affects the 
improvement in the cooling effectiveness. Because increases at a 
constant coolant mass addition (i.e. M a constant), the velocity ratio 
must decrease. This means at constant M, the highest coolant density 
dictates a velocity ratio which is characteristic of a minimum disturbance 
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to the boundary layer. Consequently, the effective transport mechanism 
between the hot gas and injected, cool film is greatly reduced (i.e. the 
effectiveness is greater). 

Additionally, lowering the coolant temperature (i.e. increasing 
Pq/p^) for constant M results in a more shallow coolant jet trajectory. 
Recall that 

V r 

I - M (vp) (9) 

v co 

and that for constant M and decreasing T c , V^/V ro must decrease. Con- 
sequently, I must decrease with a concomitant jet trajectory that keeps 
the coolant fluid nearer to the wall. 

11,2.4. Injection Geometry Considerations 

In addition to the coolant injection angle, there are a few other 
basic geometric factors that are intrinsic to a particular cooling scheme. 
Such items as the spacing between coolant holes, use of multiple rows of 
coolant holes in various configurations of staggered or in-line arrange- 
ments, hole length- to-diameter ratio and the thickness of the boundary 
layer at the point of injection are all necessary considerations when 
attempting to predict the overall cooling capability of coolant injection 
from a row of holes. 

The effects of the coolant hole spacing on the cooling effective- 
ness are manifest in the degree to which the hot mainstream can freely 
move between and interact with the individual coolant jets. Li ess and 
Carnel £16) investigated hole spacing-to-diameter ratios (S-j/dg) of 2.22, 
3.33 and 4,0. Their results appear in Fig. 11. The adiabatic cooling 
effectiveness is shown as a funr+ion of the dimensionless downstream 
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distance (x'/dg) for a blowing ratio, M ~ 0.85. Spanwise variations of 
n ADW are stlown for eac ^ S-j/dg W1 ^ va ^ UGS 2 / s -r (dimensionless spanwise 
distance between two holes) ranging from 0.5 (midway between adjacent 
holes) and 0 (in-line downstream with a hole). 

Very near injection (x'/dg ~ 1.0), closer spacing of the holes 
(S-j/dg = 2.22) yields higher film effectiveness than for S^/d Q of 3.33 
and 4.0. Further downstream, for x'/dg in the range of 2 to 10, n ADW 
along the jet centerline temporarily increases as S-j/dg increases due to 
the three-dimensional effect of jet spreading and subsequent coalescence. 
However, for x'/dg >10, lower values of the spacing ratio yield higher 
effectiveness for all remaining stations downstream as well as all 
lateral positions across the span of the test surface. A uniform, two- 
dimensional film is a desirable characteristic for any application of 
film coolii.g. As can be seen from Fig. 11 for M = 0.85, the data for 
the lowest value of the spacing ratio show essentially no lateral vari- 
ation for distances beyond 13 hole diameters downstream. In contrast, 
the data for the larger values of the spacing ratio show significant 
lateral non-uniformity for distances up to 25 hr<le diameters. Hess and 
Camel also noted that larger downstream distances are required for uniform 
spanwise film effectiveness if the blowing ratio is decreased from M = 

0.85. 

Practically speaking, the film effectiveness from a single row of 
holes falls to such a low value at such a rapid rate that multiple row 
cooling configurations may be necessary to cool any appreciable surface 
area. Arrangements of multiple row configurations also can compensate 
for the spanwise film non- uniformity from a single row, as staggered hole 
arrays fill in the gaps between holes downstream. 
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Louis et al. (23) demonstrated the effect of adding just an addition- 
al row of holes incorporating both 10° and 20° injection. The injection 
of additional mass and the modification of the effective decay charac- 
teristics of the film showed the double row of holes to be nearly as 
effective as a slot. Hiroki and Katsumata (25) measured approximately a 
20% increase in the local film effectiveness by increasing a two row 
array to a four row configuration for cooling rates of 1% to 1 ,5" of the 
external gas flow. 

Larger arrays of holes have been studied by Choe, Kays and Moffat 
(4). They used 11 rows of film cooling holes at hole array spacings of 
5 and 10 hole diameters. The overall effect of the cooling scheme was 
defined as a reduction in the local measured Stanton number due to cooling 
over that for the uncooled condition. For the 10 diameter spacing con- 
figuration minimal decreases in the local Stanton number were recorded 
in comparison to the uncooled wall. However, reducing the spacing to 5 
diameters caused the Stanton number to fall by 50% over the 10 diameter 
geometry. 

Metzger, Takeuchi and Kuenstler (26) used 10 rows of 90° injection 
holes at a 4.8:1 pitch- to -diameter ratio to study the effects of stag- 
gered and in-line multiple row arrays. Rows of holes with all holes 
arranged in-line downstream proved inferior to staggered arrangements. 
In-line arrays prevented upstream cooling jets from filling in the 
spaces between the cooling holes on the next downstream row. Although 
multiple, staggered rows of holes (full coverage cooling) would seem 
necessary for cooling an extensive area, it should be noted that regions 
exist, such as in the turbine vane leading edge, in which high heat 
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fluxes are concentrated over short distances. Consequently, single row 
cooling configurations capable of effectively protecting such critical 
regions need to be investigated. Results from single row studies, as 
will be presented for this investigation, may also serve as "building 
blocks" to predict the overall cooling effect of a superposition of many 
similar rows of holes. 

The injection scheme adopted in light of all its geometric, hydro- 
dynamic and thermodynamic characteristics must also be viewed in relation 
to the nature of the local boundary layer at injection. Thicker momentum 
boundary layers allow the coolant injection jets to penetrate farther 
into the external mainstream for identical relative mass addition and 
momentian flux. Eriksen and Goldstein (7) and Goldstein, Eckert, Eriksen 
and Ramsey (2) have shown that increasing the boundary layer thickness at 
injection with constant mass addition leads to a reduction in the film 
effectiveness for a single hole. However, the effect of injecting at 
thicker boundary layer conditions is less severe for a row of holes, 
showing even a slight increase in the effectiveness for a 50 % increase 
in displacement thickness-to-hole diameter ratio. Li ess (27) reports 
varying the displacement boundary layer thickness from S-j/dg = 0.1 to 0.6 
for a row of holes injected at 35° to the surface and in the streamwise 
direction. Figure 12 shows the variation in with displacement 

boundary layer thickness-to-hole diameter ratio at M = 0.3 to 0.43. The 
effectiveness variation with 6^/d Q is shown for three downstream locations 
of x*/d 0 = 10, 30 and 50. For the given M, was shown to decrease by 
as much as 50% from S-j/dg » 0,1 to 0.6. Figure 12 shows that the effec- 
tiveness decrease with increasing boundary layer thickness holds at least 
as far as 50 hole diameters downstream from injection. 


Dimensionless Boundary Layer Displacement 
Thickness, 

Figure 12, Variation in Cooling Effectiveness 
Wrth Increasing Boundary Layer 
Displacement Thickness (27} . 
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II. 2. 5. Effects of Freestream Conditions 
Though the momentum flux ratio, angle and point of injection, and 
configuration geometry may be viewed as primary considerations in eval- 
uating a film cooling scheme, a number of other parameters remain that 
may collectively result in a considerable perturbation to the expected 
film coolant performance. 

For instance, the transport of heat and momentum are greatly enhanced 
in the presence of a mainstream containing large fluctuating velocity com- 
ponents. High levels of freestream turbulent kinetic energy (as manifest 
in the turbulence intensity, Tu ra ) will lead to enhanced mixing between 
the coolant jets and mainstream provided the external turbulence has 
energy of the same order as that created by coolant injection into the 
boundary layer. Such a condition may exist in the leading edge region 
of a turbine airfoil. Sutera, Header and Kestin (28) and Sahed et al. 
(29,30) have shown that, if the oncoming turbulence scale of the free- 
stream is larger than approximately three times the Hiemenz boundary 
layer thickness, considerable amplification of the turbulent fluctuations 
takes place in the stagnation region. Figure 13(a) demonstrates that 
at £ = 12.12, where 

a 1/Z 

% = (f) x (10) 

where a = 4V to /D 

V w - freestream velocity 
D = characteristic leading edge dimension 
v = kinematic viscosity 
x = coordinate along surface 
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the intensity of the fluctuating velocity component in the direction of 
flow reaches up to 45% based on the freestream velocity at the outer edge 
of the boundary layer. Traci and Wilcox (31) present a stagnation region 
turbulence amplification model that also predicts dramatic effects. 
Depending on the freestream Reynolds numbers the turbulent energy ap- 
proaching the stagnation region can increase from 30% to 300%. 

Ramsey and Goldstein (32) report turbulence intensities for 90® in- 
jection at blowing ratios of 1.5 and 2.0, Figure 13(b) shows that maximum 
values in ti. - near injection region are of the same order as intensities 
generated in the stagnation region due to freestream vortex filament 
stretching. At lower blowing ratios associated with optimum effective- 
ness levels, the injection generated turbulence should diminish, increas- 
ing the dependency of the cooling effectiveness on the relatively high 
incident freestream turbulence levels. 

Far removed from stagnation-region amplification effects, freestream 
turbulence levels should not significantly enhance the transport mechanisms 
already dominated by injection effects. Nominal turbulent boundary layer 
intensities are significantly lower than those turbulent conditions created 
by discrete coolant injection. Junkhan and Serovy (33) and Kearney et al. 
(34), under zero pressure gradient and strongly accelerated conditions 
respectively, show little or no effect on the local heat transfer through 
a turbulent boundary layer exposed to freestream intensities up to 10%, 

The three-dimensional vortex stretching phenomenon near a stagnation region 
appears to be the only freestream, turbulence-related condition that could 
possibly influence the mixing behavior characteristic of the highly 
energetic regions of the turbulent injection jets. 
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The effects of freestream acceleration due to large imposed pressure 
gradients on turbulent boundary layer heat transfer without film coolant 
injection have been well documented, Kays et al . (35, 36, 37, 38, 39) 
show a decrease in the local heat transfer for moderate values of the 
acceleration parameter, K, where 


where 


_ v 


CO 


dv 

o 

dx" 


v = kinematic viscosity 
a local velocity 
x = coordinate along surface 


( 11 ) 


Apparently, with increasing levels of K the boundary layer turbulence is 

suppressed, reaching a lower average turbulence intensity than exists 

without acceleration, yet still retaining a turbulent layer appearance. 

At K > 3.3 x 10 a steep drop in the Stanton number occurs, approaching 

a "laminarization" condition in which the heat transfer nearly approaches 

values characteristic of a laminar boundary layer. 

As applied to film cooling applications, the effects of increased K 

on the overall heat transfer to a cooled surface were noted by Kays et al, 

(30). Depending upon the blowing ratio, increasing K can either increase 

or decrease the local heat transfer. However, for high transpiration 

blowing ratios (M ~ 0.0062), the Stanton number increased as much as 25$ 

-fi 

in varying K from 0 to 0.77 x 10 , For film injection from discrete 
locations (i.e. holes) much larger blowing ratios are generally maintained 
than those noted in the transpiration study. Consequently, the Stanton 
number should rise even higher than noted for the transpired condition 
of M = 0.0062. Heat transfer measurements of a film cooled, accelerated 



rUIEI 


51 

boundary layer reflect both the boundary layer structure changes due to 
blowing as well as any acceleration induced mixing that may occur in the 
outer edges of the boundary layer due to large local velocity gradients. 

Adiabatic wall temperature measurements for film coolant injection 
under varying degrees of freestream acceleration would not reflect 
acceleration-induced boundary layer perturbations (as do heat transfer 
measurements) but would give an indication of the degree of mixing which 
takes place between the coolant and mainstream due to imposed velocity 
gradients at injection. Changes in the adiabatic wall temperature 
(reflecting the degree to which the hot gas mainstream is reduced in 
temperature near the wall) due to acceleration were reported by Carlson 
and Talmor (40). They measured acceleration effects on in the 
presence of freestream turbulence for 2-D slot film cooling. They found 
that favorable pressure gradients lowered the cooling effectiveness, 
except for the case in which high freestream turbulence intensity levels 
existed at injection. Initially the favorable pressure gradient tended 
to suppress the high turbulence level, and, hence, reduce the effective 
transport mechanism between the coolant and the mainstream. However, a 
condition was established once the turbulence was diminished in which the 
acceleration led to a deterioration of the film effectiveness. 

Acceleration affects on for three-dimensional hole injection 

were reported by Liess (27). At x’/dg of 10 and M - 0.5, was reduced 

nearly 50?i over the unaccelerated case for K near 0.4 x 10" . At M = 1.0, 

£ 

^ADW was rec * ucsc * 22SS for a variation in K from approximately 0.14 x 10 
to 0.4 x io“ 6 . 
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In addition to pressure gradient effects, surface curvature intro- 
duces additional complications into determining film cooling behavior in 
accelerating flows. Nicholas and LeMeur (41) attempted to separate curva- 
ture effects from pressure gradient effects. They found that for film 
cooling on a convex surface the effectiveness was greater than on a flat 
plate of zero pressure gradient, but slightly lower than on a flat surface 
with a pressure gradient similar to the convex case. However, their 
results indicating an increase in the effectiveness for Injection under a 
favorable pressure gradient over a flat plate condition are in disagree- 
ment with the majority of cooling studies conducted in the presence of 
freestream acceleration, For a concave wall, the effectiveness was 
generally lower than the zero pressure gradient results, except at very 
high blowing ratios (M - 2.0). The concave wall had a pressure distri- 
bution essentially equivalent to a flat plate. The effectiveness reduc- 
tion due to concave curvature was as much as 30% far downstream from 
injection and negligible near injection for a blowing ratio of 1.0. 

II. 2. 6. Models for Film Cooling Effectiveness 
Numerous models and correlations for both slot and hole cooling have 
been developed without achieving applicability under all injection con- 
ditions and configurations. Instead of basing predictions on overall 
injection similarity parameters, the present trend of film cooling re- 
search is to understand the micro-nature of injection. Only in estab- 
lishing the dependency of the mass, momentum and energy transport pro- 
cesses on injection conditions can a truly universal model be developed. 
Direct solution of the boundary layer equations, using various boundary 
conditions and transport models that are characteristic of individual 
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injection schemes, is necessary to achieve the degree of universality 
that is required in modern applications. 

Until knowledge of the injected boundary layer mi crostructure is 
complete, familiarity with existing correlations is essential from a 
design viewpoint. These correlations are almost entirely two-dimensional 
but can provide a basis for understanding the trends and importance of 
injection parameters common to both two and three-dimensional cooling. 
Appendix A summarizes the majority of film cooling models and correlations 
for slot and hole cooling with their representative assumptions and appli- 
cability. The purpose of Appendix A is twofoldr (1) to provide a 
foundation for the design of a film cooling configuration for a specific 
application under given constraints, and (2) to emphasize the need for 
a more fundamental approach to modeling film injection. 
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III. MODELING CRITERIA FOR GAS TURBINE FILM COOLING SIMULATION 

As emphasized in Chapter I, the subject investigation addresses it- 
self to the problem of cooling critical components of a gas turbine engine. 
Specifically, attention is focused on the high heat-flux, leading edge 
region of an inlet turbine vane. An injection geometry is proposed which 
is compatible with leading edge curvature constraints but must be quali- 
fied as to its film cooling effectiveness. Chapter II introduced the 
many pertinent parameters previously found to directly influence the cool- 
ing effectiveness, but offered no basic analytical justification of their 
relevance. The mathematical development of relationships which express 
the importance of characteristic similarity variables for coolant injec- 
tion is a necessary exercise if the subject experimental study is to 
provide film cooling information directly transferable to high- temper- 
ature, high-pressure gas turbine applications. 

The present Chapter will examine the features of gas film cooling 
using Fig, 14 as a basic model. A three region model is defined to char- 
acterize a film cooled surface: (1) the flow along the surface prior to 

injection, (2) the area of intensive coolant-jet and boundary layer inter- 
action, and (3) the regime downstream from injection within which the 
disturbed, blown boundary layer seeks to re-establish normal similarity 
both hydrodynamics'! ly and thermodynamically. By examining the governing 
equations characteristic of the three regions, the governing parameters 
are identified such that an experimental study can be conducted at 
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reduced flow conditions that closely approximates the heat transfer 
environment typical of high-pressure and high-temperature turbine oper- 
ation. 


Ill . 1 . Heat Transfer Modeling at Reduced Floy/ Conditions 
Equations (12) ana (13) are the non-dimensionalized boundary layer 
equations for x-momentum and energy (42) representing the initial region, 
prior to injection, of the three-region model. 


p , ( U , |;+v 


3U‘\ 

3y J 
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+ G p , 0 , CQS 1 + _L 
9x Re 2 Ke 3yf 2 
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ay 


3x 


Re 


$ 


(12) 

(13) 



G = Grashof number, p 2 g3 0 (T T - T) L 3 /p 2 

co 

Re = Reynolds number, pV^x/u 
E = Eckert number, vf/C (Tj - T^) 
i = angle between surface and direction of gradient of 
gravitational potential field 
$ = viscous dissipation function 


It will be shown that by reproducing the same local Mach number and 
Reynolds number distribution along a geometrically similar surface, the 
dimensionless heat transfer will be similar. A relation exists between 
the total pressure and total temperature levels of an actual turbine 
environment and those at reduced conditions that will insure similarity 
in the local heat transfer. 
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For the special case of flow over a turbine airfoil, the ratio of 
inertial forces to bouyant forces is large (Re » G), and Eqs. (12) and 
(13) reduce to 



Solution of Eqs, (14) and (15) implies the functional dependence of the 
velocity and temperature fields, 

77- 83 rr (x‘, y‘> z' 5 Re» Pr, E) (16) 

U oo 

T t - T (T t - T) 

0' = J - y 1 ** — jf (x* , y', z', Re, Pr, E) (17) 

T co ^ 

For conditions in which the Prandtl number does not vary signifi- 
cantly, and the kinetic energy available in the freestream for dissipation 
is small compared to the thermal energy, the heat transfer phenomena at 
reduced flow conditions will simulate those at elevated temperatures and 
pressures if equivalence in the Reynolds number is maintained. For the 
remainder of this chapter attention will be focused on maintaining 
similarity between actual turbine engine corJitions (E) and those at 
reduced temperature and pressure characteristic of turbine modeling (M). 
Consequently, for similarity in heat transfer 

Re^ = Re^ (18) 

where M s model 

E = engine 

provided similar geometries are considered. 
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By definition. 

Re = ^ 

(19) 


x u 


and 

« MaV* = Ma(Yg RT) 1/2 

V 


V 

(20) 


Using the compressible flow relations to define the freestream 
stagnation conditions external to the boundary layer. 


T = 
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and 
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( 21 ) 




( 22 ) 


RT t (1 + Ma ) 

The local Reynolds number can now be written as 
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Requiring equivalence in the local Reynolds numbers between the model 
and the engine results in 
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Provided the local Mach number can be simulated at geometrically 


similar locations 

such that 



< ■ < 

(25) 

and 


(26) 


then Eq. (24) can be rewritten using Eq. (25) and (26) as 
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where 


f i - 


y E - 2 (Y E +1)/2(V E -1) 

[1 + Lrl (Ha/) ] 

H + X -f 1 (Ma x H ) ] 


(28) 


and 


l m 

C tJ jp = geometric modeling factor = -g 

^ char 


I L 

char 5 characteristic dimension of engine component 

M 

^ char ^ c * ial, ' aci:er ^ s ^ 1 * c dimension of similar model component 

For large differences between engine and model conditions with 
Ma v ^ = Ma JS f, does not significantly deviate from unity. For examples 
for T t = 1922°K (3460°R) and an ASTM-A-1 fuel-to-air ratio of 0.03, 

r 00 M 

y L - 1.263. Using T T = 700°K (1260°R) and a natural gas fuel-to-air 

oo 

ratio of 0.01 (typical for reduced conditions simulating the engine 
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M 

environment), r = 1.358. Under these specifications fj ranges from 
0.9895 at Ma - 1.0 to 0.9998 for Ma = 0.10. Therefore Eq. (27) can be 

A X 

approximated as 
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(29) 


Equation (29) states that for modeling a given engine condition at a 
particular geometric sc^e factor, the reduced operating pressure is 
fixed once a total gas temperature is selected. It is recognized at 
this point that the development of the similarity criteria has neglected 
chemical modeling. For operation at high turbine inlet temperatures, 
the effects of the combustion processes may indeed be important. How- 
ever, at reduced temperatures and pressures, simulation of chemical 
phenomena will not be attempted herein. 


III. 2. Identification of Pertinent Film Cooling Parameters 
Simulation of the near injection region requires a consideration of 
the coolant-jet and hot gas boundary layer interaction. The degree of 
cooling achievable in this region is primarily dominated by the jet- 
trajectory. Immediately following injection, the coolant jets are acted 
upon inertially rather than in a viscous manner characteristic of turbu- 
lent mixing. The initial dynamic interaction of the boundary layer with 
the cylindrical profile of the coolant jet establishes the form drag 
forces that determine the jet trajectory. Keeping the jet axis near the 
wall, wholly within the boundary layer, is necessary to realize signifi- 
cant cooling effects. 
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In the following discussion, a simplified analysis of a jet in a 
crossflow is presented. This mathematical description of the trajectory 
of the jet, representative of injection from a film cooling hole, shows 
that the injection angle, coolant momentum, and boundary layer thickness 
at injection are the major parameters governing the jet trajectory. 

Figure 15 shows a control volume for the analysis of an issuing 
coolant jet. It is assumed that the jet momentum along the y axis is 
constant and equal to the initial value, 


^ [sin X pVjdA] = 0 


JET 

which can be integrated to give 


sin X J pVjdA = p c V 2 A Q sin X Q 
JET 

2 

where p^Vq = initial coolant momentum flux 
A q = area of coolant hole 
X Q = initial injection angle 

The x-momentum equation for the jet can be written as 


(30) 


(31) 


j f o vm ird n 

dx ^ cos x J pV J dA ^ " cos x = F x,drag t 32 ) 

JET 

H 

where M is the mass flux entrained from the freestream and F is a 

c A 

drag force. 

Assuming the entrainment is initially small, Eq. (32) can be 
approximated as 

37 [cos x J = F x,drag 

.ll-T 


( 33 ) 





63 


•-rtra'iHfiL 


EWKODTOi ilUT? Ot? 1 . 
nnMMlt-Si i *J'.h ■■ ^ ■'■'•''■ 




The cylinder-in-crossflow drag force can be estimated as follows 

F* = \ Pj&O tan X C D (34) 

0 

where p ra V^ = freestream momentum flux 
d Q - coolant hole diameter 
C Q - cylinder-in-crossflow drag coefficient 

Equation (34) represents the drag force per unit length on the jet core 
assuming no spreading and slug type flow. To account for the jet spread, 
Abramovich (43) suggests a linear relation 

d = d 0 (1 + a (35) 

where a = empirical spreading constant 

Also a correction factor, C*, for non-slug type flow of an actual 
boundary layer can be used to adjust the drag force term. Equation (34) 
can now be expressed as 


F x ■ ? PJ& D C*d 0 (1 + a tan X 
Thus Eq. (33) reduces to the form 


(36) 


dx 


[cos X 


pVjdfl] 


JET 


\ Pj& D C*d 0 (1 + a f) tan X 


(37) 


To simplify Eq. (37) further* the integral term may be expressed 


as 


cos X 


’ » v j dA * W s1n x 1 <*5* 

JET JET 


(38) 
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and by using Eq. (31) in the following form 

cos X | pVjdA = S1 ' n x o cot * 

JET 

Eqn. (37) can be rewritten as 

0C V C A o sin X o 3x (cot X) = \ P - V » C D C * d 0 (1 + a ^ tan X 


( 39 ) 


where 


=otX = f 


(40) 


(41) 


Upon substitution of Eq. (41), Eq. (40) becomes a second-order, ordinary 
differential equation that can be integrated twice to yield, 

ct*x' 2 + x' + /a *x ' + — 
i tan 

y ' = — hi [ 

1 , vfc*a 

tanS: 4 

o 

4C D C* 

where c* = - !si - -^ 

5 * Pc V M 

y' = y/d Q ; x' = x/d Q 

The ratio of the coolant momentum flux to the freestream momentum 
flux (I), together with the initial injection angle (X Q ) and character 
of the boundary layer at injection (C*) are shown in Eq. (42) to be 
important parameters in modeling the near injection regime. 

Some caution must be exercised in using this jet trajectory model 
to fully explain film cooling effectiveness trends with variations in 
certain injection parameters. Freestream mass entrainment was neglected 




r 4 


(42) 
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in the development of Eq. (42). “1 though Eq. (42) may provide a reason- 

able estimate of the general jet location with respect to the wall, the 
film cooling effectiveness is expected to be strongly influenced by a 
mixing term (as represented by the jet entrainment). Thus, while Eq. 

(42) provides a reasonable representation of some of the more important 
film injection hydrodynamic parameters, an evaluation of the cooling 
effectiveness would require a more accurate treatment of the important 
mixing phenomena. 

The final region of the film cooling model characterizes the mixing 
process between adjacent jets and the boundary layer. Up to this point, 
each jet development was analyzed individually and assumed to be indepen- 
dent of the proximity of neighboring jets (i.e. a row of holes). The 
third region of the film cooling model is now discussed in which it is 
assumed that all adjacent jets have coalesced and mixed to a uniform 
degree. Figure 16 is a sketch showing a control volume for analysis of 
the mixing between the coolant and the mainstream. The important coolant 
and freestream parameters that can be expected to govern the film cooling 
performance are deduced from a basic heat sink analysis following that 
presented by Goldstein (44). A characteristic temperature, representative 
of the average, mixed temperature of the film cooled boundary layer can 
be defined by an application of the conservation of mass and energy for 
a station, x, downstream from the near injection region. The mass flow 
in the boundary layer can be expressed as 


( 43 ) 
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where 


ROTOBOC-PJlLmr OF W7 
ORMili-J A j i :u, 

m c B film coolant mass flow injected 

m = boundary layer mass flow entrained from the freestream 


If the wall over which the fluid flows is adiabatic, the energy equation 
yields 


m C C pC T C * 


m m C n T = 

co n co 


™out C p T x 


(44) 


And combining Eqs. (43) and (44), assuming a constant specific heat 
(i.e. Cp C » sr cp, the average temperature characteristic of the 
completely mixed, film cooled boundary layer is obtained. 


T = 
x 


* M 4 __ 

m r T r + m T 

L L » co 

(n> c + "U 


(45) 


T approximates the film cooled adiabatic wall temperature which governs 

A 

the heat flux to the wall for the case of a cooled (non-adiabatic) 
surface (44), 

The mass flow of the mainstream into the boundary layer control 
volume, m^, can be approximated by assuming a 1/7 power velocity profile 
and p 5 ; p M . 


m - 

03 


r° r 1 v 1/7 v 
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where 

Equation (46) then reduces to 


L = span of cooled surface 
**P 


(46) 


where 


4. ■ 7 /8 P-Wsp 

5 = boundary layer thickness at x 
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(47) 
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The coolant flow rate can be written as 


Trd; 


m. 


0 


C ■ W *hol, * l»cV " -f 

where N = number of film coolant holes 

The average boundary layer temperature may then be written. 


( 48 ) 
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For L fs S-jN, Eq. (49) can be expressed as 
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(50) 


7(o x / d 0 )(S 1 /'d 0 ) 

Equation (50) has been developed simply to bring out the important 
features of coolant injection that are responsible for lowering the 
temperature of the mainstream boundary layer downstream of the coolant 
holes. Equation (50) indicates that as M -*■ 0 (i.e. no coolant injection 
into the hot gas boundary layer) T = T . Also, as the coolant film 

A w 

moves farther from the injection site, for any M, (6 X / d Q ) will increase 
resulting in diminished effects of injection and T approaching T^. 

Drawing on the information derived from the three- region film injec- 


tion model 3 the following guidelines for proper experimental modeling 
were observed: 
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(1) to maintain similar velocity and temperature fields at reduced 
flow conditions, a unique relation between the total gas tem- 
perature and pressure is required, 

(2) the momentum flux ratio, freestream mass entrain- 

IJ 

ment (M Q ), the injection angle (X.) and the state of the 
momentum boundary layer thickness at injection (C*) are 
significant parameters for the near-injection region, and 

(3) the film coolant blowing ratio (M), the hole spacing (S^/d Q ), 
the boundary layer thickness ratio (6 x /dg), as well as the 
coolant- to-gas temperature ratio (Tg/Tj determine the down- 
stream film cooling performance. 

To achieve the objective of determining the cooling effectiveness of a 
proposed turbine vane leading-edge injection configuration, it is impor- 
tant to properly model each of the items described in the aforementioned 
film cooling model such that results obtained at reduced flow conditions 
(model) are applicable to film cooling under realistic engine conditions. 

To maintain similarity between an engine cooling condition and a 
model simulation, for a typical hole spacing configuration and injection 
angle, equivalence of the mass and momentum flux ratios must be maintained. 


fiI M = ( p cV _ (p C V C )E 
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(52) 
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Equations (51) and (52) can only be satisfied simultaneously when 

V r M V r E 

(r> = (sr) 


CD 


(53) 


Furthermore, Eqs. (51) and (53) require that 


, P C 


M Q r E 




(54) 


Equation (54) can also be expressed as 


P r M P. E 

(r 'T") (r x ) 
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(55) 


Rearranging Eq. (55) , one finds that 


m R E P r E P M M 

< R » T / ■ <or> <r> (pr> ( R cV 

C C “ c 


(56) 


With the mass flux ratio and momentum flux ratios referenced to the exit 
conditions of the coolant holes. 


P r E P r M 

</> = (/-) = i 


(57) 


Eq. (56) reduces to 


m RT E m 

(V-)" ■ ( R c T c> 


(58) 


E E 

To simulate turbine cooling conditions (T m , T c ), a hot gas static 

M 

temperature at injection, T w , can be determined from Eq. (58) by 

M 

specifying a reduced coolant temperature, . Knowing the Mach number 
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at Injection (assumed to match the prototype) a total gas temperature* 

T t M > can be determined. Using Eq. (29), a reduced total pressure, P* M , 

03 MM w 

is uniquely specified. With P T , Ty and the local Mach number dis~ 

CO l co 

tribution specified, similarity in heat transfer between the model and 
engine conditions is established. 

E M 

By requiring an equivalence of I and 1 , sizing of the injection 
hole diameter relative to actual hardware can be accomplished. To main- 
tain proper dynamic interaction between the coolant jets and the boundary 
layer, the ratio of the drag forces exerted on the issuing jets to the 
change in the coolant momentum flux from ‘ s plenum to the exit of the 
coolant hole must be maintained. Assuming that the coolant jet trajectory 
is influenced by drag forces indicative of the local momentum- thickness 
boundary layer, the following can be written 


.(pjd^d,,) E 

L "5 5 J 

_w_., 

(pcW 

where 0.. = momentum boundary layer thickness at injection 

C M 

However, since I b = I , Eq. (59) can be expressed 

(0,/d/ ■ (V d o )M 

At similar locations along the film cooled surface. 


(59) 


(60) 


( 61 ) 


Substituting Eq. (61) into (60), determines the model coolant hole 
diameter as 



” C MF d O 


( 62 ) 


Under the condition stated in Eq. (62) and assuming that the 
viscosity of the air varies as 


(rpM = (r 1 -)” < 63 > 

*Vef ref 

M E 

over the range T c < T < 5 the coolant Reynolds numbers based on the 

coolant ho\e diameters are equivalent. 
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For a simulation in which L ' /d Q remains close to the actual wall thick- 
ness-to-hole diameter ratio for a turbine vane, Eq. (64) establishes the 
necessary condition for correctly modeling the dimensionless pressure 
drop and heat transfer through the coolant holes. 

A mathematical justification has been presented for those film 
cooling similarity parameters previously noted to be important In model' 
ing gas film cooling configurations (Chapter II. Literature Survey). 

T he attention of the reader is now shifted to a description of the sub- 
ject experimental investigation in which the concepts of Chapter III. 
(modeling) are applied to assess the cooling effectiveness of a proposed 
turbine-vane, leading-edge, cooling configuration. 
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IV. THE EXPERIMENTAL INVESTIGATION 
IV. 1. Experimental Apparatus 

One of the major objectives of the present investigation was the 
simulation of the high pressure and high temperature turbine environment 
(as specified by Chapter III) suitable for modeling the heat transfer to 
a film cooled turbine vane, A geometric modeling factor was chosen such 
that the test vane was a 3X size representation of a prototype vane 
design for a high temperature and high pres p ure gas turbine engine. 

Figure 17 shows the prototype vane contour described by the tabulated 
nondimensional coordinates. The model leading edge radius was 0.0152m 
(0.6 in) with a vane span of 0.114m (4.5 in). A two-dimensional test 
section was constructed housing the suction and pressure surfaces of the 
vane model. The test section was located downstream from a hot-gas 
generator, capable of matching the test section inlet Mach number, 
Reynolds number and turbul ence 'intensify with tnose characteristic of 
the prototype turbine inlet conditions. The details of the vane model 
were such that film coolant could be injected through a single row of 
holes with three different locations of the row relative to the stagna- 
tion point (3 locations on the suction surface and 3 on the pressure 
surface) . 

The following sections present a complete description of the 
capabilities of the experimental system, a description of the overall 
flow system, the turbulence generating inlet section, the two-dimensional 


R = 0.0015m (0 . 6in) 
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Figure 17. Dimensionless Coordinates of Vane Model 
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vane channels and a review of the various measurement techniques em- 
ployed. 


IV. 1.1. General Flow System 

Figure 18 is a simplified schematic of the overall flow system 
used In conducting this investigation. A high pressure air supply 
delivered air at pressures from 9.7 - 11.4 x 10 5 N/m 2 (139.7 - 164.7 
psia) with mass flow rates up to 7.73kg/sec (17 Ibm/sec). To provide 
control of the temperature to the test channel, the supply air was 
split into two j. ths: a by-pass line and a supply line for a natural 
gas combustor. Combustor exit air temperatures up to 1033°K (1860°R) 
were available with mixing of the combustor exit and by-pass air to 
provide a wide range of operating temperatures . The mixed heated air 
flowed through an S-shaped duct anc. a straight section 6.10m (20 ft) 
long * 0.254m (10 in) diameter, and then through a turbulence generation 
assembly before entering the vane channel test section. The flow into 
the vane test section was split into three legs: two .152m (6 in) by- 
pass Tines and a ,208m (8 in) mainftow line that we re indi vidually 

throttled before exhausting to a collection manifold. The manifold air 
was throttled prior to exhausting intc an altitude chamber capable of 
back pressures from .254m (10 in) of mercury to atmospheric conditions. 
Figure 19 illustrates the flow system upstream from the test section. 
General features (Indicated by the numbered arrows) include the natural 
gas combustor (1), the dilution-air supply line (2) and the S-duct lead- 
ing to the test section (3). Figure 20 shows the turbulence generator 
section (1) and the vane channel (2) with the side cover-plate removed. 






Photograph Showing Upstream Combustor (1) 
Dilution Air Supply (2) and S-duct (3) . 
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IV. 1.2. Turbulence Generation Assembly 

The turbulence generation assembly was designed to permit variation 
of the vane channel inlet freestream turbulence intensity from low to 
relatively high levels by utilizing screens of varying sizes located at 
selected positions upstream from the vane channel. Turbulence intensity 
measurements, using a hot film probe and anemometer, were made at the 
test section inlet as a function of screen size and eight upstream 
locations. Based on the screen wire diameter, d^, the nondimensional 
distance from the screen to the hot film probe, x/d^, could be varied 
from x/dj,j - 30 to x/d^ = 335. Using the theory for decay of isotropic 
turbulence downstream from screens (45, 46), the turbulence intensity 
level (Tu^) created by the turbulence generation assembly was estimated 
to range from 1.4% to 12%. 

Two types of screens were used in this investigation: (1) a screen 
with 1.59 x 10~ 3 m (0.0625 in) diameter wires (cy at 7.94 x 10" 3 m 
(0.3125 in) center- to-center spacings (M s )> and (2) a 3.18 x 10" J m 
(0.125 in) diameter wire mesh with 6.36 x 10 m (0.250 in) center-to- 
center spacing. The full details of the design of the turbulence section 
can be found in Appendix B.l. 


IV. 1 .3. Test Section 

The 3X size vane test section consisted of the two-dimensional 
channel formed by the suction and pressure surfaces of two adjacent tur- 
bine vanes. Figure 21 is a close-up illustration of the vane test 
section shown as item (2) .in Figure 20. Figure 21 shows the test section 
(with cover plate removed) comprised of three flow channels: two by-pass 

legs (1, 2) and a main flow channel (3) formed by the suction surface of 
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one vane and the pressure surface of the adjacent vane. This channel 
design was chosen for its ease in instrumenting the vane test surface. 

As shown, the back sides of both vanes were open to the surroundings to 
facilitate access to all vane surface instrumentation. The too by-pass 
channels were throttled independently in order to split the test section 
inlet flow in the right proportions for proper control of the stagnation 
point at the vane leading edge. 

Three main items constitute the test section when assembled: the 

main housing, the individual vane surfaces, and the side cover-plate. 

With the cover-plate removed (see Fig. 21), the vane surfaces were bolted 
to the opposite wall of the housing and to the exit ducts of the three 
flow passages. The too surfaces were positioned with their chord lines 
at 44.5° with respect to the mainstream flow direction, matching the 
orientation of the prototype vane. The cover plate was installed by 
bolting it to both the housing and the vane surfaces. Removal of the 
cover plate allowed immediate access to the vanes for both inspection 
and repair. 

The stainless steel, main housing had a rectangular inlet area of 
0.144m x 0.244m and, at the vane leading edge, was divided in the ratio 
of 1:1:2 for the too by-pass ducts and the vane channel, respectively. 
Three water-cooled access ports (0.0127m (0.5 in) diameter) for intro- 
ducing traversing probes were located at the bottom and side walls of 
the housing, 0.0702m (3 in) upstream from the vane leading edge. 

Figure 22 illustrates the construction details for the vane suction 
surface and Fig. 23 presents similar information for the vane pressure 
surface. The test vanes were constructed of contoured stainless steel 



Figure 22. Construction Details For Suction Surface. 



1.588x10 m (0 . 0625in) 
OFHC Copper Skin 



Figure 23. 


Construction Details For Pressure Surface. 
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inner structures over which copper skins {1.588 * 10 m (0.0625 in) 
thick, oxygen-free, high-conductivity copper) were attached. The copper 

_3 

skins were furnace brazed to the inner bodies along 1.588 x 10 m 
{0.0625 in) wide lands that separated water coolant channels running 

_3 

along the span of the vanes. The coolant channels were 1 .588 x io m 
(0.0625 in) deep and 0. 0127m (0.50 in) wide (along the surface contour) 
and carried a water flow of 0.063 to 0.126 liter/sec (1 to 2 gal/min) 
per channel. The coolant channels were closely spaced along the surface 
contours of both vanes and were individually throttled to provide 
variation of heat removal from the copper test surfaces. This design 
insured the capability of maintaining a uniform wall -temperature for 
the heat transfer experiments and, at the same time, permitted the vane 
surfaces to operate at a temperature level that simulated the freestream 
gas**to-wall temperature ratio of the prototype vane. Moderate wall 
temperatures were also necessary for proper operation of the heat flux 
sensors selected for use in this study. 

Figure 24 illustrates another basic feature of the test vane design. 
Both the suction and pressure surfaces had removable leading edge 
sections. Three leading edge pieces were designed for each of the 
suction and pressure surfaces such that film coolant could be injected 
at three different locations relative to the stagnation point. Figure 
24 shows typical cross-sectional views of the leading edge pieces, 
illustrating the film coolant supply line, the film coolant plenum 
chamber, the water coolant channels, and the stainless steel inner 
structure. The leading edge pieces for the pressure surface had film 
coolant plenum chambers located at 0.0054m (0.214 in), 0.020m (0.783 in) 
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and 0.0315m (1.24 in) from the stagnation point. The configuration for 
film coolant injection at = 0,0054m is shown in Fig. 24, The leading 
edge pieces for the suction surface had plenum chambers located at 
0.017m (0.669 in), 0.0307m (1.209 in) and 0.047m (1.831 in) from the 
stagnation point. 

IV. 1.4. Instrumentation 
IV. 1.4.1. Mass Flow Measurements 

Flowrates for the by-pass air, fuel, and combustion air and the hot 
gas flowrate through each leg of the vane channel were measured using 
thin plate orifices constructed to ASME standards. Film coolant flow- 
rates were measured with small diameter venturi meters for low and 
intermediate coolant flowrates and turbine flowmeters for higher coolant 
flowrates. 


IV. 1.4. 2. Inlet Flow Conditions 

Inlet flow conditions to the vane channel were measured in a single 
plane perpendicular to the flow 0.0762m (3 in) upstream from the vane 
leading edges. Vertical and horizontal traverses in that plane were 
made with wedge-shaped, pitot-static and total temperature probes to 


map the total pressure, static pressure, and total temperature. The 


vane test section inlet total pressure (Py ) and total temperature (Ty ) 

co co 


were defined as the mass average of the measured values averaged over 


the sector of the inlet plane which represented the flow viewed by the 
central vane test channel. This "averaging sector" extended from 


0.3 aiy/H ^ 0.7 and 0.2 z/W ± 0.8 in the measurement plane upstream 
from the test vanes (where H is the duct height and W is the spanwidth). 
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Measurements of the vane channel inlet turbulence intensity were 
made prior to the initiation of the vane film cooling test program. A 
hot film probe (Thermo -Systems, Incorporated (TSI) 1220-20) was used to 
measure inlet turbulence intensity levels, generated by upstream grids, 
over a range of inlet Reynolds number. Additional details of the 
measurement inlet flow conditions are presented in Appendix B. 

IV. 1 .4.3. Test Surface Temperature Measurements 

Wall temperature measurements were made using small diameter, 
swedge-type copper constantan thermocouples. The thermocouples were 
furnace brazed into designated hole locations on the vane test surface 
and were finished flush to the wall. Exposed-ball or bead-type thermo- 
couples were used to measure coolant temperature in the film coolant 
plenum chamber. The exposed couples were positioned with the junctions 
at the center of the plenum chamber. Further details of the wall tem- 
perature instrumentation are included in Appendix B. 

IV. 1.4. 4. Test Surface Wall Pressure Measurements 

Vane surface static pressure measurements v/ere made with flush 
mounted copper tubes brazed into the vane suction and pressure surfaces. 
Ports to monitor film coolant plenum pressure consisted of holes drilled 
through to the plenum and sealed on the back side with larger diameter 
stainless steel extension tubing. Details of the static pressure 
instrumentation are included in Appendix B. 
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IV. 1.4. 5, Test Surface Heat Flux Measurements 

Direct measurements of the vane surface heat flux were made using 
miniature, Gardon-type, thin-foil heat flux gages (Model No. 2000, 
Thermogage, Inc.). The basic operating principle for the thin-foil 
Gardon transducer is illustrated in Fig. 25. A thin, circular constantan 
foil is fused to the end of a tubular, copper heat sink, forming a con- 
tinuous cop per- co ns tan tan thermocouple junction along the circumference 
of the foil. On the underside of the foil, a fine copper wire is joined 
to the center of the foil, forming a second thermocouple junction. Under 
operating conditions, the incident heat flux travels radially along the 
foil to the copper heat sink, establishing a parabolic radial temperature 
profile with a temperature difference from the center to the edge of the 
foil that is dependent on the magnitude of the incident heat flux (47). 
Due to the physical property variations of the constantan foil, and the 
change in the junction emf between the two copper- foil junctions, the 
differential emf between the two foil junctions becomes linearly depen- 
dent upon the magnitude of the incoming heat flux (48). 

In order to minimize vane surface irregularities with gage instal- 
lation and, at the same time, approach point-source measurements, a 
miniature gage design was chosen with a sensing surface dimension small 
in comparison to the characteristic dimensions of the vane test surface. 
Figure 26 is a detailed drawing of a heat flux sensor and the installa- 
tion technique used in this investigation. The basic sensor was pressed 
into a 0.0037m (0.145 in) threaded (5-40 UNC) oxygen-free, high-conduc- 
tivity copper plug. The sensor area itself was 5.08 x 10~\ (0.02 in) 
in diameter with a foil 2.54 * 10"^m (0.001 in) thick. A 5.08 x 10” 5 m 
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Figure 25. Details of Thin-Foxi, Gardon 
Heat Flux Transducer. 
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(0.002 in) diameter copper wire was laser-welded to the center of the 
foil from beneath. This wire and a 1.27 * 10 m (0.005 in) diameter 
copper, ground-wire (attached to the heat sink) were both coated with 

_ r 

5.08 x 10 m (0.002 in) of teflon, fully shielded, and brought out the 
bottom of the threaded plug. The heat flux gage was installed into the 
vane test surface by threading into an oxygen-free, high-conductivity 
copper plug mounted in the test surface wall. Each mounted gage was 
located within a water coolant channel, providing the maximum in heat 
sinking capability. The threaded portion of the gage was long enough to 
provide 10 to 11 threads of contact with the water-cooled plug. 

With the sensor designed for a maximum temperature difference 
between the center and edge of the sensing foil of only 25°F, the maximum 
output level of the gage was in the microvolt range. Each gage, there- 
fore, had an individually calibrated operational amplifier with a gain 
near 1000, to convert the transducer microvolt signal to a more conve- 
nient millivolt signal. The gages were calibrated with the amplifiers, 

-3 

and the gain adjusted to obtain a sensitivity of 8.64 * 10 mv/joul e/sec 
- m^ (10 mv/BTU/sec-ft^) . 

Due to the extremely low level of the sensor output, caution was 
exercised in isolating the entire system from such environmental factors 
as random noise and changes in the ambient temperature. All solder 
connections were made with low-noise, cadnium-tin solder to eliminate 
any additional emf sources between the sensors and the amplifiers. 

The amplifiers have a characteristic drift of 3 millivolts per 10°F 
change in the ambient temperature (a zero shift and not a calibration 
change). Therefore, all amplifiers and associated electrical connections 


were thermally isolated from ambient conditions using a sealed, styrofoam 
enclosure. 

Calibration of the sensors followed the technique depicted in Fig. 

27 and was performed by the sensor manufacturer. Thermogage, Inc. Each 
sensor was mounted on the center axis at the entrance to one of two, 
back-to-back (symmetrical) black-body cavities. The black body was a 
graphite assembly, tapered to achieve a uniform inside wall temperature 
with resistance heating. The heat flux source-temperature was monitored 
on the opposi J , „ side of the apparatus in a similar black body cavity 
using an optical pyrometer, with a calibration traced to the National 
Bureau of Standards. The length-to-diameter ratio of the cavity was 
large enough, coupled with the high emissivity of the walls, to simulate 
an emissivity of near unity in the vicinity of the sensor. A view factor 
between the black body opening and the sensor surface can be easily 
calculated and, knowing the radiative properties of the previously 
coated sensor, an accurate calibration was performed (49). Output from 
a typical calibration test is shown in Fig. 2B. The calibration was 
accurate to within ± 3 % of full scale with repeatability within 1 / 2 % 
and a typical sensor response time of 2.4 milliseconds. 

Interpretation of the sensor output requires some qualification. 

In order for the calibration to be valid, the sensor body temperature 
must be greater than 283°K (510°R) and loss than 506°K (910°R). Outside 
this range, the change in emf of the thermocouple junctions is not com- 
pensated by the physical property changes of the foil to produce the 
typical linear calibration curve. Secondly, the gages must operate in 
a field of heat flux uniform over an area that is much greater than the 
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sensor area. These two restrictions posed no difficulty in conducting 
this investigation. 

IV. 1.4.6. Instrumentation Locations 

Except for the heat flux sensors, all vane surface instrumentation 
was located at the lands separting the water coolant channels. At 
various stations along the vane surface contours, pressure taps and 
thermocouples were installed to determine the two-dimensionality of the 
flow and the uniformity of the wall temperature boundary condition. The 
heat flux sensors were positioned between lands at the center of the 
coolant channels. The requirement of maximum heat sinking for the 
sensors determined the spacing between gages in the chordwise direction. 
No more than one gage was installed per coolant channel to minimize 
pressure loss in the water coolant channel. Sensor mounting plugs were 
approximately one half the width of the coolant channel. The sensors 
were staggered with respect to the center span line to minimize the 
possible disturbance of local heat flux due to the influence of an up- 
stream gage. 

Tables 1, 2, and 3 list the locations of all suction and pressure 
surface instrumentation giving the coordinate x/L s relative to the 
stagnation line for the three film coolant plenum configurations em- 
ployed. Lateral (spanwise) spacing of the pressure taps, thermocouples 
and heat flux gages, z/L s , are referenced with respect to the center 
span line. Film coolant pleniffn instrumentation is also referenced to 
the center of the plenum. 
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0.78 

0.14 





14 

0.84 

0.11 





0.78 

0.00 





15 

0.84 

0.00 

*z is referenced to cents rspan 

0.78 

-0.14 










PIENC-1 IKSTRiT43m^nCN 






1 


0.24 


0.26 




0.24 


0.26 



2 


-0.03 


0.05 




-0.03 


0.05 



3 


-0.24 


-0.26 




-0.24 


-0.26 




Table 2 


Pressure, ^temperature and Heat Flic: Rfeasureirent locations For Suction 
and Pressure Surfaces For Oonfigurations of x-jyts = 0.150 and x^/Lp = 0.121 


SUCTION SURFACE ( y^/L ^ = 0,228, = 0.020m) PRESSURE SURFACE C x^/Lp = 0,193, Lp = 0.016m) 



Pressure 

Terrpsrature 

Heat Flux 

Pressure 

Tenperature 

Heat Flux 


X/Lg 

z/%* 


z/Lg 


Z/Lq 

y/Lp 

z/Lp 


z/Lp 


2 Aj? 

1 

0.03 

0.16 

0.03 

-0.16 

0.32 

0.00 

-0.05 

0.20 

-0.05 

0.09 

0.32 

0.00 

2 

0.03 

-0.06 

0.03 

0.07 

0.39 

0.02 

-0.05 

-0.08 

-0.05 

-0.21 

0.40 

0.03 

3 

0.11 

0.05 

0.11 

-0.05 

0.46 

-0.02 

0.04 

0.26 

0.04 

-0.07 

0.48 

-0.03 

4 

0.11 

-0.02 

0.18 

-0.04 

0.53 

0.05 

0.04 

0.06 

0.14 

-0.05 

0.57 

0.06 

5 

0.18 

0.15 

0.18 

0.09 

0.67 

-0.05 

0.14 

0.19 

0.14 

-0.11 

0.74 

-0.06 

6 

0.18 

0.03 

0.36 

0.13 

0.81 

0.00 

0.14 

-0.04 

0.36 

0.16 



7 

0.18 

-0.15 

0.36 

0.00 



0.14 

-0.19 

0.36 

0.00 



8 

0.28 

0.02 

0,36 

-0.13 



0.27 

0.02 

0.36 

-0. 15 



9 

0.43 

0,11 

0.50 

0.00 



0.44 

0.14 

0.53 

0.00 



10 

0.43 

0.00 

0.64 

0.00 



0.44 

0.00 

0.69 

0.00 



11 

0.43 

-0.11 

0.77 

0.00 



0.44 

-0.14 

0.87 

0.00 



32 

0.57 

0.00 

0.91 

0.00 



0.61 

0.00 





13 

0.72 

0.00 





0.78 

0.14 





14 

0.84 

0.11 





0.78 

0.00 





15 

0.84 

0.00 





0.78 

-0.14 





16 

0.84 

-0.11 











17 

-0.03 

0.0 *2 

is referenced 

to canterspan 
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1 


0.24 


0.26 




0.24 


0.26 



2 


-0.03 


0.05 




-0.03 


0.05 



3 


-0.24 


-0.26 




-0.24 


-0.26 




Table 3. Pressure, Tanperature and Heat Flux Measurement Locations For Sucticr, 

and Pressure Surfaces For Configurations of xj_/Lg = 0.228 and xj/Lp = 0.193. 
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IV. 1.5. Measurement Systems 

All pressure measurements were made using two methods: strain-gage 

type transducers or a digital automatic multiple pressure recording 
system. Each thermocouple was referenced to a 338°K (610°R) oven. 

The output from each heat flux sensor was amplified using operational 
amplifiers with a gain near 1000. 

All data was digitally sampled and linked to an on-line computer 
as well as a central IBM 360 processing center. The on-line computer 
provided continuous information to a visual display unit for purposes 
of monitoring current flow conditions. 

A detailed description of the measurement system including raw 
data conversion techniques, data monitoring, acquisition, storage and 
processing is presented in Appendix B. 

IV. 2. Description of Experiment 

The reduction in the local heat flux due to film cooling from a 
single row of laterally angled holes in the leading edge region of the 
model turbine vane was measured under reduced flow conditions simulating 
a high temperature, high pressure turbine environment. The test surface 
was a 3X scale model of a prototype turbine inlet vane designed for 
operation at 1922°K (3000°F) inlet total temperature and 3.04 * loSl/m 2 
(441 psia) inlet total pressure. Simulation of these high temperature 
and pressure conditions was achieved in the subject investigation at a 
model inlet total pressure and total temperature of 3,1 x 10 5 N/m 2 
(44.67 psia) and 698°K (1240°R), respectively, with a vane surface Mach 
number distribution that closely duplicated the prototype design. An 
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inlet Mach number of 0.22 and Reynolds number of 8 * 10 4 (based on the 
vane leading edge radius (0.0152m (0.6 in)) were employed throughout the 
film cooling experiments. The model inlet turbulence intensity (TuJ was 
held constant at 6 . 6 % by placement of a fine mesh screen (Mg/d^ - 5.0) at 
a location 0.533m (21 in) upstream from the model vane leading edge. The 
vane wall temperature was maintained at approximately isothermal conditions 
with T t /T w = 2.15, thus approximating a boundary layer with large density 
variations. Using ambient cooling air the dimensionless coolant temper- 
ature, BjL, was near 1,03. With T c at ambient conditions (294°l< or 530°R), 
the inlet temperature (T^ ) of 693°K (1240°R) resulting from Eq. (58) 
satisfies the modeling requirement that M and I be maintained equivalent 
between the model and the prototype vane. 

The film coolant injection geometry was chosen for possible applica- 
tion to future turbine designs. Coolant hole diameters were specified at 

- Q 

1.17 x 10 m (0.046 in) in an attempt to keep the local momentum boundary 
layer thickness -to -hole diameter ratio ( G^/cIq) similar to expected engine 
conditions. Typically, at high temperatures and pressures, may vary 
from 5 x io** 6 m (1,97 x IQ"' 4 in) near x/L s - 0 (stagnation) to 2 x 10“^m 
(7.87 x TO -4 in) at x/L s = 0.20 (50). For a typical coolant hole diameter 
of 3.8 x 10~ 4 m (0.015 in), 0|/d G may range from 0.013 to 0.0525 in the 
leading edge region. Length -to-diameter ratios (L‘/d 0 ) of the coolant 
hole passages were specified at 3.0 for 3 = 35° and 4.4 for 3 - 18°. A 
surface spacing-to-hole diameter ratio S^/dg of 4.0 was selected to avoid 
highly non-uniform, spanwise film coverage. 

The film cooling effectiveness of the subject cooling configuration 
was determined by comparing the measured heat flux with film cooling. 
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qpQ* to that without cooling, qjj, at dimensionless distances from injec- 
tion, (x-x.j)/d 0 , ranging from 10.45 to 126,11. The heat flux reduction 
was determined as a function of the blowing ratio, M, ranging from M = 

0.1 to M = 2.2. To determine M (M = p^Vq/p^V^ the coolant mass flux was 
assumed uniform across the span of the holes and equal to the total 
coolant mass flow divided by the total area of the coolant holes (19 
coolant holes across a span of 0.114m (4.5 in)). The hot gas mass flux 
was determined from density and velocity values at injection, calculated 
from total gas temperature, total pressure and wall static pressure 
measurements. 

Two injection angles were investigated: 0 = 18° and 0 = 35°. 

Injection was from a single row of holes on both the suction and pressure 

surfaces with three different locations for the holes being investigated. 

On the suction surface, injection occurred at x/L s values of 0.0833, 

0.150 and 0.228 corresponding to 0^/dg values of 0.013, 0.036 and 0.0514 

with acceleration parameter, K^, values of 7.9 x 10" 6 , 0.723 x io“ 6 and 

0.656 x 10" , respectively. On the pressure surface, injection occurred 

at x/Lp values of 0.0333, 0.122 and 0.199 corresponding to 0./dg values 

of 0.0111, 0.0514 and 0.0827, and K,. values of 6.6 x lo" 5 , -6.0 x 10' 5 
-6 

and 4.267 x lo , respectively. These coolant hole locations were 
selected to represent possible regions of considerable curvature in 
which streamwise injection schemes may prove difficult. 

Table 4 lists in detail a direct comparison between the frees tream 
flow and film coolant parameters for a high temperature, high pressure 
engine environment and those conditions maintained in the subject model 
study. 
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PARAMETER 

ENGINE 

MODEL 

T , °K 

1922 

698 

P , N/m 2 

CO 

3.04xl0 6 

3.1xl0 5 

V ° K 

800 

294 

Ma x 

Design 

Design 

T T / T W 

! CO 

2.0 

2.15 

Tu„, * 

10-20% 

6.6% 

Ma 

ilct INLET 

0.22 

0,22-0.23 

L.E. Radius, m 

0.0051 

0.0152 

| Re LER, INLET 

8xl0 4 

8xlQ 4 


INJECTION PARAMETERS 


(PV) C /(PV) 

0 2.5 

0 + 2.5 

Injection Angle 

15° -> 90° (a or 0) 

18°, 35° (0) 

s a/ a o 

3.0 ■+ 8.0 

4.0 

d Q , m 

2 * 3.8xl0 -4 

1.17X10" 3 

Lyd 0 

3,0 h- 5.0 

3.0, 4.4 

V a c 

0.013 0.0525 

0,011 ->• 0.0827 

X i /L S,P 

0 + 0.25 

0.083 *> 0,228 


Table 4. Comparison of Main Plow and Film Cooling 
Parameters For a Gas Turbine Engine With 
Those of the Subject Investigation* 


OBlGMMi £ 


Vvjjj ' 














































The actual mechanics of recording the film cooling data are dis- 
cussed in Appendix B.4, The essential features of the procedure include: 
(1) a horizontal and vertical traverse of the inlet plane to define the 
freestream total and static pressures and total temperature for each heat 
flux survey, and (2) a survey of the measured wall heat flux with the 
probes retracted to eliminate the effects of the probe wakes. 

IV. 3, Data Analysis 

Heat flux measurements were made downstream of injection at locations 
ranging from (x-x.j)/dQ of 10.45 to 126,11 for the suction surface and 
(x-x^J/dg of 10.35 to 93.06 for the pressure surface. Heat flux levels 
with film cooling at these locations were compared to measurements at 
the same locations with no coolant flow. From this comparison a frac- 
tional reduction in the heat flux due to coolant flow was calculated and 
presented as a measure of the cooling effectiveness. 

Under normal boundary layer conditions, the local heat flux to a 
surface can be written 

"o " h o - V W 

where qjj . - local heat flux without film cooling 

hg = local convective heat transfer coefficient 
T t = total gas temperature 

CO 

T^ = local wall temperature 

With the addition of film cooling, the local heat flux can be expressed 
in a similar manner. 


* T 
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q FC = h FC < T T " T W,FC* 

CO 

where 

qji c = heat flux with film cooling 

hp C - local heat transfer coefficient with film cooling 

T w p C “ local wall temperature with film cooling 

Equation (66) combines all effects of the film injection phenomena into 
the single parameter, hp C , which can be expected to depend on the follow- 
ing variables 

hpc " hp C (M, V c /V„, |3, 6^, Tu to , K^, etc.) (67) 

This is in contrast to the use of an adiabatic wall temperature to char- 
acterize film cooling (5,44) in which the local film temperature is 
presented in terms of an adiabatic film cooling effectiveness while the 
heat transfer coefficient reflects only the hydrodynamic effects of 
coolant injection. 

The film cooling results are presented as a fractional reduction in 
the local Stanton number where the local Stanton number is defined as 



Since the frees tr earn parameters do not change with film cooling (i.e. 
p M , V w , C remain constant with or without coolant injection), Eq. (68) 

rcg 

can be used to express the local Stanton number reduction due to film 
cooling in the following manner 


C 69) 


st fc <>fc (T t„ - W 

ST o - qj 


Under isothermal conditions, g = T^ pg (which holds for the 
highly cooled walls in this investigation), the heat flux data give the 
Stanton number ratio directly. Thus, 


ST FC _ q FC 
ST 0 


(70) 


The isothermal film cooling effectiveness, a is defined as 

n - q o - q rc 
n iso “ ~ ? 0 


which gives the Stanton number reduction as 


(71) 


ST 

ST 



n IS0 


(72) 


Having described the objectives, methodology and details of the 
subject investigation. Chapter V. will now present a summary of the 
experimental results. These results are divided into three categories: 
(1) documentation of the turbine environment, (2) non-film cooled (dry- 
wall) vane heat transfer, and (3) actual film cooling heat flux reduc- 
tions, A detailed discussion is also included for each of the above 
items. 



106 


V* RESULTS AND DISCUSSION 

V.l. Documentation of Modeling of Turbine Environment 

The experimental model was employed to simulate both the inlet free- 
stream conditions and the local vane surface conditions representative of 
a high temperature, high pressure gas turbine environment. At the inlet 
to the model vane channel the following conditions were maintained: (a) 
freestream Reynolds number (based on vane leading edge radius) of approx- 
imately 8 x 10 4 , (b) freestream Mach number of 0.22, and (c) an inlet free- 
stream turbulence intensity of 6.6^. Along the vane contour, the vane 
local conditions were as follows: (a) Mach number distribution approx- 

imating the prototype vane design calculations, (b) nearly uniform wall 
temperature with T^. /T^ near 2.15 (representative of density variations 

CO 

through the boundary layer for a gas turbine vane), and (c) approximately 
constant dimensionless coolant temperature, 6£, of 1.03 (representative 
of cool ant- to-gas density ratio typical of high temperature turbine 
operation). 

Simulation of the inlet conditions for the vane prototype was achieved 
at a reduced inlet total pressure and total temperature of 3.1 x io® N/m^ 
(44.67 psia) and 689°K (1240°R), respectively. Under these conditions the 
following inlet velocity and temperature profiles were measured. 

Figure 29 shows the nondimensional inlet velocity profile in the 
vertical direction. The inlet velocity, as compared to the sector-average 
velocity, is shown as a function of the duct height. The sector average 
velocity, representative of the flow through the vane channel, was defined 


Dimensionless Inlet Duct Height, y/H 



0.7 0.8 0.9 1.0 1.1 1.2 1.3 

Dimensionless Inlet Velocity, v / v se02QR ^VE 


Figure 29. Variation of Inlet Velocity 
{Vertical Traverse) . 
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as the mass average velocity for the sector, 0.3 ^.y/H £.0.7 and 
0,2 <. z/VJ £ 0.8. Within the sector limits, the velocity did not vary 
more than B% from the average. All vertical traverses were made at a 
fixed position spanwise in the channel at z/W =0.5. 

Figure 30 shows the variation of the nondlmensional inlet velocity 
in the horizontal direction across the span of the duct. The profile is 
less flat in the spanwise direction than shown in Fig. 29 possibly due 
to a contraction in the lateral direction (2,22:1 area ratio) as the 
flow transitioned from a 0.254m (10 in) 0D circular pipe to the 0.244m 
0,114m (9.6 in * 4.5 in) test channel. Within the ’'sector-average" 
limits defined above, the velocity is shown to vary by as much as 10" 
from the sector average. This variation did not have a significant 
effect on the local vane surface measurements such as wall heat flux or 
calculated local velocity. The total and static pressure at the inlet 
plane was uniform within 1/2% of the calculated mean values. 

Figure 31 presents the variation of the measured inlet total tem- 
perature, nondimensionalized with respect to the sector mass average 
total temperature, as a function of the inlet duct height. For the data 
shown, the inlet total temperature was uniform to within 1/2% of the 
sector average. Figure 32, presents similar data for the horizontal 
temperature variation, which shows some nonuniformity due to the lateral 
contraction from the initial circular duct. However, the deviation of 
tiie temperature in the spanwise direction (for the sector limits defined 
above) was never more than 1/2% of the sector average. 

Based on the sector average values of the inlet velocity and total 
temperature shown previously, the experimental values of the inlet 
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Reynolds number and the inlet Mach number ranged from 8.0 to 8.4 * 10 , 
and 0.28 to 0.23, respectively. The prototype vane row is designed to 
operate at an inlet Mach number of 0,22 and an inlet Reynolds number of 
8 x 10^ (based on leading edge radius) with an inlet total temperature 
of 1922°K (3460°R), inlet total pressure of 3.06 x 10 6 N/m 2 (441 psia), 
and 65.67 kg/sec (122.5 Ibm/sec) total air flow through the 36 vane ring. 

To complete the documentation of the vane channel inlet conditions, 
the inlet turbulence intensity, Tu^, was measured 0,076m (3 inches) up- 
stream of the vanes as a function of various screen types and positions. 

A hot film probe was oriented perpendicular to the flow at the midspan 
of the channel and a vertical traverse of the inlet turbulence intensity 
was made. A sector average intensity was calculated from those measure- 
ments for the sector between 0.3 < y/H <. 0.7. 

The sector average turbulence intensity at the vane channel inlet 
are presented in Fig. 33 as a function of distance from the upstream 
screen. Intensity values of 6.6% to 10.55J were measured at inlet Reynolds 
numbers of 4.34, 6,51, 9.98, 11.29 x 10^ (corresponding to inlet veloc- 
ities of 30.5, 45.7, 70.1 and 79.2 m/sec). As suggested by Von Karman 
(51), the intensity date- was found to correlate with the nondimensional 
distance from the grid based on the grid wire diameter rather than the 
grid mesh length. Three significant conclusions were drawn from the 
turbulence study; (1) the absolute level of the measured intensity is 
higher than values reported by Baines and Peterson (52) and Eckert, 

Sparrow and Newman (53) for turbulence decay downstream of a grid as 
shown in Fig. 33, (2) the variation of the gradient of the one-dimensional 
turbulent kinetic energy is representative of typical decay downstream 
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Figure 33. Grid-generated Inlet Turbulence Intensities. 
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from generating screens, and (3) the turbulence level is independent of 
the inlet Reynolds number over the range investigated, 

A least-squares curve fit of the data (excluding the 70.1 m/sec 
results) in Fig. 33 is represented by the solid curve. The agreement 
with correlations reported for similar investigations (52,53) is poor, 
especially in the far downstream regions from the generating grid. 

The only significant difference betv/een the present data and the 
cited references is the character of the turbulent field prior to inter- 
action with the generating grid. Clear tunnel intensities (no grid 
present) were generally reported from M!% to \% (52,53). In the present 
study, clear tunnel intensities near 7% were measured. Since no provi- 
sions were made in this test facility to minimize clear tunnel conditions, 
it is reasonable to expect high overall values of turbulence. The close 
relationship between initial conditions and decay downstream from a grid 
was clearly demonstrated by Tsuji (54). Turbulence produced by a grid 
was shown to be superimposed on the high wave number range of that tur- 
bulence incident upon the grid. The decay of the spectral components 
after the initial stages most nearly resembles that of the initial 
intensity conditions. Such a progression of events seems to be sub- 
stantiated in Fig. 33 as the intensity level for x/d w > 100 appears to 
have approached that of the clear tunnel. 

The abnormally high values of Tu at an inlet velocity of 70.1 m/sec 
(230 fps) were explained by a hot film probe resonance due to a vortex 
shedding phenomenon. At 70.1 m/sec, the Reynolds number based on the 
hot film probe wire diameter is 225 (a region in which vortex shedding 
is significant). At this Reynolds number, the Strouhal. number 


(Str = w'd/V w , where to 1 is frequency) is near 0.19, yielding a shedding 
frequency of about 40,000 Hz. During these measurements the turbulent 
signal was monitored with an oscilloscope. At 70.1 m/sec, the signal 
developed a large scale periodic pattern with a frequency of 20K Hz. 

At 79.2 m/sec, a similar regular signal appeared at 25K Hz, but of a 
much smaller amplitude. 

The cantilevered probe in an oscillatory flow constitutes a non- 
linear vibratory system in which harmonics will naturally exist. The 
phenomenon of subharmonic response shows that the system may oscillate 
with a frequency of 1/2, 1/3, ..., 1/N of the applied forcing frequency. 

It appears that the forcing frequency due to vortex shedding at 40,000 Hz 
was able to induce and sustain subharmonic resonance of the probe at 1/2 
of its forcing frequency. At vortex shedding frequencies greater than 
approximately 50K Hz, the forcing function is unable to sustain any 
induced vibration of significant amplitude that may affect the RMS signal. 
It was concluded that the data shown for a velocity of 70.1 m/sec is not 
representative of the grid generated turbulence. 

Although the turbulence intensity levels are significantly higher 
than measured in other investigations, the gradient of the one-dimensional 
turbulent kinetic energy with respect to distance from the grid in the 
initial decay regime (just downstream of the grid) is well within the 
range of published values. Figure 34 shows the gradient of the turbulent 
kinetic energy in the main flow direction as a function of the dimension- 
less distance from the screen. The present data is represented by the 
solid line, as calculated from the least squares relation shown in Fig. 33 
In addition to the data from references cited in Fig. 33, the data for 
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Dimensionless Distance From Screen, 

Figure 34 . Gradient of One-Dimensional Turbulent 

Kinetic Energy With Distance From Screen 
{Comnarison to Published Results) . 





Kestin and Wood (55) and Van Der Hegge Zijnen (56) are presented for 
comparison. Figure 34 shows that all the cited references, including the 
present investigations report a change in the turbulent kinetic energy 
gradient that is related to approximately the negative third power of the 
dimensionless distance from the grid. The solid line for the subject 
investigation, was determined from the least squares quadratic represen- 
tation of the data in Fig. 33. The minus three power is not representa- 
tive of isotropic decay in the final period, but may describe the influ- 
ence of any grid initial conditions, 

n 

All turbulence data were obtained at a total pressure of 1.46 x 10 

o 

N/m (21 psia) and a total temperature of 294°K (70°F). Figure 33 repre- 
sents data at these conditions for velocities of 30.5, 45.7, 79.2, and 
70.1 m/sec. Within ± 10% of the least squares curve, there is no signif- 
icant influence of the inlet Reynolds number on the magnitude of the 
generated turbulence intensity. 

With the inlet total pressure and total temperature prescribed to 
properly model the engine conditions, similarity in the vane heat transfer 
through matching of the local Reynolds number along the vane contour, can 
be achieved by reproducing the expected prototype Mach number distribution 
over the vane. Figure 35 presents the local Mach number variation with 
distance from the stagnation point for both the pressure and suction 
surfaces of the test vane. Data from a typical test point show that the 
scale model Mach number variation through the test channel matched the 
prototype design values reasonably well with a maximum deviation of 12% 
for the suction surface and 20% for the pressure surface. The large 
variation from design (20%) for the pressure surface occurred at 


Local Mach Number, Ma. 
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Figure 35. Local Mach Number Distribution Along Vane Contour. 
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x/Lp “ 0.15, just downstream from an adverse pressure gradient region. 

In the discussion on dry wall (non-film cooled) heat transfer results, 
to be presented in a later section, evidence points toward possible flow 
separation in this adverse pressure gradient region. 

The uniformity in the wall temperature boundary condition is shown 
in Fig. 36, The gas- to -wall temperature ratio is shown as a function of 
the nondimensional distance from stagnation for both the suction and 
pressure surfaces. The suction surface wall temperature varied no more 
than 5 % from its mean value. The pressure surface wall temperature 
varied, at most 1 OSS from its mean value. This 10 % deviation occurred at 
the pressure surface trailing edge. At this location, the wall thermo- 
couple was actually imbedded in the vane, stainless steel inner structure 
and indicated a higher temperature than those thermocouples within the 
water cooled copper skin. 

Figure 37 illustrates the variation in the film coolant temperature 
both during a given test run (Run Number) and over a range of measure- 
ment periods. The coolant temperature measured in the film coolant 
plenum was corrected for heat pick-up through the vane wall. (See 
Appendix D for a discussion of the coolant plenum temperature calibra- 
tion.) The variation in the dimensionless coolant temperature, for 
a particular day is shown as the solid vertical bands. The band length 
varies from day to day (Run Number to Run Number) due to the range of 
coolant flows investigated. For a wide range of M, the coolant tempera- 
ture shows the greatest variation due to heat pick-up through the vane 
wall. During the entire film cooling study* the dimensionless coolant 
temperature did not deviate more than 7 % from a mean value of 1 .03. 


Gas- Wall Temperature Ratio 









Research Run Number 


Figure 37. Variation in Dimensionless Coolant Temperature With Run Number 
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For a mean gas-to-wall temperature ratio (T^ /Ty) of 2.15, this corre- 

CO 

sponds to an average coolant-to-gas density ratio (p c /p m ) of approximately 
2,18 characteristic of the film cooling experiments reported herein. 

V.2. Non-Film Cooled (Dr.v-Hall) Vane Heat Transfer 

Having matched the appropriate modeling parameters for both the 
vane channel inlet and vane local conditions, heat transfer data without 
film cooling were taken on both the suction and pressure surfaces of the 
vane to: 

(1) compare the measured dimensionless dry-wall heat transfer 
Stanton number with the results from an existing analytical 
boundary layer prediction method, 

(2) investigate the effect of the freestream turbulence intensity 
(Tu^) on the vane dry-wall heat transfer, and 

(3) demonstrate that the heat flux measurement technique v/as 
reliable and durable under severe operating conditions. 

V.2.1. Analytical Prediction Method 

The boundary layer prediction method used for comparison with the 
subject experimental data was developed by W. M. Kays, et a!., of 
Stanford University (57) . A finite-difference technique is used to 
solve the conservation equations, including the turbulent kinetic 
energy equation (TKE) . The basic numerical scheme in the solution of 
the conservation equations follows the Spalding and Patankar method (58). 
The prediction technique allows for the insertion of a variety of models 
for transport processes (with the freedom to modify these models using 
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empirical correlations) or to obtain results directly from the simulta- 
neous solution of the turbulent kinetic energy equation and the other 
conservation equations. 

The present version of the boundary layer program makes use of the 
concept of an eddy viscosity and thermal conductivity, which is helpful 
computationally. A mixing length hypothesis is used for the boundary 
layer turbulent transport in which the length scale varies linearly out- 
side the viscous sublayer. Within the viscous region, the linear re- 
lation is modified using the Van Driest damping function with appropri- 
ately chosen empirical constants reflecting experimentally substantiated 
effects on the sublayer due to external effects such as acceleration, 
blowing, and suction. 

Typical program inputs that are required include initial profiles 
for the velocity, enthalpy and turbulent kinetic energy. Appropriate 
wall boundary conditions are specified on mass addition and enthalpy. 

The mass addition capability can appropriately handle transpiration 
cooling; however, discrete film cooling (through holes) has not been 
modeled. An external pressure gradient can be specified by inputing an 
arbitrary velocity distribution external to the boundary layer. Laminar 
flow can be forced to transition to turbulent conditions by specifying a 
transition momentum- thickness Reynolds number, Re Q Finally, values 
of the freestream stagnation enthalpy, stagnation pressure and specific 
heat must also be supplied. 

A typical program output can provide boundary layer profile pre- 
dictions at each integration station. Hydrodynamic quantities are also 
computed such as the boundary layer thickness (6), the displacement 


thickness (6^), shape factor (H* ) , momentum thickness Reynolds number 
(Re Q ), skin friction coefficient (Cp) and acceleration parameter (K). 
Additionally, the values of the local Stanton number (ST^), heat transfer 
coefficient (h Q ) and energy thickness Reynolds number (Re & ) can be 
obtained. 

The program can be operated in either of two modes. Mode (1) solves 
all the conservation equations with the exception of the turbulent 
kinetic energy equation. Laminar flow is assumed to develop until a 
specified value of the momentum thickness Reynolds number, Re 0 is 
reached (an input variable for Mode (1) with suggested values from 200 
to 300). At this point there is a smooth transition to turbulent flow. 
Mode (Z) accomplishes the simultaneous solution of the conservation 
equations including the turbulent kinetic energy equation, (TKE). 

V.2.2. Suction Surface Results 

To determine the local Stanton number v' iation along the suction 
surface for comparison to predicted value r , measurements of the local 
heat flux were made at 8 locations along the vane contour. By maintain- 
ing the local Mach number distri* . tier similar to the vane prototype 
over a range of total pressures, he«i> flux measurements were recorded 
over a wide range of the local Reynolds number. Six total pressure 
levels were imposed, maintaining a constant Tj /T^ . The inlet freestream 

C3 

turbulence intensity was varied from 6 * 6 % to 10.5/i over the range of 
total pressures investigated. Table 5 summarizes the inlet and local 
conditions for which dry wall heat transfer data were obtained. 


Condition Held Constant 


T « 689°K {1240°R) 

03 

? T /T = 2. IS 

CO 

Design Ma x Distribution 


Conditions Varied 

P„ = 2.43X10 5 4.17X10 5 N/m 2 

03 (35-60psia) 

Tu = 6.6% (Fig. 38) 

8.0% (Fig. 39) 

10.5% (Fig. 40) 


Table 5. Summary of Conditions For Dry-Wall 
Heat Transfer Measurements. 
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Figure 38 shows the variation of the Stanton number along the suction 
surface with the local Reynolds number for an inlet freestream turbulence 
intensity of 6.62. The data shown are represented by a least squares 
curve fit, designated by the solid line. The broken line is the predic- 
tion of the local Stanton number using the boundary layer program. Mode 
(2) (use of TKE equation) was employed in the prediction assuming negli- 

•k 

gible freestream turbulent kinetic energy. 

The data are shown to be distributed approximately ±102 about the 
least squares curve. A comparison between the predicted and least 
squares representation of the measured Stanton numbers shows nearly an 
identical dependency on the Reynolds number, following approximately a 
-0.3 power of Re^. 

The predicted level of the Stanton number is approximately 102 
higher than the least squares fit of the data. .Two explanations are 
offered for this discrepancy: (1) the boundary layer program does not 
include the effect of surface curvature on the development of the hydro- 
dynamic and thermal boundary layers except for ax i -symmetric problems 
and (2) the vane surface may not have been smooth in the initial start- 
ing length region. 

Curvature (with dR'/dx > 0, R' being the body radius) would affect 
the development of the hydrodynamic boundary layer in the same manner 
as a strong positive pressure gradient. Experimental evidence (38,39) 
has verified that sufficiently large values of K (local acceleration 

- 

Predictions using Mode (2) with input freestream turbulence inten- 
sities from 02 to 102 showed no significant effect on the Stanton number 
which is in agreement with previous experiments (33) (34). 
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Figure 38* 
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Surface Stanton Number Variation With Reynolds 
Freestream Turbulence Intensity. 

Number 






128 


t 


parameter) can depress the local Stanton number as the thermal boundary 

layer penetrates the hydrodynamic layer- Values of of the order 

“6 00 -6 
of 10 existed along the suction surface contour. Values of K near 10~ 

have been shown experimentally to depress the local Stanton number 

approximately 3% to 4%. It appears that surface curvature may account 

for a portion of the difference between the measured and predicted 

Stanton numbers* 

In addition to curvature effects, the measured ST^ could be indica- 
tive of a boundary layer that had transitioned to turbulent flow prior 
to the location which is characteristic of a smooth initial starting 
length. To predict such a behavior, the prediction curve would be 
shifted to the left on Fig. 38, resulting in a lower ST^ for a fixed 
Re^. The combined effects of premature transition as well as surface 
curvature would appear to bring the data in closer agreement with the 
boundary layer prediction. 

Figure 39 indicates the measured Stanton number variation (suction 
surface) with Re^ for an inlet turbulence intensity of 8.0%. The data 
is represented by a least squares curve within a ± 15% band. The 
measured STj, correlates with the -0.15 power of Re^ and appears less 
sensitive to the Re^ than the prediction. Agreement of the data with 
the prediction improves at the higher Reynolds numbers, after being 20% 

5 

lower near Re^ - 3 * 10 . 

Figure 40 presents the STj, variation with Re^ for an inlet Tu w of 
10.5%. The data is ± 20% about a least squares curve fit. The measured 
STj^ appears to depend on the -0.15 power of Re^ as in Fig. 39. 


Local Stanton Number , ST 



Figure 39 . Suction Surface Stanton Number Variation With Reynolds Number 
at 8.0% Frees tream Turbulence Intensity. 
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Figure 40. Suction Surface Stanton Number Variation With Reynolds Number 
at 10.5% Freestream Turbulence Intensity. 
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In comparison, Figs. 38, 39, and 40 show that a least squares curve 
fit of the Stanton number data follows approximately the same Re^ depen-* 
dency as the prediction (ST V * Re!! where n = -0.15 to -0.30). However, 

A A 

the overall level of the measured ST^ is 10% lower than predicted for a 
Tu ot of 6.6% and 20% lower, in the low Re^ range, for Tu ro of 8.0% and 
10.5%, With increasing Tu^, scatter in the data appears to increase. 
Within the limits of the data scatter, there appears to be no definite 
change in the local Stanton number level for a variation in Tu^ from 
6.6% to 10.5%. Since the Tu w level effect was not marked and higher 
Tu to levels indicated greater uncertainty in the heat flux measurement, 
the decision was made to conduct all subsequent film cooling experiments 
at the lowest value of Tu of 6.6%. 

CO 

V.2.3, Pressure Surface Results 

Figure 41 shows the data for the pressure surface, with the local 
Stanton number presented as a function of the local Reynolds number for 
the conditions listed in Table 5. The data encompasses all measurements 
of pressure surface dry-wall heat flux at all three Tu ot levels. The 
locations of the first three heat flux sensors are shown along the 
abscissa by their corresponding range of Re^, for the various total pres- 
sures imposed. 

The subject data exhibits the general trend of a decreasing Stanton 
number with increasing Reynolds number. In the region 7*10^ <. Re^ <. 1.2 x 

5 

10 , the Stanton number ranges from approximately 0.003 to 0.0023. A 

C 

step increase in the Stanton number occurs near Re^ - 1.25 x 10 with a 
subsequent general decrease with increasing Reynolds number. A slight 


:Ocal Stanton Number ST, 


OitwiO CTwoc^ Mode (2) 

*- ■ * — ■ x — XMode (l)r Re =300 
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\ A, 
v A 


A#° 
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N/m (psia) 
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A 2.43x10° (35) 

O 3.10xi0 5 (44.67) 

□ 3.82xl0 5 (55) 
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Figure 41. Pressure Surface Stanton Number Variation 
With Reynolds Number (Tu^ = 6. 6 10.5%). 
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depression in the Stanton number with an apparent recovery region is 
noted for the range 2.3 x 10 5 it Re^ < 6.2 x 10 5 . 

An attempt to predict the local Stanton number for the pressure 
surface was made using Mode (2) of the boundary layer program, with a 

c o 

representative total pressure of 3.1 x 10 N/m (50 psia) and negligible 
freestream turbulent kinetic energy. The resulting prediction is illus- 
trated in fig. 41 by the circled curve. The data follow the general 
trend of the prediction showing a decrease in STj, for increasing values 
of Re^, with the exception of the Reynolds number range 7 x 10^ ± Re^ <_ 

5 

1.2 x io . In this lower Reynolds number range the measured Stanton 
number values were nearly 50/5 lower than the Mode (2) prediction. If 
one were to extend the Mode (2) prediction past the region of apparent 
transition (-Re^ = 5.5 x 10^), it would appear that the data encompassed 
in the area marked “gage #1“ might be indicative of an extended laminar 

C 

boundary layer. The sudden rise in the measured ST^ at Re^ = 1.25 x 10 
may appear to be a subsequent transition to turbulent conditions. 

In an effort to account for the lower measured ST^ level in the 
"gage #1" region, Mode (1) of the prediction technique was used. Recall 
that Mode (1) corresponds to the solution of the conservation of mass, 
momentum and energy equations only. In this mode, the extent of the 
laminar boundary layer can be controlled by specifying a transition, 
momentum- thickness Reynolds number, Re 0 y fi . To predict what may appear 
to be measured, laminar ST^ values out to Re^ - 1.2 x io , a transition 
Reynolds number, Re 0 y^, of 300 was specified. The results of this 
prediction appear as the crossed curve. A comparison of the data from 
"gage #1" with the Mode (1) prediction shows that, at first, the 
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predicted trend of the Stanton number indicates that agreement with the 
data may be good. Note that up to Re^ - 5.5 x ]o 4 the slope of the Mode 

(1) ST^ vs Re^ is nearly identical to the trend of the measured data. 
However, for ftej, > 6 x 10 4 the difference between the predicted values 
and the data increases significantly. In fact, at Re^, = 9.6 x 10 4 
(which corresponds to Re^ ^ 200 along the surface) separation is indicated 
by the Mode (1) prediction as the calculated wall shear stress becomes 
vanishingly small. The measured ST^ values in the region of "gage #1" 
do not appear representative of extended laminar flow. From this infor- 
mation the following conclusions were drawn from which a subsequent 
analysis could be initiated to explain the trends of the data: 

(1) For Tu^ levels that were imposed during the measurement of the 
pressure surface ST X , it seems unlikely that a laminar boundary 
layer would persist out to Re^ = 1.2 x 10 5 . This is substan- 
tiated by the predicted transition Re^ of 6.9 x io 4 using Mode 
( 2 ), 

(2) By forcing laminar flow to a transition Re Q of 300, using 
Mode (1), there still appears to be a discrepancy between 
measurement and prediction, 

(3) If the data between 7 x ]o 4 < Re^ < 1.2 x io 5 is really char- 
acteristic of laminar flow then the data between 7 x 10 4 < Re^ 

< 2.2 x .0 represent a laminar transition to turbulent flow 
conditions. However, the Re^ range through which the data 
would indicate a transition is too abrupt, and 

(4) Laminar boundary layer separation is predicted at Re & = 200 
(Re^ = 9.6 x IQ 4 ). The data may actually be representative of 
separated flow. 
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From these observations. It appears most likely that separation occurred 
in a region near "gage #1". With this assumption, a separation bubble 
length model was used to predict reattachment as described in the follow- 
ing paragraphs. 

Horton (59) suggests that a separation bubble consists of two 
regions: (1) a laminar shear layer of length &-j, and (2) a reattachment 
zone of length & 2 * T,ie separation bubble is then the sum of the lengths 
of these two regions. The laminar shear layer length is described in a 
modified fashion by Dunham (60) as 


■ U .?IEJ. = 5000 v EXP (10.463 - 1.54 TuJ 5 ' 629 ) 


(73) 


The reattachment length, & 2 is defined according to Horton as 


O - 


RE 


*) 6c 


V 


u SE p' u sep 


u 


0.011233 (t#~) - 0,003033 
U SEP 


(74) 


where 


'RE 


SEP 


SEP 


= freestream velocity at reattachment 
= freestream velocity at separation 
= momentum thickness at separation 


For the conditions in Fig. 41, computations using Mode (1) yields, at 
separation 


U SEP " m /sec (340 ft/sec) 

v = 1.9 x 10"' 5 m 2 /sec (2.1 x 10*^ ft 2 /sec) 
0 SEp = 5,08 x 10" 5 m (0.002 in) 


Using a mean turbulence intensity of 8*02* assumed representative of the 
data, &i is calculated from Eg. (73) to be approximately 0.00612m (0.241 in). 


7 

* i 
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Before Jig can be calculated, the local velocity at reattachment 
must be estimated. Figure 42 shows the local Mach number variation with 
non-dimensional i zed distance along the pressure surface. The solid curve 
indicates the calculated design Mach number distribution. The locations 
of the heat flux gages are also indicated. The point of separation was 
predicted at x/L p = 0.125 (obtained from Mode (1) with P T = 3.47 * 10 5 
N/m (50 psia)). For the range 0.125 & x/L p 0.275* there appears to 
be little variation in the local velocity for the data shown. Thus, the 
value of U R g is not sharply defined. A value for the reattachment 
velocity of ~ °7.5 m/sec (320 ft/sec), corresponding to a local Mach 
number near 0,20 was selected as representative of the data in Fig. 42. 
Therefore, using the predicted momentum thickness and separation veloc- 
ity, Jig was calculated from Eq. (74) to be 4.06 x lO^ni (0.016 in). A 
total separation bubble length, (Jl-j + P„ 2 ), was then computed to be 
0.0066m (0.26 in), representing approximately 4% of the total pressure 
surface length. The selection of the reattachment velocity for use in 
the £ 2 calculation was somewhat arbitrary. However, since & 2 contributed 
less than 10% to the total bubble length, had little influence on the 
total length, Jl 1 + Jig. 

Erlich (61) presents another approximate method for the calculation 
of a separation bubble length based on flow visualization studies on an 
airfoil with a rapid change of curvature at the leading edge. Tests 
indicated that 

(P-j + Jlp) 

— n — = constant (75) 

U SEP 



Local Mach Number , Ma. 


P T = 2 .43-4 . 17x10 N/m {35-60Psia) 
T T ~ = 689°K (1240°R) 


Static Pressure 
Measurement Stations 




— »• * 9 — 

\a 


Calculated Design 
V“ Distribution 


% 

W 1 

\ 

q w 

2 

I 

0.1 

0.2 


Heat Flux Measuring 
Stations 


0.2 0.3 0.4 0.5 0.6 0.7 

Non dimensional Distance From Stagnation, x/L p 


Figure 42. Variation in Local Mach. Number At Various Pressure 
and Heat Flux Measurement Locations. 
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For conditions in which 200 <. Re Q S£ p < 500, this constant was found to 

* 

be 120. Recalling that G^p was 5.08 * 10 m the predicted bubble length 

-3 

from Eq. (75) is 6.1 x 10 m (0.24 in) or approximately 3.7% of the pres- 
sure surface length. This is in good agreement with the preceding 
results from Horton's analysis. 

It will be recalled from the discussion of Fig. 41, that separation 
was predicted at x/L p = 0.125 using Mode (1) with Re Q ~ 300. Adding the 
bubble length (&, + 18 Lp), ^attachment was predicted to occur 

at x/Lp = 0.165, or at a value of REATTACH = 1,27 * lo5 > which 13 32Ji 
greater than the Reynolds number at separation. Using Mode (2) of the 
prediction scheme, assuming the boundary layer reattaches at Re^ p^jYACH 
with a subsequent turbulent behavior, the local Stanton number variation 
was predicted with the result shown as the broken line in Fig. 41. The 
data appear to agree with the trends predicted oy the separation-reattach- 
ment model. The measured Stanton numbers are higher than predicted by 
laminar theory (Mode (1)) in the region of "gage #1", indicative of 
separated flow. Following reattachment, noted in the region marked 
"gage #2", the measured Stanton number level is representative of an 
established turbulent boundary layer. With the exception of the regime 
2.2 x io < Re^ < 7 x io , the measured Stanton number values appear to 
be in reasonable agreement with the results predicted for a reattached 
turbulent boundary layer. 

The trend of the experimental data for the region between 2,2 x io 3 
< Re^ < 7 x 10 , however, still needs to be qualified. It appears that, 
following reattachment, the measured Stanton number is initially depressed, 
but recovers to a level which is predicted for the reattached turbulent 
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boundary layer. The degree to which ST^ is depressed appears to be pres- 

d o 

sure dependent. At the lowest pressure level (2.43 x 10 N/m , 35 psia), 
indicated in Fig. 41 by triangles (a), the Stanton number is only slightly 
depressed, recovering to the predicted reattached turbulent boundary 
layer trend (broken line) near Re^ = 4.5 * 10 . At the highest pressure 

r O 

level (3.82 x 10 3 N/m , 55 psia), represented by squares (O), the largest 
depression of the Stanton number was noted. At this high pressure level, 

r r 

ST^ did not begin to recover until Re^ - 4 x 10 . At Re^ = 7 x 10 a 
full recovery was noted (i.e. a return to the trend of the broken line 
predicted by Mode (2) for reattachment). At an intermediate pressure 

c n 

(3.1 x to N/nr, 44.67 psia), shown as circles (o), depression and 
recovery of ST^ occurred approximately midway between the values for the 
two extreme pressures. The solid vertical ellipses correspond to pres- 

C O 

sures between 2.43 and 3.1 x 10 N/m (i.e. between the corresponding 
pressures of the triangles and circles) while the open ellipses represent 

c n 

pressures between 3.1 and 3,82 x io N/m (i.e. between the corresponding 
pressures of the circles and the squares). At first glance these results 
seem somewhat inconsistent with the reattachment prediction; however, the 
separation-reattachment model is capable of explaining this "depression- 
recovery" region of the local Stanton number by a close examination of 
the boundary layer after reattachment which follows below. 

Returning briefly to Fig. 42 the local Mach number is shown to 
remain relatively constant between x/Lp = 0.20 and 0.30. In fact, a 
slight deceleration may be indicated by the measured data even up to the 
location of the third heat flux sensor. A fairing through the wall 
pressure data indicates that the start of a favorable pressure gradient 
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may not occur until x/Lp * 0.325. The presence of freestream flow accel- 
eration or deceleration will certainly affect the development of the 
reattached boundary layer and, therein, provides an explanation for the 
trends observed in Fig. 41. 

C 

With reattachment predicted near Re^ - 1.27 x 10 , the measured 
values of ST X shown in Fig. 41 in the region marked "gage #2" are prob- 
ably representative of the reattached flow striving to establish itself 
as a normal turbulent boundary layer. However, at reattachment the 
boundary layer may be relatively thick as a result of the previous 
separation phenomenon. In such a case it would not be unreasonable to 
assume that the hydrodynamic boundary layer (characterized for this 
discussion by the momentum thickness, 8) is much thicker than the thermal 
boundary layer (characterized by the boundary layer enthalpy thickness, 
A). Since nr appreciable favorable pressure gradient is established 
until the near region of "gage #3", the reattached boundary layer cannot 
be thinned. In fact, it appears to thicken even more as evidence by the 
drop in ST^ between 2 x 10 < Re^ < 4 x 10 . 

As indicated from Fig. 42, a noticeable freestream acceleration 
begins in the region of "gage #3" (2 x 10 5 < Re x < 4 x lo 5 ). The onset 
of this acceleration directly affects the hydrodvnamic boundary layer 
development. Due to the assumed condition of the boundary layer at 
reattachment (i.e. 8 > A), this acceleration-dependent change of 0 in 
relation to A can offer a reasonable explanation for the pressure depen- 
dent "recovery" region. 

The effects of acceleration are characterized by the magnitude of 
the nondimens ional acceleration parameter, K. As the acceleration 


increases, 1C increases. Recall that K is defined as 


v dU 

k „ j? ® 

* " U Z dx 


( 76 ) 


In the region marked "gage #3", Fig. 4], the pressure-dependent recovery 

of STj- was noted, which corresponds to the initiation of local free- 

stream acceleration (Fig, 42). With increasing K, 9 thins and rapidly 

approaches the local value of A, This amounts to concentrating the 

effective transport properties of the momentum boundary layer within 

the limits of the thermal boundary layer. Consequently, the Stanton 

number will increase as 6 approaches A. The inverse situation, in which 

ST^r decreases for increasing K for A initially greater than e, has been 

experimentally verified by Kays, Moffat and Thielbahr (39), 

An explanation for the recovery of the Stanton number in the region 

following reattachment has been offered. However, the influence of the 

total pressure level on the initiation of recovery must be discussed. 

The data for the lowest pressure level (A symbols in Fig. 41) recovers 

sooner than the data for the higher pressure levels (o, □ symbols) . The 

delay of ST^ recovery with progressively higher pressure levels can be 

explained by referring to the definition of the acceleration parameter, K. 

As the pressure decreases, v OT increases with a subsequent increase in 1C. 

Consequently, at the same sensor location (as shown in Fig. 41 within the 

limits specified by "gage #3) any increase in the local Stanton number 

(i.e. recovery) due to acceleration-related thinning of the boundary 

layer should be noticed first at the lower pressure levels (i.e. higher 

local value of K). This is substantiated by the data (A symbols) between 
5 5 

2 x 10 < Rej, <7x10 in which ST^ at the lowest pressure level 


recovers quickly to near the predicted normal turbulent level* Data for 
the highest pressure level (□ symbol) shows ST^, recovery delayed until 
a much higher Reynolds number. All the data between the two extreme 
pressure levels are consistent with this interpretation. 

V.2.4. Summary 

The previous sections have been concerned with emphasizing two main 
points: (1) an environment can be maintained at reduced flow conditions 

which simulates the high heat flux levels characteristic of high temper- 
ature and pressure gas turbine operation and (2) the thin-foil, Gardon- 
gage, heat flux measurement technique incorporated in this investigation 
can provide reliable heat transfer data under a large gas-to-wall tem- 
perature ratio (Ty /T^ = 2.15) representative of the significant density 

CO 

variations through the boundary layer of a high temperature inlet tur- 
bine vane. Having both established the proper modeling environment and 
qualified the basic measurement technique, the film cooling study was 
initiated, the results of which are presented in the following section. 

V.3, Film Cooling Results 
V.3.1. Presentation of the Data 

The results of the subject film cooling study involving injection 
of cool air from a row of s panwise angled holes will now be presented. 

It is the sole purpose of this section to catalogue the recorded film 
cooling data in a convenient manner. Any comments concerning the data 
will be of a descriptive nature, restricted to pointing out general 
trends, noting various parameter levels of importance and making 
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appropriate comparisons. Discussion and analysis of the data will be 
thoroughly handled in the following section. 

The data to be presented represent local heat flux measurements 
which were made at various locations on the vane surface downstream from 
the row of injection holes. Each figure which will follow represents 

i 

the results from a particular heat flux sensor at a designated downstream 
position* (x-x^)/d Q . The data are presented in terms of the local 
Stanton number reduction, STp C /STg, and are shown as a function of the 
coolant blowing ratio, M, for each (x-x^/dg. Values of ST^/STg were 
calculated from heat flux measurements giving the local Stanton number 
with no film cooling (STg) and the Stanton number for specified film 
cooling conditions (specified value of M). 

Unless otherwise noted, all STp C /STg values which appear as points 
on the following figures are mean values representing heat flux measure- 
ments that were approximately t 2% to ± 6% from the mean. Larger devi- 
ations generally occurred near injection at high blowing ratios and 
these cases are documented with vertical bands designating the percent 
deviation from the mean. A few representative points in the lower blow- 
ing ratio regime are also banded showing the typical ± 2% to ± 6% 
variation. Figure 43 is a typical heat flux gage output visicorder 
trace for a sensor nearest the injection holes. A comparison is made 
showing the variation of the signal in relation to its mean for: (1) no 
blowing, (2) low blowing, M = 0.08, (3) moderate blowing, M = 0.84, and 
(4) strong blowing, M - 1.42. To eliminate extraneous data, the calcu- 
lated mean Stanton number ratio (computed from mean heat flux levels) 
had to meet a specified qualification criterion*, namely, the Stanton 


number without film cooling (STq, reference condition) did not change 
significantly during a series of film cooling measurements. To determine 
if this constraint was met, the value of ST Q was calculated from measure- 
ments before and after a range of coolant flow had been set and STp C 
measured. If STq had changed no tsjre than 2 to 3%, an average STq was 
used to form the ratio, ST fc /STq, for each coolant flow value. If STq 
changed more than 3 % during the experiments for a range of coolant flow, 
the data was discarded. The reader is referenced to Appendix C for a 
more complete discussion of the data qualification. 

Table 6 outlines the manner in which the data for the suction sur- 
face is organized herein. Six cooling configurations were investigated: 
two injection angles 3 = 18° and 3 = 35° with three locations (for the 
single row of holes) relative to the stagnation point of x.j/L s = 0.083, 
0.150 and 0.228. Cases I through III correspond to the results for 18° 
injection for the three injection locations shown, with heat flux gages 
located at nondlmensional distances from injection of (x-x^)/d Q as shown 
in Table 6. Cases IV through VI are for 35° injection at the same injec- 
tion locations (x^/L s ) as for 18° and for the same range of heat flux 
gage locations, (x-x.)/d Q . Figures for each value of (x-x.j)/d 0 are 
presented as each case number is dealt with separately. Table 7 is a 
summary of the injection configurations for the pressure surface. Again, 
six cooling configurations were investigated: two spanwise injection 

angles of 18° and 35° at three injection locations relative to the 
stagnation point of x^/Lp - 0.033, 0.121 and 0.193. Cases VII through 
IX represent results for 18° injection at the three injection locations 
for the range of heat flux gage locations (x-x^J/dg shown. Cases X 
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Case No. 

I 

II 

III 

IV 

V 

VI 

Injection Angle, 0 

18° 

18° 

18° 

35° 

35° 

35° 

Injection Distance 
From Stagnation 

x i/ L S 

0, 083 

0.150 

0.228 

0.083 

0.150 

0.228 


Dimensionless Distance From Injection (x - x^ 

>/ d 0 


10.45 

11.54 

16.29 

10.45 

11.54 

16.29 


22.20 

29.80 

28.24 

22.20 

29.80 

28.24 


42.55 

41.76 

41.07 

42.55 

41.76 

41.07 


53.50 

54.59 

53.03 

53.50 

66.54 

53.03 


66.33 

66.54 

100.85 

66.33 

90,46 

100.85 


78.23 

90.46 


78.28 

114.37 



102.20 

114.37 


102.20 




126.11 



126.11 




Table 6. Summary of Cooling Configurations For Suction Surface Injection 











Case No. 

VII 

VIII 

IX 

X 

XI 

XII 

Injection Angle, B 

18° 

18° 

18° 

35° 

35° 

35° 

Injection Distance 
From Stagnation 

V l p 

0.033 

0.121 

0.193 

0.033 

0.121 

0.193 


Dimensionless Distance From Injection (x - x^) /d^ 


10.35 

10.87 

17.17 

10.35 

10.87 

17.17 

*■ 

22.09 

27.17 

29.35 

22.09 

27.17 

29.35 


39.48 

39.35 

40.22 

39.48 

39.35 

40.22 


51.65 

50.22 

51.96 

51.65 

50.22 

51.96 


62.52 

61,96 

75.76 

62.52 

61.96 

75.76 


74.26 

85.76 


74.26 

85.76 



98.06 



98.06 




Table 7. Summary of Cooling Configurations For Pressure Surface Injection. 
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through XII are for 3 = 35° at identical values of x^/Lp and (x-x^)/d 0 
as in Cases VII to IX. 

A least-squares curve fit is included with the data for each heat 
flux gage. In some cases (for Use locations nearest injection), two 
least-squares curves were used to represent data within two distinct 
regimes. A smooth transition from one regime to the other was assumed 
consistent with the intersection of the two least-squares representations. 
A table will be presented summarizing all the least-squares equations 
representing the data. 

Before the data are presented it would be helpful to briefly review 
those parameters which were held fixed during the course of the entire 
film cooling investigation. An inlet Reynolds number (based on the vane 
leading edge radius) of 8.0 x 10^, an inlet Mach number of 0.22 and a 
freestream turbulence intensity of 6 , 6 % were maintained throughout the 
experiments, A gas-to-wall temperature ratio (Tj /T^) of 2.15 was 
maintained as was the prototype vane local Mach number distribution. 

The film coolant was supplied at ambient conditions and the coolant tem- 
perature, corrected for heat pick-up through the vane surface, resulted 
in an average dimensionless coolant temperature, g£, of 1.03. 

The data for Case I are presented in Figs. 44 to 51 corresponding 
to (x-x^J/dp values ranging from 10.45 to 126.11 for constant values of 
3 = 18° and x./Lg « 0.0833. Figure 44, at (x-x-J/d^ = 10.45, shows that 
STpq/STq decreases with increasing M up to M - 1.3. For M > 1.3, STp^/ST^ 
begins to increase with further increases in the blowing ratio. Near 
M = 1.3 large amplitude oscillations in the measured heat flux were noted 
(as shown with a band about the mean). As M was increased beyond a value 
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Figure 46. Film Cooling Stanton Nurrber Ratio Variation 
With M, Case I, (x-x. )/cL = 42.55. 
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CASE I 


Figure 

(x-x.)/a 0 

% Reduction in 
ST 0 

44 

10.45 

35 

45 

22.20 

C12 at M=2 . 15) 

46 

42.55 

12 

47 

53.50 

13.5 

48 

66.33 

10 

49 

78.28 

12.5 

50 

102.20 

7.5 

51 

126.11 

8.0 


1.3 

Not Defined 


% Deviation Prom 
Least-Squares 


+5 

+5.5 

+5 

+5 

+10 

+10 

+5 

+5 


Table 8. Summary of Details From Figures 44 to 51 
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of 1.3, the oscillations diminished in magnitude until at M ** 1.92 they 
represented only ± 2% of the mean STp C /ST Q . The minimum STp C /ST Q , 
occurring at M = 1,3, amounts to a 35% reduction in the local Stanton 
number from the non-film-cooled, dry-wall condition. However, for 
values up to M - 2.15, the value of STp C /ST Q was always less than 1.0, 
showing approximately a 9 % reduction in the heat flux for 1.8 < M < 2.15. 
The minimum Stanton number ratio which is observed for each (x-x^ )/d Q 
occurs at a particular value of the blowing ratio, designated hereafter 
as the "optimum blowing ratio", M^. The term M Q p T may be intro- 
duced in later discussions to describe the optimum blowing ratio found 
at the heat flux measurement location nearest to the coolant holes 
(i.e. smallest (x-x.)/d Q ). 

Two least-squares representations of the data are shown in Fig. 44 
as the solid line for the regions 0 < M < 1 .3 and 1 .3 < M < 2.14. Two 
curve fits were used to more accurately describe the data in each of 
these regions. Both curve fits represent the averaged data within ± 5%. 

At (x-x^ )/d Q = 22.20, Fig. 45 shows that STp C /ST Q decreases in a 
monatonic fashion with increasing M. The lowest Stanton number ratio 
occurred at M = 2.15, the highest blowing ratio investigated for Case I. 
The STp C /ST Q level out to M - 2.15 is higher than at (x-x^J/dp = 10.45 
for all M values being some 40% higher at M = 1.3 and 10% higher at 
M = 2.15, The least-squares curve fit represents the data within ± 5.5%. 

Similar results for (x-x^ )/d Q values from 42.55 to 126.11 are given 
in Figs. 46 to 51. In each of these figures, the value of STp C /ST Q 
decreases with increasing M to an optimum value of M(M«py) corresponding 
to the lowest recorded film cooled Stanton number. All the locations 
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from (x-x-)/dg = 42.55 to 126.11 exhibit a minimum STpg/STg at a value 
of M similar to that at (x-x^/dg = 10.45. Following Mg p y a subsequent 
rise in STpg/STg occurs with increasing M but never exceeds the value of 
1.0 for values of M up to 2.15. Table 8 presents a summary of the sig- 
nificant details of Figs. 44 to 51 indicating the value of M at which 
film cooling was most effective (Mg pT ) and the percent reduction from 
STg at Mgpj. Between (x~x.j)/dg of 42.55 and 126.11 minimum values of 
STfc/STq at Mgpy corresponded to a 7.5% to 13.5% reduction from STg. A 
general trend of lower percentage reductions from ST Q with increased 
distance from injection was noted. Least-squares representations of the 
data were within ± 5% or ± 10%. 

The results for Case II are shown in Figs. 52 to 58, representing 
heat flux measurements at (x-x-)/dg from 11.54 to 114.37 for film coolant 
injection at x./L s = 0.150 with 8 = 18° (see Table 6). Figure 52 shows 
STpg/STg values as a function of M at (x-x.)/dg = 11.54. The Stanton 
number ratio decreases with increasing M up to M s 1 .05* thus defining 
the value of Mg p y at which STpg/STg was a minimum. This value of Mg p y 
was established by obtaining least-squares curve fits of the data using 
two separate linear functions. A smooth transition between the two 
functions was constructed consistent with the trend established by the 
data to give Mg p y » 1.05. In the region M > 1.05, STp C /ST Q increases, 
reaching a level of 0.97 at M - 2.37. Within this range of M, mean 
values of STpg/STg represent actual levels that are + 9% of the mean. 

The minimum value of STpg/STg, occurring at Mg p y - 1.05, amounts to a 
38% reduction from STg due to film cooling. The least-squares curves 
represent the average data within ± 10%. 
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CASE II 


Figure 

(x-*.)/d 0 

% Reduction in 
ST q 

M 

OPT 

% Deviation From 
Least-Squares 

52 

11.54 

38 

1.0 

+10 

53 

29.80 

11.5 

1.2-1. 5 

+10 

54 

41.76 

10.5 

1.2-1. 5 

+5 

55 

54.59 

10.0 

1.0-1. 2 

+10 

56 

66.54 

8.5 

1. 0-1.2 

+10 

57 

90.46 

9.0 

1. 1-1,3 

+10 

58 

114.37 

9.5 

1.4-1. 6 

+5 


Table 9 . Summary of Details From Figures 52 to 58 
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Table 9 lists the percentage reduction from ST Q at each Mq PT for 
(x-x^/dg between 29.8 and 114.37 with the results for these remaining 
downstream locations of Case II presented in Figs. 53 to 58. For these 
figures STp C /STg decreases to a minimum value for increasing M up to Mg pT . 
M opT for these locations has been defined within a range rather than a 
specific value due to the gradual trend of the data. For (x-x.*)/dg 
increasing from 11.54 the percentage reduction of ST Q due to film cooling 
at M = Mg P j decreases from 38% to approximately 9%. The range of values 
of M QpT noted for (x-x^J/dg * 29 to 114 bracket the value of Mg p ^ deter- 
mined for the location nearest injection (Mg P j WEAR ). The range of M QpT 
for each (x-x^J/dg > 11,54 was determined from the least-squares curve 
fit which represented average STpg/STg values within ± 5% and ± 10%. 

The data for Case III with injection at x^/L^ = 0.228 and p = 18° 
are presented in Figs. 59 through 63 with a summary of the details listed 
in Table 10. Values of STpg/STg are shown as a function of M for 
(x-x-)/dg from 16.29 to 100.85. Figure 59 shows that at (x-x^ )/d Q = 16.29 
the Stanton number ratio decreases with increasing M up to Mg p y = 0.50. 

At this value of the blowing ratio the value of STp C was 16% lower than 
STg. For M > Mg P j» STp C /STg increases with M and exceeds unity at an M 
value of approximately 1.55. Thus, for M > 1.55, the effects of film 
coolin'] were to increase the heat flux. As the blowing ratio exceeded 
MgpT = 0.5, the measured heat flux began to oscillate with peak-to-peak 
levels which were larger than usual. The magnitude of these fluctuations 
is indicated in relation to the mean value of STpg/STg by the vertical 
bands. The oscillations caused STpg/STg to range from ± 7.5% to ± 12.5% 


of the mean. 





Figure 60. Film Cooling Stantcn Nunfcer Ratio Variaticn 
With M, Case HI, = 28.24, 
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Figure 62. Filin Cooling Stanton Nunfcer Ratio Variation 
With M, Case III, (x-x^/dg = 53.03. 
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CASE III 


Figure 

(x-x. )/a n 

% Reduction in 

M 

OPT 

% Deviation From 



STq 

Least-Squares 

59 

16.29 

16 

0.50 

+5 MOI 
+25 H)M““ 

60 

28.24 

7 

0. 8-1.0 

+11 

61 

41.07 

4.5 

0.50 

+6 

62 

53.03 

8 

0. 8-1.0 

+5 

63 

100.85 

7.5 

0. 8-1.0 

+2 


Table 10. Summary of Details From Figures 59 to 63 
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Figure 60, at (x-x^ )/dg = 23.24, shows STpg/STg decreases up to M 
between 0.8 and 1.0. In this narrow M range, film cooling reduced the 
Stanton number 7 %. For M > 1.0, STp^ increases to within 2 % of the value 
with no cooling at M = 1.67. 

Farther downstream at (x-x.j)/dg = 41,07, Fig. 61 shows a greater 
dependency of STpg/STg on M in the higher blowing ratio regions (M > 0.5) 
than at the previous location. A slight decrease in STp C /STg was noted 
for M increasing to 0.5 at which point the ratio began to increase. At 
M = 0,5 a 4.5/5 reduction from ST Q was measured. The rise in the Stanton 
number ratio for M increasing beyond 0.5 is rapid, yielding a value of 
1.11 at M = 1.67. 

The two measurement locations farthest from injection are presented 
in Figs. 62 and 63. At both locations a gradual STp C /STg decrease with 
increasing M is shown out to a value near M = 1.0. At the highest M 
investigated (M = 1.67) a 5 % reduction from STg was measured at (x-x.)/dg 
= 53.03. At the same value of M, a 7 % reduction from STg was noted for 
(x-x^/dg = 100.85. 

Similar data with an injection angle of 3 - 35° were obtained for 
“ 0.0833, 0.150 and 0.228. The data for Case IV with 8 = 35° and 
x.j/L s = 0.0833 are presented in Figs. 64 to 71 for ( x-x^ )/d Q from 10.45 
to 126.11. It should be noted that the maximum value of M as shown in 
these figures varies from 0.88 to 2.30. The criterion for data quali- 
fication (see Appendix C) was not met in certain ranges of M for some 
heat flux gage locations and consequently the data for those ranges were 
discarded. 
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Figure 64. Film Cooling Stanton Nirrber Patio Variation 
With M r Case IV, (x-XjJ/<3q = 10.45. 



Figure 65, Film Cooling Stantcn Nurrber Ratio Variation 
Vlith M, Case IV, (k-Xj. )/d^ =* 22.20. 
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Figure 64, at (x-x^)/^ = 1°*45, shows a gradual decrease in ST FC /ST Q 
for increasing M with no clear Hq PT apparent. Up to M * 2.3, the local 
Stanton number was reduced by 1355 from ST Q . All the data is represented 
by a least-squares curve within ± 555. 

Similar results were found for Figs. 65 through 71 with Table 11 
summarizing the important features for (x-x.)/dg from 10.45 to 126.11. 

Up to the maximum M investigated, ST FC /ST Q shov/ed a gradual decrease, 
with absolute levels dependent upon the value of (x-x..)/dg. Consequently, 
the maximum value of M is reported since no clear Mgp^ was observed. For 
those locations in which ST FC /STq values were qualified up to M = 1.78, 

the Stanton number ratio varied from 0.93 at (x-x.)/dg = 22.2 and 66,33 

to values near 0.96 and 0.95 for (x-x.)/d 0 = 102.20 and 126.11. 

The results for Case V at 8 - 35° and x./L<. = 0.150 are illustrated 
in Figs. 72 through 77. Figure 72 shows a reduction in the Stanton 
number ratio for increasing M at (x-x..)/dg = 11.54. The minimum ST fc /STq 
value occurred at M QpT = 0.85. At M opT , a 32 % reduction from ST Q was 
measured. For M > M opT , ST FC /STg rises to within 9 % of ST Q at M = 2.1. 

Two linear, least-squares representations of the data are illustrated 
with the data + 555 about the line for M < 0.0 and ± 1055 for M > 0.8. 

Table 12 summarizes the results presented in Figs. 72 to 77. The 

figures show a general trend of a decreasing ST fc /STq with increasing M 
for all locations greater than (x-x^J/dg = 11.54 with no obvious M™. 

The lowest STj-g/STg occurs near the largest value of M investigated. 

At M = 2.1 the percent reduction from STg due to film cooling ranged 
from 1855 to 7.555, decreasing as (x-x^/dg increases. 
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Measurements of the local heat flux were taken with and without 
film cooling for Case VI of the suction surface and Cases VII through 
XII for the pressure surface. The variation of STq for these cases 
during the course of investigating a range of coolant conditions was 
such that the data qualification criterion was not satisfied. Conse- 
quently, no film cooling results are presented. The reader is referred 
to Appendix C for further discussion concerning data qualification. 

Table 13 lists all of the least-squares, curve-fit equations used 
to represent the data of Figures 44 to 77. These equations will be used 
in a subsequent section for further analysis of the data. 

V.3.2. Analysis and Discussion of Film Cooling Results 

The previous section presented the results of the study of film 
cooling of a turbine vane leading edge for three different locations of 
coolant injection and two injection angles (see Figs. 44 to 77). These 
results will now be analyzed and discussed to isolate the effects on 
STpq/STq due to changes in three primary cooling parameters: (a) the 

blowing ratio, M» (b) the angle of injection, 0, and (c) the location 
of injection relative to the stagnation point, x^/Lj. It will L-- shown 
that STpq/STq decreases with increasing M for all the subject configu- 
rations (i.e. Cases I - V). However, the value of M at which a minimum 
value of STp C /STg occurs (most effective cooling) depends upon the injec- 
tion angle and the location of injection (indicative of the thickness of 
the boundary layer at injection). The minimum level of STp C /ST 0 depends 
upon the location of injection (which characterizes the local degree of 
freestream acceleration (K.)). It is postulated, judging from an 
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CASE I 



FIGURE 

A 


B 

C 

44 

“rfjoPT 



-0.298 

1.009 

44 

-0.509 


1.960 

-1.072 

45 

OPT 

0.013 


-0.061 

0.950 

46 


0.074 


-0.178 

0.995 

47 


0.079 


-0.206 

0.998 

48 


0.060 


-0.174 

1.023 

49 


0.105 


-0.271 

1.047 

50 


0.039 


-0.092 

0,997 

51 


0.064 


-0.163 

1.021 



CASE II 



52 

h<m opt 

m>m opt 



-0.426 

1.007 

52 



0.287 

0.310 

53 

0.048 


-0.134 

0.979 

54 


0*046 


-0.121 

0.974 

55 


0.098 


-0.193 

0.993 

56 


0.045 


-0.096 

0.967 

57 


0.047 


-0.118 

0.984 

58 


0.036 


-0.113 

0.992 



CASE III 


59 

M<M™ 



-0.355 

1.000 

59 


0.083 


-0.054 

0.881 

60 

0.070 


-0.123 

0.987 

61 


0.115 


-0.118 

0.995 

62 


0.070 


-0.131 

0.981 

63 


0.068 


-0.144 

1.001 



CASE IV 



64 


0.036 


-0.108 

0.950 

65 




-0.0021 

0.973 

66 




-0.029 

0.980 

67 




-0.050 

0.993 

68 




-0.039 

1.000 

69 




-0.027 

0.993 

70 




-0.0075 

0.979 

71 



CASE V 

-0.025 

0.990 

72 

m>m opJ 



-0.389 

0.965 

72 

73 



0.174 

-0.054 

0.535 

0,933 

74 




-0.043 

0.960 

75 




-0.064 

0.953 

76 




-0.009 

0.949 

77 




-0.017 

0.953 


Table 13 
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analysis of the present data, that large velocity gradients at injection 
(K. ~7 x 10 ) can cause a significant deterioration in the cooling 
effectiveness in the form of enhanced turbulent mixing at the outer edges 
of the boundary layer. Specifying shallower injection angles in regions 
of severe acceleration can minimize the interaction between the coolant 
film and the mainstream flow by maintaining the low jet trajectories 
necessary to avoid penetration into the outer regions of the boundary 
layer. 

V.3,2,1, Effect of Blowing Ratio, M 

To evaluate the effect of the blowing ratio on STp C /ST Q , the results 
of Figs. 44 to 77 will be presented in a manner which can illustrate the 
influence of increasing coolant mass addition at all locations downstream 
from injection. The figures presented in this section were generated 
from the least-squares curves representing the data in Figs. 44 to 77. 

The method for constructing the figures in this discussion is out- 
lined in Fig. 78. For a particular configuration (Cases I - V) a specif- 
ic value of M was selected, H-j . Figures 73(a), (b) and (c) show a 
typical least-squares representation of the variation of STp C /ST Q with 
M for three locations from injection (same form as Figs. 44 - 77). For 
Ml the values of STp C /ST 0 at x-j, x 2 , and x 3 were noted. Figure 78(d) 
shows the resulting transformation depicting the STp C /STg variation with 
(x-x.j)/clQ at M - M-j. Similarly, the procedure is repeated at any value 
of M, resulting in a family of curves which illustrates the overall 
effect of M and (x-x^J/dg on the film cooling performance. 

Since the data will be represented by the least-squares values of 
STpg/STg for each M» the smooth curves generated (i.e. Fig. 78(d)) will 
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be a function of the least-squares equations. The measured STp C /STg 
values were shown in Figs. 44 - 77 to exhibit a definite amount of scatter 
in varying degrees; therefore, the least-squares curves will reflect this 
uncertainty to some extent. Because of this transformation procedure, 
the figures generated will show a certain amount of waviness (as illus- 
trated in Fig. 78(e)). Generally, the wavi ness occurs in regions far 
from injection (40 or 50 hole diameters downstream) where the cooling 
effect has decayed to such an extent that uncertainties in the data are 
approaching the level of the cooling effect. The waviness should not 
be given any more significance than the fact that the STpg/STg level was 
quite high (i.e, cooling effectiveness fairly low). Rather than curve 
fit the transformed data (essentially a second curve fit), it was 
decided to maintain the characteristic transformed waviness and simply 
draw a smooth faring through the points for selected values of M. 

The question may arise as to what effect the nature of any lateral 
film nonuniformity may have on the data. It would seem plausible that 
such nonuniformities may tend to produce the waviness previously mentioned. 
It was not possible to experimentally determine the lateral variation 
(across the span, i.e, hole-to-hole), of STp C /ST 0 at any (x-x.)/cIq. 

However, it is not expected to be sufficient to account for the wavy 
nature of the transformed data. Recent investigations at the NASA Lewis 
Research Center (62) indicate that for a similar geometry as the subject 
investigation, the cooling effectiveness at (x-x^J/dg =11 varied less 
than 15% across the span. As (x-x^/rig increased, the cooling effective- 
ness was reduced but became more uniform. 
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The effect of the blowing ratio on STp C /ST Q is illustrated in Figs. 
79 through 87. This series of figures is divided into two groups: (1) 

injection angle of T8°, and (2) injection angle of 35°. Within each 
group the effects of M on STp C /ST Q for the various injection locations 
will be compared. For each configuration within a group, two figures 
will generally be presented: the first figure describing the effect of 

increasing M up to the value corresponding to a minimum ST fc /STq at 
(x-x^/dg nearest injection (Mgpy and the second characterizing 
the dependency of STp C /ST Q on M for blowing ratios greater than M Q p T 
Presentation of the data in this manner allows for the identification of 
four main features: 

(1) the magnitude of the blowing ratio effect on STp C /ST Q in the 
near injection region, (x-*x..)/d 0 < 30, 

(2) the degree of persistence of the cool film for various M 
values by observing the decay of STp c /ST Q in the near region, 

(3) the shifting of the overall level of STp C /ST Q to a lower 

value as M increases to M QpT and 

(4) the shifting of the overall level of STpg/STg when M surpasses 

M 0PT,NEAR’ 

Figures 79 and 80 describe the STp C /ST Q variation with (x-x^/dg 
for {3 » 18° and injection at x^/Lg a 0.0833 (Case I). Figure 79 shows 
a family of curves for M ranging from 0.2 to M QpT NEAR = 1.3. Figure 80 
includes the higher blowing ratio range for M from 1.3 to 2.0. 

In Fig. 79, the most significant reduction in STpg/STg for an 
increase in M from 0.2 to 1.3 occurs in a region 0 < (x-x^J/dg < 25, 

The lowest Stanton number recorded occurred at M = 1.3, amounting to a 



Figure 79- Effect of Blowing Ratio on Film Cooling Stanton Number 
Ratio For s = 18°, X./L s = 0.083 (M < >t 0 p T/SEaE > . 
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35% reduction in STq. Near M - 1.3 the deterioration of the cooling 
effect downstream is severe showing that most of the benefit of increas- 
ing M is lost after approximately 25 hole diameters downstream (as shown 
by the close proximity of the curves for M = 0,8, 1.0 and 1.3). However, 
farther downstream, for 50 < (x-x..)/dQ < 90, there appears to be an over- 
all, beneficial shifting of ST^/STq to lower values with increasing M, 
amounting to approximately an 11% decrease in the Stanton number ratio 
from M = 0.2 to M = 1 .3. 

A reversal of the trend shown in Fig. 79 occurs with M > 1,3. 

Figure 80 illustrates that for (x-x.)/d Q < 20, STp C /STg increases as M 
is increased from 1.3 to 2.0. Between 20 < (x-x.)/d 0 < 30, however, the 
Stanton number ratio remains nearly constant. For (x-x..)/d 0 > 30, 
ST fc /ST 0 increases with M, but to a lesser degree than near injection 
for the same change in M. A blowing ratio of M = 1.3 appears to be the 
value at which STpg/ST 0 is lowest for all downstream locations. 

Figures 81 and 82 illustrate the blowing ratio effect on STp C /ST Q 
for 5 = 18° at x^/Lg - 0.150 (Case II), Since was identified 

as 1.05 (see Fig. 52) for this case, Fig. 81 presents a family of curves 
for M = 0,2 to M = 1.05. Figure 82 includes blowing ratios from 1.05 
to 1.8. 

The region which experiences the greatest benefit from increasing 
M appears to be concentrated, as shown in Fig. 81, in the range 
(x-x^/dg ^ 30. At M = 1.05, the lowest value of STp^/STg amounted to a 
38% decrease from ST Q at (x-x.J/dg = 11.54. For (x-x.)/d 0 > 30, STp C /ST 0 
decreases by approximately 5% for an increase of M from 0.2 to 1.05. 



Figure 81, 


Effect of Blowing Ratio on Film Cooling Stanton Number 
Ratio, 3 = 18% Xi /L s = 0.150 {M < M QpTfNEAR ). 
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Figure 82 shows a reversal of the effects of increasing M past 
Mqpj similar to that noted in the previous configuration (Case I), 
Near (x-x^O/dg - 11.54, the value of STpg/STg increases with increasing 
M. At M = 1,8, however, the local Stanton number with film cooling was 
still 17.5% lower than STg. For (x-x^J/dg > 30, STp C /ST Q seems insen- 
sitive to increases in M from 1.05 to 1.8. The overall level of STpg/STg 
stayed approximately constant for (x-x^)/dg > 30. (No significance is 
assigned to the wave near (x~x.)/dg - 55 for K - 1.8.) It appears that 
the optimum blowing ratio for this configuration is near 1,05, 

The results for Case III for 8 = 18° and x.j/L s = 0.228 appear in 
Figs. 83 and 84 wherein Mgp^ has been identified as approximately 
0.5. Figure 83 presents the Stanton number ratio as a function of 
(x-x.)/dg for the range 0.20 < M < 0.50. Figure 84 includes a range of 
M from 0.50 to 1.3. An overall decrease in STpg/STg is noted in Fig. 83 
with an increase in M up to a value of 0.5. Near injection at (x-x^/dg 
c 16.29, a 16% decrease in STp C from STg was measured at M - 0.5. For 
(x-x.)/d Q > 40 and M “ 0.5, the reduction in STp C from STg was approx- 
imately 5%. 

Figure 84 shows the variation of STp C /STg for M > 0.5. The blowing 
ratio appears to affect the Stanton number ratio differently for three 
distinct regions: (1) (x-x^)/dg < 30 in which STp C /STg increases with 
increasing M, (2) 30 < (x-x^J/dg < 60 in which a substantial "peaking" 
and subsequent decrease of STp C /STg is shown as a function of increasing 
M, and (3) 60 < (x~x^)/dg < 100 in which the Stanton number ratio is 
reduced slightly ("3%) over the value at M a 0.5 but remains essentially 
constant for M increasing from 0.8 to 1.3. In the initial region, the 
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local Stanton number at M — 1.3 was still 555 lower than that without 
film cooling. For the second region, the local Stanton number at M = 1.3 
increased to a level 42 greater than STg, but then fell to a value 
nearly TOSS lower than STg. In the final region, nearly a constant value 
of ST^/STq is illustrated, representing a 7 - 102 decrease from STg 
for M between 0.8 and 1.3, 

Figures 85 through 87 illustrate the effective Stanton number reduc- 
tion with increased blowing for Cases IV and V with injection at 3 = 35°. 
Figure 85 gives results for x^/L^ = 0.0833, and Figs. 86 and 87 describe 
injection at x../Lg = 0.150. 

Figure 85 describes the Stanton number reduction due to increased 
blowing which, in contrast to the previous configurations, exhibited no 
optimum M for blowing up to M = 1.8. Near injection ((x-x_j)/dg « 10) 
the Stanton number ratio decreases with increased blowing. An increase 
in M from 0.2 to 1.8 only results in a reduction of STp C from STq of 72 
to 132, respectively. In the region 35 < (x-x.j)/dg < 70, for the same 
increase in M, STp C /ST Q is reduced, but only in the amount of 22 to 42. 

For Case V, shown in Figs. 86 and 87, injection takes place at 
x^/Lg = 0.150 for 6 = 35°. Figure 86 illustrates the reduction of 
STpQ/STg due to increased blowing from M - 0.2 to M = 0.85. An optimum 
blowing ratio of 0.85 near injection {{x-x^/dg ~ 11) is noted. The 
greatest benefit of increasing the coolant flow is realized in the 
region (x-x^J/dg < 30, with a 322 reduction from STq at {x-x^J/dg - 11.54 
for M = 0.85. The cooling capacity of the film decays rapidly near 
injection showing only a 7.52 reduction from ST Q at (x-x^)/d Q a 50 for 
M = 0.85. For 90 < (x-x^/dg < 114, the value of STp C converges to 








Figure 86. Effect of Blowing Ratio on Film Cooling Stanton Number 
Ratio, 8 = 35”, Xi /L s - 0.150 CM < M opTjHEM? ) . 






Figure 87 . Effect of Blowing Ratio on Film Cooling Stanton Number 
Ratio, 8—35 , — 0.150 (M > ^qft,NEAR^ • 
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approximately a 5 . 5% reduction from STg for all blowing ratios up to 
M = 0.85, 


Figure 87 includes the higher blowing ratio range for Case V (0.85 < 
t<i < 1.8). Three regions can be clearly distinguished from this figure: 

(1) (x-x^/dg < 20 in which STj- c increases for M > 0.85, (2) 20 < (x-x^/dQ 
< 80 showing a decrease in ST^/STq (improved cooling) for M > 0.85 and 
(3) (x-x*}/dQ > 80 in which STp C /STg appears insensitive to changes in the 
blowing ratio. 

The significant features indicated by the data presented in Figs. 79 
to 87 may be summarized as follows: 

(1) For the shallow injection angle (3 “ 18°), a systematic de- 
crease in the value of STp C /ST Q for increasing M up to M Q p T WEAR 
was noted at all locations downstream from injection, while for 
M > Mgpj an increase in STp C /STg with increasing M was 
observed. The value of Mg p ^. increased as the injection 
location was moved closer to stagnation. 

(2) For the steep injection angle (3 = 35°) the trends varied with 
the location of injection relative to stagnation. For injection 
at = 0.083 (Fig. 85), the data exhibited no optimum M 
value. The value of STp C /ST 0 decreased continually as M in- 
creased for all locations downstream from injection. For injec- 
tion at x.j/L s = 0.150 a value for Mg p j was defined. How- 
ever, only near injection, (x-x^/dg £ 20, did the trend of 
STp C /STg vs M follow that outlined in (1) above. Further down- 
stream, the value of STp c /STg decreased progressively as M was 
increased. 

The data for 3 - 18°, for all three injection locations, follow the 
trends expected for STpg/STg vs M. As M is increased, the film cooling 


performance increases (i.e, lower values of STp^/STg) due to the increased 
thermal capacity of the coolant flow. However, at some condition, further 
increases in M enhance the mixing between the coolant and the mainstream 
such that mixing effects predominate over thermal capacity effects. There- 
after, increasing the value of M results in a decrease in the cooling per- 
formance (i.e. STp C /STQ increases). The data for 3=35°, x./Lg = 0.150, 
near injection, are also consistent with the aforementioned trenss. How- 
ever, for the most part, injection at 3=35° resulted in some decrease in 
the value of STp C /ST 0 as M was continually increased. The change in the 
level of STpq/STq was relatively insensitive to M for M values greater than 
approximately 1.0. 

V.3.2.2. Correlation of STp^/STg With Two-Dimensional Model 

In order to make an assessment of the relative effect on STp^/STg of 
geometrical configuration changes, such as the injection angle (3) and the 
injection location (x^/L s ), the data were correlated with the blowing ratio 
(M) and the dimensionless distance from injection in a manner suggested by 
a two-dimensional slot cooling model. By placing each configuration on a 
common basis in relation to M and the dimensionless distance from injection 
the effects on STpg/STg for changes in 3 from 18° to 35° and x^/Lg fv'om 0.0833 
to 0.150 can be illustrated (Case III, 3 = 18° and x^/Lg = 0.228 will not 
be included in the comparisons due to the uncertainty in M opT m EAR )* 

As the spacing between coolant holes for a single row injection con- 
figuration approaches zero, two-dimensional slot-type cooling is achieved 
in the limit. Since models for three-dimensional (hole) injection are 
few (Eckert) proposes one based on a point heat source moving in a semi- 
infinite medium which requires knowledge of the eddy diffusivity) , it 
is convenient to correlate the data as if the film coolant was issuing 
from a slot of total area and span-width equal to that of the row of 
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holes. Since lateral variations in STp C /ST Q are expected to be minimal 
for the spanwise injection angles of the subject investigation, use of 
the slot cooling correlating parameters seemed justified. 

A two-dimensional model, correlating the effects of the blowing 
ratio and distance from injection, was developed by Klein and Tribus 
(63). Assuming a line-source of enthalpy input, the adiabatic cooling 
effectiveness was expressed 



U °- Z C D X-X "°- 8 

» 4.62 Pr 2 ' 3 Ro c °- 2 $) (/') 

Pco 

(78) 

where 

Re^ a coolant Reynolds number 


Defining 

X = ( Re c -°- 25 


(79) 

then 

VL 0,2 C p r “0.8 
^ADW “ 4 ' 62 fr ' * 

(80) 


Pc 


In the present investigation, the values of Pr, u /iv and C n /(L 

u P c p c3 

were essentially constant. Therefore, to correlate the effects of (x-x^) 
and M variations, the data are presented as a function of the parameter A. 

Figures 88 through 92 present STpg/STg data for Cases I - V corre- 
lated with A. The data shown in each figure represent STpg/STg values 
for M < Mgp T The effects of boundary layer injection-related tur- 

bulent mixing possibly associated with values of M > Mgp-j. would not 
be expected to correlate with the two-dimensional model. Consequently, 
the data for M > ^oPT a NEAR are not inc ^ uded in the correlation at this 
point, but will be presented in later figures for comparison purposes. 
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FIGUHC 

CASE 

M K2NGE 

IEaST-SQOMES EQUKTCCNS 

% IEVIMTCN FBOM 
IEaST-SQUMES 

88 

I 

0 - 1.3 

-0.096 Dog X) 4 + 0.738 (log A) 3 - 2.058(log A) 2 
cm + 2.525(log A) - 0.254 

0 

±8 

89 

II 

0 — 1 » 05 

ST 

^ = “0.083 (log A) 4 + 0.720 (log A} 3 - 2. 318 (log A) 2 
S2? 0 + 3. 322 (log A) - 0.869 

+8 

90 

III 

0 - 0.5 

STpc 

__ = “0.050 (log A) 2 + 0.275 (log A) + 0.609 

STg 

±6.5 

91 

IV 

0 - 1.8 

ST FC 2 

= “0.016 (log A r + 0.100 (log A) + 0.840 

ST 0 

±5 

92 

V 

0 - 0.85 

OT FC a 3 2 

« “0.032 (log A) 4 + 0.344 (log A) - 1.339 (log Ar 

ST Q + 2.285 (log A) - 0.518 

±6 


Table 14, Suamary of Details Fran Figures 88 - 92. 
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Figure 88 , 

Correlation of Measured ST^/STq For 8 = 18° 
and 2?j_/Lg = 0.083 (M < M o?0? r NEAR^ * 
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Figure 09. 

Correlation of Measured ST FC /ST 0 For 0 
and Xj,/Lg - 0.150 (M < M 0PT/MEAR^ * 


Stanton Number Reduction, ST^/ST 
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Figure 90, 


of Measured ST FC /ST 0 For 3 
0,228 (M < M 0 p TfNEAR ) * 
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= 18° 


Stanton Number Reduction, ST^/ST 
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Figure 91 . 

Correlation of Measured ST-^/ST- For 

FC' 0 

» = 3S ° and X i/ L s = 0.083 {M < M f!AX ) . 



Stanton Number Reduction, ST„-/ST 
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Figure 92. 

Correlation of Measured ST^/STq For 3 = 35° 
and X-j/Lg = 0.150 (M < m opT/NEAR ) * 
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Figure 92. 
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Table 14 summarizes the pertinent details for the data correlations 
of Figs. 88 to 92. A least-squares curve fit of the data was developed 
for each of the figures. A fourth-order polynomial curve fit is shown 
for the data of Figs. 88 and 89, representing the data within ± 8% in 
both cases. Cases III and IV are illustrated in Figs. 90 and 91, showing 
quadratic, least-squares curves that represent the data within ± 6.5% and 
± 5%, respectively. The data of Case V are within ± 8% of a fourth-order 
polynomial , least-squares curve. The data are now in a convenient form 
in which the effects of changing various injection parameters other than 
M and (x-x^) can be illustrated. 


V.3.2.3. Effect of Injection Angle, & 

By comparing the least-squares representations of STp C /ST Q as a 
function of X for M < M QPT near > the effect of the injection angle, B, 
on the Stanton number ratio can be assessed at both x./L s = 0.0833 and 
x./L-e = 0.150. The reduction of the local Stanton number due to film 

T o 

coolant injection at x./L s “ 0.0833 is compared in Fig. 93 for injection 
angles of g - 35° and 8 = 18®. The broken line is the least-squares 
curve fit for the data at the 35° injection angle for 0 < M < 1.8 (this 
is Case IV in which M QpT NEAR « shown in Fig. 91). The solid line 
(least-squares curve from Fig. 88) represents injection at an angle of 
18° for 0 < M < 1.3. The open, circular symbols designate the Stanton 
number ratio for values of the blowing ratio greater than m qPT*NEAR ^ or 
8 = 18°; they are not well correlated with X but are included for com- 
parison and completeness. 

The predominate feature which is illustrated in Fig. 93 is that the 
shallower injection angle (8 - 18°) is more effective than the 35° 


Stanton Number Reduction, ST„ r /ST 



Figure 93, 


Effect of Coolant Injection Angle For 
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0,083. 
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injection angle in lowering STp^/STg. Near X - 4, the 18° injection 
angle shows a Stanton number ratio of 0.670 while a value of 0.895 is 
shown at 3 = 35°, a difference of 25%. As the blowing ratio is decreased 
(M < 1.3) and/or the distance from injection increases, X increases, and 
in the range 10 < x < 600 the difference between STp C /ST 0 for the two 
angles narrows considerably. At high blowing ratios (M > 1.3), as rep- 
resented by the symbols for 3 = 18°, a shaded region is shown indicating 
the data for (x-x.)/d 0 ~ 11. The effects on STp(,/ST 0 of large scale 
mixing resulting from M increasing past 1.3 are characterized by the 
direction of the large vertical arrow. For M > 1.3, with increasing 
(x-x-)/d,-j (shown by direction of large horizontal arrow), STp C /ST Q data 
is quite scattered when correlated with X. 

Figure 94 compares the data for 18° and 35° injection for x^/Lg = 
0,150. STp C /ST 0 is shown as a function of X for the too least-squares 
representations of the data. The broken line corresponds to 3 - 35° for 
0 < M < 0.85 (from Fig. 92). The solid line represents 8 = 18° for 
0 < M < 1.05 (from Fig. 89). The circular symbols represent STp C /ST Q 
values for M > open symbols for 8 = 35° and 0.35 < M < 2.10, 

and the solid symbols for 8 = 18° and 1.05 < M < 2.37. 

Figure 94 indicates that for M < Mq PT ^ r , injection at x^/L s = 
0.150 for 8 = 18° and 8 = 35° produces essentially the same cooling 
effect (i.e. reduction in STp C /ST Q ). It appears that injection at 8 = 
35° Is slightly more effective than 8 = 18°, however the difference is 
at most 4% and within the limits of the least-squares representations 
of the data. As M increases past for both configurations 

(i.e. 8 = 18° and 3 = 35°), the STp^/STg values are represented by the 




Stanton Number Reduction, ST„„/ST 
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circular symbols. Since M Q p T for 3 = 18° is greater than that for 
35° (1.05 compared to 0.85), a small range of M will exist in which the 
results for 18° are still represented by the least-squares curve while 
the values of STpg/STg for 35° are shown as open symbols. For 0.85 < M 
< 1.05, the 18° configuration (still shown by the solid line) is more 
effective than 8 = 35° (now shown as open symbols). For 1.05 < M < 2.37, 
the data for both injection schemes are illustrated by symbols. The near 
injection region ((x-x^J/dg * 11), characterized by the shaded region, 
shows a steep rise in STp C /STg with increasing M (illustrated by the 
direction of the large vertical arrow). At the high blowing ratios for 
large distances downstream (large horizontal arrow) the data is quite 
scattered, generally showing reduced cooling (STpg/STg) increasing). 

Figures 93 and 94 indicate that increasing the injection angle from 
8 » 18° to 3 ® 35° results in the following: 

(1) a decrease in the value of Mgpy with increasing 

(Mgpx ^ar decreases from 1.05 to 0.85 for 8 inc casing from 
18° to 35°, at x./L s « 0.150) and 

(2) a decrease in the cooling effectiveness (i.e. larger STpg/STg) 
as injection moves closer to the stagnation point. 

Item (1) above can be explained at x^/L s » 0.150 by analyzing 
Mgp.p for both the 18° and 35° in relation to the normal component 
of the coolant momentum flux. It is postulated that the optimum value 
of M is associated with coolant penetration into the boundary. There 
is a characteristic penetration distance (y*) corresponding to the value 
of Mgpj is assumed now that y* is proportional to the normal 

component of coolant momentum flux, expressed as 



2 

Referencing Eq. (81) to the constant freestream momentum flux, p m V„, 
the following is obtained 

* p r Vp 

_L = k * Jd- s in 3 « k*i sin (3 (82) 

pX pX 

where k* - proportionality constant. Equation (82) can be expressed in 
terns of M as 

-d = k* (Jp) M Z sin 6 (83) 

P X Pc 

It is assumed that should correspond to a similar penetration 

distance for each injection angle. In the present experiments, since 

o 

(p m /p c 5 as well as pj/^ remained constant at a particular x^/L s , Eq. 

(83) indicates that Mq P j and {3 should be correlated as follows, 

m opt,near 1 sin B 1 * h opt,near 2 sin h ~ conStant 

Using the results at x^/L s = 0,150 and Eq. (84) one obtains, 

0 = 18°, ^0PT,NEAR 13 M sin 0 “ 0.341 

3 = 35° s ^ 0 PT,NEAR = 0,85j ^ sin 8 * 0,414 

The agreement with Eq. (84) suggests that the normal component of the 
coolant momentum flux may be useful to correlate the coolant blowing 
conditions for optimum film cooling performance. 
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Using Eq, (84) for the results obtained at x./Lg = 0,0833, 3 = 18°, 
Mqpt - 1.3, gives M sin 3 ® 0.522, indicating that the value of 
the constant varies with x^/Lg. Although the data for 3 = 35° did not 
exhibit an optimum H for values up to M « 1.8, there was no significant 
change in STp^/STg for M greater than J.Q. This would correspond to 
M sin 3 = 0.574 which is in reasonable agreement with the data for 18°. 
It is apparent that the film cooling performance varies with the loca- 
tion of injection, x^/Lgj or more importantly, the local freestream 
acceleration at injection. The following section, illustrates the inter- 
dependence of the injection angle and the local freestream acceleration 
as reflected in the film cooling performance for different injection 
locations. 

V.3.2.4. Effect of Local Freestream Acceleration at Injection 

For the values of injection angles investigated in this stud/ (18° 
and 35°), the effect of incorporating shallower injection angles on the 
film cooling effectiveness (i.e, 35° to 18°) does not appear to be 
significant in itself. However, under freestream conditions that exist 
for injection near the stagnation region, the cooling effect of the 
injected film shows an interrelation between the local value of the 
freestream acceleration (characterized by K-) and the angle of injection, 
3 . 

Figure 95 shows the variation of STp C /ST Q with X for B « 18° at two 
injection locations, x../Lg = 0.0833 and x./Lg “ 0.150, The least-squares 
curves from Figs. 88 and 89 representing the data correlations are shown: 
the broken curve is indicative of injection at x^/Lg = 0.150, while the 
solid curve represents injection at x^/Lg = 0.0833. By moving the 







injection site nearer to the stagnation point, two pertinent freestream 
parameters are changed: (a) the local' freestream acceleration parameter, 

increases, and (b) the ratio of the momentum boundary layer thickness 
to the hole diameter ratio, O-j/dg, decreases. At - 0.150, the 

acceleration prediction scheme) is 0.72 x 10“ 6 with G^/d Q =0.013. At 

x^/Lg - 0.083, the values are = 7.9 x 10~ 6 and G-j/dg = 0.013. 

Figure 95 illustrates that up to X « 30 film coolant injection 

nearer stagnation (i.e. x./L*. = 0.0833) is less effective in reducing 
ST fc /ST q than for injection at x^Lg = 0.150. As X is increased from 
7 to 30, the difference in STpg/STg between the two injection locations 
narrows from 22% at X = 7 to zero at X = 30. The difference in cooling 
effectiveness may be attributed to the effects of local acceleration. 

More than an order of magnitude increase in K.. (0.72 x 10 to 7.9 x 10 ) 
is represented in Fig. 95. Goldstein and Jabbari (64) report a 10% de- 
crease in the cooling effectiveness from a double row of streamwise 
angled holes (3-diameter spacing) for M between 0.5 and 1.0 for changes 
in K.j from 0 to 1.0 x io" 6 . Kays, Moffat and Thielbahr (39) measured 
a 20 to 25% increase in the local transpiration cooled Stanton number 
(at a strong transpiration blowing ratio of 0.0062) for an increase of 
K from 0 to 0.77 x 10“ , The range of K. exhibited for the subject 
investigation extends to values nearly an order of magnitude greater 
than the cited references. Consequently, significant effects of K. on 
STpg/STg near stagnation would appear to be consistent with published 
data. 

In addition to the change in 1C. for the data shown in Fig. 95, 

9^/dg also changed significantly, from 0.036 to 0.013 for x./L s = 0.150 


to x^/L^ ~ 0,0833. The effect of a decreasing boundary layer thickness 
should result in an improved cooling effect, due to a more shallow tra- 
jectory of the coolant jet. In reference to Fig. 95, the influence of a 
decreasing value of should result in the broken curve both shift- 
ing slightly downward and extending to smaller values of A. However, 
decreasing Q-jAIq accompanied with increasing K. resulted in just the 
opposite change in STpq/STq. It appears then, that the effect of in- 
creasing K.j may be responsible for the differences in ST^/STq at the 
two injection locations. In view of the opposing trends of K* and 0-j/dQ* 
the absolute level of the effect of t(. may be even greater than illus- 
trated. 

As the injection angle gets steeper, the effects of become more 
prominent. Figure 96 presents STp^/STg as a function of A for 3 - 35° 
at the same injection locations analyzed in Fig. 95. The solid line is 
the least square.- representation of the data for x./L s = 0.0833, while 
the broken line is the least squares representation of the data for 
x./Lg = 0.150. 

At A = 7, the value of STpq/STq is 48% higher for x../Lg - 0.0833 
than for the configuration farther from stagnation {x./Lg = 0.150). 
Between A = 7 and 50 the difference in STpQ/ST 0 between the two configu- 
ration narrows. For A > 50 a constant difference of approximately 4% is 
observed. The effect of increasing K. appears to be more significant 
for the steeper injection angle (3 = 35°) than previously observed in 
Fig. 95. The opposing trends of K. and e^/d^ can be argued similarly as 
for Fig. 95, indicating that the acceleration influence is even more 
pronounced than indicated in Fig. 96. 
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In comparing Figs* 95 and 96, it appears that the effects of accel- 
eration for injection at 3 = 35°, corresponding to a rnore penetrating 
jet trajectory, shows a greater dependence of STp C /ST Q on the local 
value of the freestream acceleration at injection than injection at 
8 = 18°, Specifying shallow injection angles would appear to be most 
important in regions of high freestream acceleration. 



VI. EXTENSION OF RESULTS TO DESIGN APPLICATIONS 


The data for Cases I - V of the subject investigation have been 
presented indicating STp C /ST Q results for a film cooling scheme which 
models a high temperature, high pressure turbine inlet environment. The 
results of this study were obtained under model operating conditions of 
T y = 689°K (1240°R) , P T - 3.1 x 10 5 N/m 2 (44.67 psia), - 294°K 
(530°R) and a specified vane geometry and orientation. The values re- 
ported for STp C /STg are characteristic of a dimensionless coolant tem- 
perature, 0£ - 1.03 or cool ant- to-gas density ratio of P c /p ro * 2.15. 
Extrapolation of the data to approximate STp^/STg for other values of 
e£ can be accomplished by utilizing the superposition principle for a 
linear energy equation which is discussed below. 

The superposition principle resulting from the energy equation 
applies when property variations (more specifically, density) within the 
thermal boundary layer may be assumed negligible. Two approaches to film 
cooling predictions, the adiabatic wall temperature method and the heat 
flux reduction method, rely on the linearity of the energy equation to 
extend experimental model results to high temperature and high pressure 
engine conditions. 

The following sections will present a discussion concerning the 
application of film cooling data to high temperature and high pressure 
conditions. A parallel development of the two models and qp^/qg) 

for predicting the local, reduced heat flux is included. A discussion 
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is also included illustrating the use of each technique. Finally, the 
limitations of each method in making film cooling predictions can be 
noted in light of recent information concerning the effects of large 
density (pq/pJ variations. 

VI. 1. Film Cooling Heat Flux Reduction Predictions 
A comparison will now be made of the adiabatic wall temperature 
model and the heat flux reduction model for predicting the heat flux 
under film cooling conditions, 

VI. 1.1. Adiabatic Wall Temperature Method 
Using the concept suggested by Eckert (5), the heat flux to a film 
cooled surface can be calculated if both the hydrodynamic and thermo- 
dynamic effects of injection can be determined. The thermodynamic effect 
is manifest in a reduced boundary layer gas temperature near the wall 
due to cool mass addition. For an adiabatic surface, this reduced 
boundary layer temperature is referred to as the adiabatic wall temper- 
ature which has frequently been measured in film cooling experiments. 

The hydrodynamic effect results from enhanced mixing due to injection, 
resulting in a more effective transport mechanism for mass, momentum 
and energy. This effect is measured by noting the change in local heat 
flux due to injection of a coolant at a temperature equal to that of the 
freestream. The combined thermodynamic and hydrodynamic effects on the 
local heat flux with film cooling can be expressed as 


q FC = h 0£=O (T ADW “ T W,Fc) 


hg i^q = h l as discussed in Chapter 1 1. 1.1.1. 
C 


( 85 ) 
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where h fl ,„ n = local heat transfer coefficient reflecting the hydrodynamic 
u tf u effect (T c => T t ) 


ADW 


« adiabatic wall temperature reflecting the thermodynamic 
effect {Tq < Ty ) 


The fractional reduction in the local heat flux due to film cooling can 
be expressed as 


q Fc _ V tT /m " t w,fc* 
q o ' h o ( t t„ - T w,o> 

where qjj = heat flux without film cooling 

hg = local heat transfer coefficient without film cooling 

Recalling the definition of the adiabatic wall effectiveness. 


(86) 


n, 


t t w “ t adw 

\ ' T c 


'ADW 

one can express in the following manner 

T ADM = T T " n ADW “ V 
such that Eqs, (86) and (88) combine to give 


(87) 


( 88 ) 


q FC _ ," 6 C=° 

Ti V h ) 


,=n < T t ■ T W,FC^ 


< t t -V 

CO w 


q 0 h 0 < T T " T W,0^ nADW ( T T " T W,(P 

CO CQ 


(89) 


If a constant wall temperature is prescribed (T^ ^ Q ) to evaluate 

the fractional reduction of the heat flux due to film cooling, Eq. (29) 


gives 


qr r h e;=0 

#=(-*£-) l-n 

q o 


h, 


'ADW S C 


(90) 
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Equation (90) states that the fractional heat flux reduction is a linear 
function of 0£. Once a film cooling geometry is fixed and a coolant 
blowing rate is specified, (h 0l= Q/h Q ) and can be determined from 
experimental data* To achieve a certain fractional heat flux reduction 
under these conditions, 0£ is then uniquely determined. 


VI. 1.2. Direct Heat Flux Measurement Method 


Instead of conducting experiments to determine n ADW and (h Q , =0 /h 0 ), 

C 

the fractional heat flux reduction can be measured directly. This method 
compares the measured heat flux at a particular location for a specified 
set of coolant conditions (i.e. 9£, M, B, S-j/dg, etc.) to the heat flux 
measured without film cooling. The ratio is expressed 


*FC _ "FC 


< t t “ t w,fc ) 


% 


TU 


00 


T W,0 J 


(91) 


where qjig/qg is a function of all the coolant injection parameters. 

Note, in contrast to Eq. (85), hp C is defined in terms of the gas-to-wall 
temperature difference, (T t - T,, rf >). 

To extend the results from experiments wherein qp C /qjj is measured 
to coolant temperatures other than those investigated, the linearity of 
the energy equation allows the following development (4). At a fixed 
position and coolant flow condition (i.e. M) the boundary layer temper- 
ature is a function of the boundary condition, 0£ and can be expressed 


T ( 0 c> = + e£ CT(9^ (2 ) - T(e£ (1 )] 


(92) 
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where e£ ^ = dimensionless coolant temperature for T c = T c -j 

9 C 2 = d ™ iensionless coolant temperature for = T c 2 

T(0p) = solution for T at arbitrary 0i, and specified x* y, 
u and M. 0 


Recalling 



HALL 


(93) 


and using Eq. (93) to operate on Eq. (92), the following results 


qji c (ei) = qp C (0Q 9l ) + 0Q “ q FC^°C»l^ 

Referring Eq. (94) to the condition with no film cooling, q". 


q FC^ 9 C^ q FC^ 9 C,1^ 
q 0 " q 0 



(95) 


Equation (95) is a linear relation in e£, with [<lpc( 0 c -j)/9q 3 as the 
intercept of the ordinate, and [q£ c (0£ -j )/9 q - qp C (0Q 2^ q 0^ as tfie 
slope. If qj c (e£ -j = Oj/qJ represents the hydrodynamic effects of 
injection (e^ -j = 0 with = Tj ), the slope can be evaluated by letting 

9 C,2 corres P° nd to the case in which or qp C (e£ A0W ) = 0. For 

the adiabatic case 



e C,ADW [sl °P^ 


Therefore, Eq. (96) gives 


slope « 


fo( 6 c -0 > 

q 0 ■ °C,ADW 


(96) 


(97) 
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Nov/ Eq. (95) can be written as 

<fo( 9 r.) qFC< e c =0 > °c iFct°> 

q 0 q 0 0 C,ADW q 0 


(98) 


Since T T " T C i 

Q I ^ 00 = * 

C,ADW T t - t adw n ADW 

CO 

therefore, Eq. (98) becomes 



“ n A0W °C^ 


(99) 


(100) 


If a constant wall temperature is prescribed (T W#FC = T w# q) to evaluate 
the effect of film cooling in reducing the heat flux, Eq. (91) shows 
(for e£"0) 


qFc (e c° 0) V 0 

q 0 h 0 


( 101 ) 


Therefore, Eq. (100) can be reduced to 



hp I = n 

~hj] [1 ' n ADW 



( 102 ) 


Equation (102) for the heat flux reduction model is now identical 
to Eq. (90) which was developed using the adiabatic wall temperature 
model . 

Since Eq. (102) is in the linear form 




( 103 ) 
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measurements of Rf C ( 9 c)A]q at onl ^ tw0 values °f 9 £ can determine A and 
B. There will be a unique set of A and B for each blowing ratio and 
dimensionless distance from injection {geometry held fixed). In general 
A and B will be different for each change in a parameter associated with 
film cooling. For similar cooling configurations and injection condi- 
tions, A and B are specified from experimental results. Consequently, 
the fractional heat flux reduction can be determined for any choice of 
the coolant temperature (9^). 

The adiabatic wall temperature and heat flux reduction models were 
shown to yield identical results in predicting the effects on the local 
heat transfer due to film coolant injection. It would be helpful, at 
this point, to illustrate the basic procedure one would follov;, using 
one of these aforementioned methods, to predict q^/qQ for high temper- 
ature and high pressure conditions. 

VI. 2, Application of Film Cooling Heat Flux Prediction Methods 

As an example of applying one of the aforementioned film cooling 
prediction schemes, consider the problem of cooling a critical area in 
the leading edge region of an inlet turbine vane. Assume that, prior to 
incorporating film cooling, the following conditions are specified: 

Ty “ freestream gas total temperature 

co 

T w 0 = wall temperature without film cooling 

Re ra - freestream Reynolds number 

x/L<j = location of the critical area relative to stagnation 

T. = temperature of the coolant inside the turbine vane 
used to maintain T^ Q . 
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It is desired to introduce film cooling to reduce the vane wall temper- 
ature in the critical area. To simplify the example, assume that the 
following film coolant flow and injection geometry parameters are speci- 
fied: 

M = coolant blowing ratio 

x./U - location of injection holes relative to stagnation 
1 ^ 

3 - injection hole angle 

d Q = injection hole diameter 

(x-x.)/dg = location of critical area relative to injection holes 
S-j/dg - hole spacing ratio. 

We ask, what value of Tg is required for film cooling to reduce the wall 
temperature by some fraction, n? 


T W,FC = n T W,0 


(104) 


Or, in terms of the dimensionless coolant temperature. 


e 




what value of 0£ is required to reduce the wall temperature by some 
fraction n? 

It is assumed that measurements of qjig/qjj were taken at reduced 
(simulated) conditions for an identical film coolant injection geometry 
and consequently, a unique relation has been determined from the simu- 
lated conditions between qpg/qjj and 0g for the specified M, (x-x^ )/dg, 
etc. for use in the prediction. 
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Prior to introducing film cooling,, the heat flux can be expressed 


as 


<0 ■ h o < t t - T w,o> " 

03 


(T t - T.) 

G3 

TiT + R W + R 1 


(105) 


where 


Ry = thermal resistance of the vane wall 

R. = thermal resistance due to convection on the vane 


inside surface. 

In a similar fashion* the heat flux with film cooling can be written as 


q FC “ h FC ( T T " n T li,0^ = , h 

03 I / 


< t t _- T i> 


*Src + Rw + Ri 


(106) 


where and are assumed to be the same as in Eq. (105). The heat 
flux ratio can be formed from Eqs. (105) and (106) to give 


q FC , h FC . _ 
"0 ' h 0 


(ER) 


0 


U-(J-) t C2R) - f 

h 0 n FC u n 0 


(107) 


where 


1 - n . 


W,0 


b = 


1 - 


'W,0 


(107a) 


^0 * h 0 * R VJ * R i 


(107b) 


Equation (107) can now be solved for hp C /hg to give the following 


h FC _ b " a 
h« “ 1 - a 


(103) 



where 


(108a) 


O/h 0 ) 

a “ 

For the subject example, b can be determined from Eq. (107a) for a 
selected value of n, the desired fractional reduction in wall temperature. 
The value of "a" can be determined from the information specified at the 
outset (Tj , T iS T WjQS Re w , x/L s ) and Eq. (108a). Thus, Eq. (108) can 
be used to evaluate h p g/hg for the specified conditions and the selected 
value of n. Consequently, Eq. (107) can be employed to determine the 
value of q£ c /qjj for the selected value of n. 

Recall that q pc /qg as a function of 0g was assumed known from the 
model experiments (under simulated conditions) for the specified film 
coolant injection conditions (M, (x-x^/dg, etc.). These results would 
be in the linear form 


q FC „ » „ fll 

qJT ' A ‘ 


(109) 


With the value of q pc /qg known from Eq. (107), the value of G£ required 
to reduce the wall temperature by some fraction n can be determined 
from Eq, (109). 

The assumption carried along during this example was that data for 
qjlg/Qg as a function of 0^, obtained from model experiments are applicable 
to problems involving large temperature differences between the coolant 
and mainstream. For example, at a value of = 1 the predicted heat 
flux ratio, qjig/q*, f° r a P rot °typs condition of Tj = 1922°K (3460 G R), 

T w FC = 961 °K (1730°R) and Tg = 961 °K (1730 G R) is presumed to be identi- 
cal to a model condition of T ? = 330°K (594°R), T w>pc = 290°K (526°R) 
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and T c - 290°K (526°R) providing identical film coolant conditions (i.e. 

M, 6, S ]/ d o* Cx-x-J/d^, etc.) are employed. For both cases Q£ = 1; how- 
ever, the respective values of density ratio (pq/p^) are 2.0 and 1.14. 

For a constant value of M, the velocity ratio (Vq/V^) will vary inversely 
with the change in the density ratio. This variation in velocity ratio 
at constant M and can cause difficulties in extending data, obtained 
for a low density ratio (-1) to predictions of ^pq/Qq ^ or high temperature 
applications where a density ratio of the order of 2 is more typical. 

Recent evidence presented by Pedersen (24) has shown that the 
laterally-averaged adiabatic film cooling effectiveness, (determined 
using a mass transfer analogy ), depends upon the coolant-to-freestream 
velocity ratio, Vq/V^, for a fixed geometry and blowing ratio. Figures 
97(a), (b) and (c) from Pedersen's results show as a function of 
x'/dg (where x'=x-x..) for blowing ratios of 0.213, 0.515 and 1.05. The 
cool ant- to-gas density ratio was varied between 0.75 and 4.17. The 
magnitude of the effect of P c /p w (as manifested by V c /V 03 ) increases with 
increasing M. Also, the dependence of n™ on the density ratio is most 
pronounced near injection ( K '/ d Q < 30). 

Figure 97(a), at M = 0.213, represents the variation in rTjy with 
x'/d 0 for the velocity ratio in the range 0.051 V c /V ra <0.282. As 
discussed in Chapter V.3.2.I., a value of velocity ratio near V c /V ro « 0.5 
appears to represent an equilibrium condition below which film cooling 


The adiabatic cooling effectiveness, n/\py, and impervious wall 
effectiveness, my, are equivalent when the turbulent Prandtl and Schmidt 
numbers are identical. Pedersen's results included a verification of 
this point (24). 


Figure 97(a). Effect of Density Ratio on Laterally Averaged 
Cooling Effectiveness at M = 0.213 (24). 
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performance increases with M (or V^V^JCdue to increased thermal capacity 
of the coolant) and above which the performance decreases with increasing 
M (due to enhanced mixing). Consequently, for the range of density ratio 
associated with Fig. 97(a), Vq/V^ never exceeds 0.23 and the mixing 
effects should be minimal. This is substantiated by the small shaded 
region which characterizes the variation hi n IW due uO the variation of 
density ratio from 0,76 to 4.17. 

For M « 0.515, shown in Fig. 97(b), the velocity ratio ranges from 
0.12 to 0,68. The variation in for P^/P ro between 1.50 and 4.17 is 
small since this corresponds to V^/V^ less than 0.34. However, for 
P c /p ro between 1.18 and 0.753, Vq/V^, varies from 0.44 to 0.68. Since 
this velocity ratio range includes and surpasses the previously noted 
value (~0.5) which may be cham£tecis±ic of significant mixing, this may 
explain the decrease in with decreasing P c /p ro shown by the data. 

Figure 97(c), for M = 1.05, shows an even wider separation between 
the measured effectiveness values for each value of density ratio. The 

higher values of density ratio (4.17 and 1.97) correspond to velocity 

ratio values of 0.25 and 0.53, respectively. The lower values of density 
ratio are characteristic of a coolant-to-gas velocity ratio well above 
0.5. Consequently, the data exhibit a strong dependence of on Pq/p„. 
However, as Pq/p^ increases above 1.50 (Vq/V^ > 0.7), the low value of 
n 9 presumably mixing dominated, is relatively insensitive to further 
increases in V c /V to . One can also expect that, for higher values of M, 
all curves begin to collapse into a single correlation since V^/V^ would 
be greater than 0.5 for all values of density ratio investigated. 
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Pedersen's findings indicate that the prediction methods illustrated 
by the previous example may prove inadequate for determining qpg/qjj under 
conditions wherein there are large differences in density ratio from the 
model data to the predicted case. The main difficulty, of course, arises 
from the necessary change in Vg/V^ to maintain the same M for increases 
in pg/p m . If no velocity-related mixing effects occurred, the extension 
of the aforementioned prediction methods to other than small density 
ratio variations would prove valid. 

To illustrate the magnitude of the velocity ratio effect. Fig. 98 
was prepared using Pedersen's results for effectiveness values as a 
function of pg/p^. The predicted fractional reduction of heat flux due 
to film cooling is shown as a function of 9g. The open symbols represent 
results for M = 0.515 at x'/dg = 20 (Fig. 97(b). The solid symbols are 
associated with M = 1.05 at the same distance from injection (Fig. 97(c)). 
These too values of M and x‘/dg = 20 were chosen since the effectiveness 
was previously shown to be sensitive to pg/p OT variations at these con- 
ditions. Circular symbols represent values of qjig/qg predicted using 
the adiabatic wall effectiveness method assuming that was not a 
function of Pg/p^. Values for n^.j were obtained from Pedersen's data 
by interpolating for p c /p ot = 1-0. The square symbols show the heat flux 
reduction predicted using the results for rj^gy from Pedersen at the 
correct value of pg/p ra . 

A gas temperature of 1922°K (3460°R) and a constant wall temperature 
of 961°K (1730°R) were assumed for the sample calculations of Fig. 98. 

The value of Tg was varied between Ty and 0.24 Ty (Gg from 0 to 1.52) 

‘■■CO co 

and n^gy was evaluated at the appropriate value of Pq/p^ in one case and 


\ r 


0.7 



Figure 98 . 


Predicted Film Cooling Stanton Number Ration Using Both 
at p-/p. * 1 and ^_ DW as a Function of P C /P OT - 


r\. 

i 
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assumed independent of Pq/p w in the other. The following relation was 
used to determine the heat flux reduction due to film cooling, 



* ( 


hrt i „n 

-hf-) <’ - "ADW 



(no) 


The hydrodynamic effects of blowing on the heat transfer coefficient 

* 

were neglected to simplify the comparison (i.e. h 0 , =o /h o = 1). 

The reader is now referenced to the open symbols illustrated on 
Fig. 98. They correspond to a blowing ratio of 0.515 at x'/dg = 20. 
Focusing on Q£ = 1.0 (a representative condition for turbine cooling 
applications), a 3.62 difference in qp C /q{J is shown between the pre- 
diction which assumes is not a function of Pq/p^ {circular symbol) 
and the prediction which accounts for the variation with Pq/p^ 
{square symbol). As 0£ increases from unity, the percentage difference 
between the two prediction methods increases. 

Examining similar predictions for M = 1.05 at x'/dg = 20, repre- 
sented by solid symbols in Fig. 98, a difference of near 1955 is shown 
between the two prediction techniques due to the effects of p~/p . As 

CO 

for the previous case (M = 0.515), the difference in the predicted 
reduction of heat flux increases significantly as 0 0 increases above 1.0. 
As 0^ is increased, P c / pancreases (maintaining Ty and Ty constant) 


It is recognized that at the higher value of M_hydrodynamic blow- 
ing effects may k e significant. However for nADW = 0ADW {pc/p<J> the 
effect of blowing will diminish as BA increases from zero since Vf/V^ 
must decrease with pq/p^ increasing jfi.e. square symbols on Fig. 98)? 
Therefore the percentage differences between the two prediction methods 
may, in fact, be even greater than shown. 



with a subsequent decrease in Vq/V^ required to maintain M constant. 
Consequently, injection related mixing is diminished resulting in an 
ever improving cooling capacity relative to the prediction which ignores 
the dependence of on Pq/p^. 

The difference between the two predictions of the reduction of heat 
flux due to film cooling is particularly important to note in relation 
to calculating qp C /qjj for turbine cooling applications. Typical values 
of 0£ for turbine conditions would be somewhat greater than unity; the 
coolant temperature would approach the wall temperature only after 
extensive use as an internal convective coolant. Therefore, extending 
predictions of the reduction of heat flux to lower coolant temperatures 
{0(1 > 1.0) would result in progressively more conservative estimates if 
the effects of V^/V^ (as manifest by Pq/p^) are not considered. 

The following section presents a discussion, in light of Pedersen's 
findings, concerning the extension of the results of the subject inves- 
tigation to other than the modeled turbine conditions. 

VI. 3. Use of Present Data 

The data for the reduction of heat flux from the subject investi- 
gation are a true representation of the film cooling capability of an 
injection configuration consisting of a single row of spanwise angled 
holes in the leading edge region of a turbine vane. Specifically, the 
measured values of qj^/q^ are characteristic of an average gas-to- 
coolant temperature ratio (T t /T r ) of 2,18 with a gas-to-wall temper- 

co U 

ature ratio (T t /T,,) near 2.15. Consequently, all of the data from 

oo 

this investigation correspond to 0 q - 1.03. Typically, turbine cooling 
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conditions are represented within a rather narrow range of 9^. A limit- 
ing case in which the coolant air has been utilized to the maximum for 
internal convection prior to ejection as film coolant occurs at G£ = 1 
(T c = T w ). However, values of 9^ up to 1.2 or 1.3 are conceivable. 
Extrapolation of the results from the subject investigation to 0 q values 
in the range of 1.3 would be useful to permit a wider application of the 
data. 

Consider extending the results from Case I herein to predict the 
reduction of heat flux due to injection of film coolant with a lower Tg, 
i.e. 6 q > 1.03. For example, let us examine film cooling of the loca- 
tion (x-x^J/dg = 10.45 using a blowing ratio of M = 1.05. For the sub- 
ject experiments, 9£ = 1.03, and the corresponding value of qp C /qQ is 
0.71 (based on the least squares curve from Fig. 44). Figure 99 illus- 
trates the extrapolation of qj^/qg to higher values of indicating 
the uncertainty in the slope of the extrapolation curve. Depending upon 
the value of M, both che slope (for a linear approximation) or the 
degree of the polynomial (for the case in which Pq/p ot effects can be 
determined) representing the variation of qp C /qJJ are unknown. The most 
significant item to note in Fig. 99 is that the single point shown is a 
true representation of qj^/qg for a density ratio, p c /p m , of 2.18. The 
major difficulties in extrapolating qj^/qjj are: 

(1) the choice of a representative value of q^/q^ at = 0 and 

(2) the determination of the proper functional representation of 

qpc /q o vs e c' 

Fortunately, the value of qp C /qg at 9^ = 0, i.e. the hydrodynamic 
effect, is well documented for a range of blowing ratios and injection 
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Figure 99. Possibilities in Predicted Film Cooling Haat Flux 
Using A Single Data Point. 




260 


geometries (20). Depending upon M and (X”X^)/dg # the value of q FC /q 0 
for Sj = 0 (i.e. T« = T T ) may vary from 0.95 to 1.3 for a typical single 

^ ^ <X> 

row configuration of S^/d Q =3,0. In the absence of density ratio 
effects, a linear extrapolation of qjig/qg from the subject data point 
to other values of 6g can be employed. Figure 99 gives an indication 
of the extremes of the linear extrapolation depending on the magnitude 
of the hydrodynamic effect. 

The largest error between the linear extrapolation and the true 
variation of qp C /qJj with 9£ will occur when the true representation is 
most nonlinear. Recalling the discussion from Fig. 98, the predictions 
of qp C /qg using Pedersen's data showed the greatest nonlinearity with 
0g at the highest values of blowing ratio (M = 1.05). In addition, the 
regions rearest injection (i.e. x'/dg « 10) were more sensitive to 
changes in P c /P ro or Vg/V^ than farther downstream. Pedersen's results 
have been used to provide a comparison between the proposed linear 
approximation and the true, nonlinear variation of qjig/qg vs. e£ as 
presented below. 

For Case I of the subject study, with (x~x.j)/dg = 10.45, e£ = 1.03, 
Tf /T w = 2.15, q“ c /qJJ has a value of 0.71 for M = 1.05. This corresponds 
closely to a case for which Pedersen presents data (i.e. x'/dg = 10, 

M = 1.05) showing the influence of Pq/p m . Figure 100 shows a linear 
extrapolation of q^/qg vs. Qg through the subject data point with the 
hydrodynamic effect (at 6g 3 0) evaluated from data presented by Eriksen 
(20). Also presented in Fig. 100 is a nonlinear representation of 
q Fc/ q 0 vs * ^ eri * vec * us ^ n 9 the data from Fig. 97. Reference to Eq. 

(100) shows that q^/qQ can be expressed as 
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Equation (111) was used to generate the nonlinear variation of q Fc/ q 0 
vs. 6 q as follows. For a selected value of ©£ (maintaining T-j- /T^ - 
2.15), a value of T c /Ty and p c /p w was calculated. For that density 

CO 

ratio, (i.e. rf IW ) was determined from Fig. 97(c) at x'/dg = 10. 
Noting that M is to be held constant at 1.05, the value of Vq/V^ was 
computed. This value of the velocity ratio dictates the hydrodynamic 
effect manifest in q^G^ = 0)/q”. A value of q jl^/q" for T c - Tj 
(i.e. G’ = 0) was obtained at the corresponding value of Vq/V^ from 
Eriksen (20) for a cooling configuration identical to Pedersen's. Eq. 
(Ill) was used to calculate qp C /q|j for selected values of 0^ with the 
nonlinear variation shown in Fig. 100. 

In Fig. 100, the region marked "turbine cooling regime" provides 
an estimate of the error in the predicted value of q^/dp due to the 
linear extrapolation for the subject case. For a known value of qj^/qQ 
at e,i = 1,03, the heat flux ratio at 0£ = 1.0 and e£ = 1.3 predicted 
from the linear relation would be in error by only 1 .455 and 10%, respec- 
tively. 

Although it is recognized that data for are not ava ^“ 

able for exactly the same cooling configuration as used in the present 
study, values obtainable for similar geometries would probably be rep- 
resentative. For a variation of q^/q^ with 8£ that is truly linear, 
small errors in the assumed value of )/ c *o will not result in 

significant errors for extrapolation in the range e£ = 1.0 to 1.3. If 
q Fc/ q 0 a i* 0 c = 0 1S in evTOr by as muc ^ as 10"» the error incurred by 


263 


extrapolating linearly from - 1.03 to e£ = 1.3 is less than 4%. Figure 
98 illustrates that the linear relationship improves for low values of M 
and large values of x'/dg and was a good approximation for that case for 
M - 0.5 and x'/d g as low as 20. 

In summary, it has been shown that the data of the subject investi- 
gation can be extended to predict the reduction in heat flux due to film 
cooling for a narrow range of coolant temperatures. A linear extrapola- 
tion of measured values of q^/qg at 9£ - 1.03 should be accurate within 
10% for 1.0 9g 1 .3. The methodology ( ,used to extend the current 
results to a wider range of turbine cooling conditions can be summarized 
as follows: 

(1) For a particular configuration {i.e. Cases i - V), a value of 
q F C / q 0 ( 1,e * STpg/ST Q ) 15 k nown f° r each M and (x-x.J/dg. 

(2) Choosing values of (x~x.)/d Q and M of interest, the appropriate 
qpg/qJJ is noted. A coordinate pair [(qjig/q^, 9g = 1.03] can 
be plotted on a (qp C /qjj) “ plane. 

(3) The value of qj^/qg at e£ = 0 for the M of interest is deter- 
mined from a published source for a similar geometry, and 

(4) A straight line between the two heat flux ratios at 0g = 0 
and 6g - 1.03 is then used to interpolate for dp C /q|J at other 
values of 6g. 
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VII. COMPARISON OF RESULTS WITH PUBLISHED DATA 

As noted in Section II. 2. (Literature Review of Important Cooling 
Parameters), there is a vast amount of information pertaining to previous 
analytical and experimental studies of film cooling. To provide some 
perspective it is desired to compare the results of the present study 
with published data conducted under similar cooling conditions. 

The author knows of no available experimental data which would 
permit a direct comparison. Two features of the subject investigation 
set the subject results quite apart from published data: (1) film cool- 

ing heat transfer measurements were made downstream from a single row of 
laterally angled holes at a surface spacing- to-diameter ratio of 4.0, 
and (2) the heat transfer measurements were taken with a cool ant- to -gas 
density ratio of 2. IP. The work reported by Pedersen (24) provides 
similar data for the adiabatic wall effectiveness with density ratio in 
the same range. Though the injection geometry is different (streamwise 
angled holes), his investigation is perhaps the best reference for com- 
parison to the present data. By correcting Pedersen's cooling effective- 
ness data to correspond to the spacing ratio used in the subject inves- 
tigation (Pedersen had S^/d^ - 3.0) a comparison can be made between the 
heat flux reduction due to film cooling from a row of holes (Pedersen) 
angled at 35° in the streamwise direction and that for injection from a 
row of holes angled at. 35° in the spanwise direction (present study). 
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Pedersen conducted a mass transfer study. His results, however, 
are analogous to heat transfer. The data are reported in the form of a 
laterally-averaged, adiabatic effectiveness, n"^. 

Equation (90), developed for the adiabatic wall temperature method 
was used to convert Pedersen's results to a heat flux reduction. Recall 
that Eq. (90) is of the form 


q" h e'=0 

q o = * h o ^ nADU 



( 112 ) 


Three items are of importance in using Eq. (11 ;> for comparison to the 
present data: 

(1) values must be adjusted to correspond to a spacing-to- 
hole diameter ratio similar to the present study, 

(2) values of n^.j must be consistent with the dimensionless 
coolant temperature (i.e. effectiveness data must be at same 
density ratio as subject study) and 

(3) the hydrodynamic "blowing" effect, (h Ql=0 /h 0 ), must be deter- 

C 

mined as a function of M and distance from injection. 

Item (1) through (3) will now be discussed in more detail in the follow- 
ing paragraphs. 

To contrast the 35° streamwise injection configuration of Pedersen 
to the 35° spanwise injection of the present study, compensation must be 
made for the additional geometric parameter not common to both studies; 
namely, the spacing-to-hole diameter ratio. Pedersen had an S-j/d^ "3.0 
in contrast to S-j/d^ = 4.0 for the subject investigation. Liess and 
Carnel (16) {as noted in the literature review) examined the spacing 
ratio effect on the cooling effectiveness for streamvrise injection at 


D .. 
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a = 35°. Spacing ratios of S^/dg - 2.22, 3.33 and 4.0 are reported for 
a single value of the blowing ratio of 0.85. From their results, 
laterally-averaged adiabatic film effectiveness values can be determined 
for each spacing ratio out to x'/dg = 25. 

To correct Pedersen's values of rj^py to a spacing-to-hole diameter 
ratio of 4.0, the following assumptions were made using the data of 
Li ess and Carnel (16): 

(1) the percentage change in due to a variation in S-j /d Q is 
independent of the value of M, and 

(2) a smooth extrapolation of from x'/dg = 25 to xVdg = 80 
for the three values of spacing ratio investigated by Li ess 
and Carnel. 

Interpolating between S^/dg - 2.22 and S^/dg - 4.0 from reference (16), 
the corresponding values of effectiveness at S-j/dg - 3.0 were compared 
to those at 4.0 at each x'/dg from 10 to 80. Generally, at S-j/dg - 
4.0 was between 14% and 20% lower than at S-j/dg ” 3.0. This correction 
factor, transforming the value of at S-j/dg «* 3.0 to S-j/dg = 4.0, 
can be incorporated into Eq. (112) in the following manner 


q FC .V" 0 


^0 


‘0 


ADW 


s* e:] 


(113) 


where S* - correction factor to adjust from S^/dg = 3.0 to 
4.0; a function of distance from injection, x'/dg. 

The present investigation was conducted with a coolant-to-gas 
density ratio of 2.18. With the gas and wall temperatures remaining 
nc^/ly constant (Tj /T w ^ 2.15), Gg was typically 1.03. For P c /p ra = 
2.18, values from Pedersen were interpolated between Pq/p w ~ 2.0 
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and 4.17. Equation (113) then takes the form 


a » ^0'=Q 

q FC „ ,_°C u 


"io * ( * 


0 


•)[1 - 1.03 (n ADW ) S*] 


2.18 


(114) 


where 


( n ADlP 


2.18 


= adiabatic wall effectiveness evaluated at 


P c / Po = 2.18. 


Eriksen (20) reports values for the hydrodynamic effect (h Q i-q/Hq) 

C 

for x’/d Q up to 100 and M ranging from 0.2 to 2.0. However, these 

results were for streamwise injection at 35° with S-j/dg = 3.0. Since 

no information was available concerning the hydrodynamic effects for 

S-j/dg - 4.0, h 0)= Q/h o was taken as 1.0 in the use of Eq. (114) to trans- 
C 

form Pedersen’s data. 

Figure 101 was prepared using Eq. (114) in conjunction with Pedersen' 
adiabatic wall effectiveness data. Figures 101(a), (b) and (c) contrast 
the reduction in heat flux due to film cooling for Case V for the subject 
investigation (solid curves) to the transformed values using Pedersen's 
data corrected to S^/dg = 4.0 (broken curves). Streamwise injection at 
35° to the surface is represented by the broken line, while spanwise 
injection at 35° is shown as the solid line. Figures 101(a), (b) and 
(c) are for a blowing ratio of 0.213, 0.515 and 1.05, respectively, for 
x'/dg ranging from 10 to 80. 

From the comparison in Fig. 101 it appears that the streamwise 
injection configuration is more effective in reducing the local heat flux 
than the spanwise injection scheme. At M = 0.213 (Fig. 101(a)) stream- 
wise injection shows qp C /qjj values approximately 1% lower at x'/dg = 10. 
However , at x'/dg = 80 the difference is less than 2%, At this low value 
of blowing ratio, the decay of the film cooling effectiveness (i.e. rate 
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of increase of qp C /qg with x'/dg) is approximately the same for both 
configurations. 

At M = Q..515, shown in Fig. 101(b), the value of q^/q^ is nearly 
the same for a ~ 35° and $ = 35° at x'/dg = 10. However, the rate of 
increase of qp C /qjj with x'/dg is more severe for spanwise injection at 
this blowing ratio, especially in the range x'/dg of 20 to 40. At the 
farthest location from injection, streamwise injection is still 4% better 
than the spanwise configuration. 

Finally, Fig. 101(c) contrasts the two configurations at M = 1.05. 

At x'/dg = 10, the two injection schemes yield essentially the same 
values for q^/qg. Again the rate of increase of Op C /qg with x'/dg for 
a = 35° is more pronounced than for 3 = 35° but roughly the same as was 
noted for M = 0.515. For spanwise injection, qj^/qg is nearly 7 % higher 
at x'/dg = 20 and 11% greater at x'/dg = 40 compared to injection with 
a = 35°. These differences are somewhat higher than at the lower blow- 
ing ratio of 0.515. At x'/dg = 80 the streamwise configuration is still 
7% more effective than spanwise injection. 

It should not be surprising that spanwise injection is less effec- 
tive in reducing qp C /q[) than streamwise injection. With no initial com- 
ponent of the coolant momentum in hhe freestream direction, spanwise 
injection would probably create a greater flow disturbance near injec- 
tion. It is noteworthy that spanwise injection appears to be nearly as 
effective in reducing qjig/qg near injection (i.e. x'/dg « 10). Overall, 
the two configurations do not appear to differ more than 11%. For lead- 
ing edge cooling, where curvature may be a constraint on the minimum 


streamwise angle that can be incorporated, the penalty for spanwise 
injection (as illustrated in Fig. 101) may not be as severe as would be 
incurred for steeper streamwise injection angles. 
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VIII. CONCLUSIONS 

The results from an experimental investigation of film cooling in 
the leading edge region of a turbine vane model were presented. A single 
row of holes, angled s panwise to the main flow at either 18° or 35° to 
the surface, was used to inject cool air ( pq / p ^ = 2.18) into the suction 
surface boundary layer at x../L s values or 0.0833, 0.150 or 0.228. The 
effects of cool mass addition on the reduction in the local Stanton 
number were measured for the six combinations of injection angle and 
location. 

During the entire investigation, all flow conditions external to 
the test surface were maintained at nearly constant values and corresponded 
to a simulation of a high temperature and high pressure environment. Gas- 
to-wall and coolant-to-gas temperature ratios were held nearly uniform 
with average values of 2.15 and 2.18, respectively. The inlet turbulence 
intensity remained at 6.655 during the course of the film cooling study. 

From the analysis of the results the following conclusions can be 
drawn concerning the effects of film cooling on the value of STp^/STg for 
the investigated configurations: 

1. For film cooling injection in the leading edge region of a tur- 
bine vane at injection angles of 8 = 18° and 8 = 35° (Cases I - 
V), the local film cooling heat flux was reduced compared to 
uncooled conditions (i.e. STp C /ST 0 < 1) by increasing the blowing 
ratio, M. 
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The value of M at which the greatest heat flux reduction occurred 
(HqPT was a function of both the angle of injection, B, 
and the location of injection, x../L s : 

a. Decreasing the angle of injection reduced the effective 

normal component of the coolant momentim, I sin 3 (where 
P p 

I - (— ) Nr). The use of the shallower injection angle 
P C 

permitted a larger film coolant blowing ratio, M, before 

the film cooling performance began to decrease due to 

excessive mixing. The critical value of blowing ratio, 

Mqpt was observed to be 0.85 and 1.05 for B ~ 35° 

and 18°, respectively (for moderate values of freestream 

-fi 

accel eration* K- - 1 x 10 ). 

b. Injecting closer to stagnation (i.e. into a thinner boundary 
layer) caused Mgpj to increase. It is expected that, 
for a constant injection angle and blowing ratio, the 
coolant will remain nearer to the wall for injection into 

a thinner boundary layer. Steep injection angles allow the 
coolant to mix more freely with the mainstream, resulting 
in decreased cooling effectiveness. For injection into a 
thinner boundary layer, a larger value of M is permitted 
before the cooling performance begins to deteriorate due 
to mixing. For B - 18°, the following trend was noted: 
Optimum Blowing Momentum Boundary Layer 


Ratio, H 0 p Tt NEAR 

e i^ d 0 

1.3 

0.013 

1.05 

0.036 

C .5 

0.051 



273 


A similar trend was noted for 3 = 35° , however, due to the 
possible influence of local freestream acceleration at this 
steeper angle, the full impact of changing G./dg could not 
be assessed. 

The increase in M Qp j due to decreasing both $ and 
9.j/d 0 can be of benefit to a designer who may be concerned with 
hot gas ingestion into the coolant holes. In the leading edge 
region, the difference between the local hot gas pressure and 
the coolant supply pressure may not be large. Consequently, 
small variations in the total gas pressure or coolant supply 
pressure could lead to back flow through the coolant holes. 

The subject investigation has indicated that optimum cooling 
in the leading edge region is achieved with high values of the 
blowing ratio which may be helpful to insure a safe margin in 
pressure between the coolant and hot gas. 

3. Values of M corresponding to a minimum in STp C /ST Q support the 
results of Pedersen (24) which indicate that the velocity ratio 
V c /V ro , may play an important role in determining how effective 
a given amount of coolant can be in reducing the local heat 
flux. Values of Vg/V^ between 0.4 and 0.5 have been reported 
(24) to yield the highest cooling effectiveness for all density 
ratios, pg/p^, for streamwise injection. For 0 = 18° at Q^/d^ 
0.013 and 0.036, the corresponding (Vg/V ro ) values for the 
observed Mgpy were 0.596 and 0.432, respectively. 
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4. For moderate levels of the frees tream acceleration at injection 
(i.e. K. * 1 x 10“ D ), decreasing the injection angle from 3 = 

35° to 3 = 18° appears to have a minimal effect on STp^/ST^. 
Except for extending the value of M Qp ^ to a higher level 
(with a corresponding improvement In STpQ/ST Q ), there is little 
difference between the cooling potential for both angles up to 
Mqpj for $ - 35°. This result may have significant 
economic implications. For regions of the vane with moderate 
and low freestream acceleration, the designer need not specify 
extremely shallow angles (i.e. <35°) to achieve reasonable 
cooling. Shallow angles are costly. Tooling becomes more 
intricate and tolerances more difficult to hold as smaller 
angles are specified. A considerable savings in fabrication 
costs may be realized if the designer is able to identify 
regions of the vane in which the local reduced heat flux due 

to cooling will be essentially insensitive to differences in 
3 below a certain threshold level. 

5. Large values of the local freestream acceleration parameter, 

“6 

(i.e. K. > 1 x 10 ) may reduce the cooling effectiveness 

markedly under certain injection conditions. The effects of 

acceleration appear to be coupled to the angle of injection. 

For injection at x^/L s - 0.0833 and K. = 7.9 x 10“ 6 (i.e. Case 

IV), the Stanton number ratio was as much as 48% higher for 

3 = 35° compared to 3 = 18°. At x^/L s = 0,150 and fCj " 0.72 x 
-6 

10 , little difference was noted between the Stanton number 

ratio for both angles. It appears that by specifying as shallow 
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APPENDIX A. rams AND CORRELATIONS i'O^ FILM EFFECTIVENESS 
A . 1 . Analvti cal Model s 

A very thorough discussion concerning the differences in film effec- 
tiveness predictions for 2-0 slot- type models is contained in a review 
of film cooling by Goldstein (44). Familiarity with those cooling para- 
meters which result from those modeling analyses is essential to the 
understanding of the vest library of film cooling daw. and its extension 
to practical applications . The following discussion on 2-D slot models 
is summarised from Goldstein (44), 

The effectiveness of an injected coolant has been evaluated analyt- 
ically by a number of investigators using simplified models to solve the 
integral energy equation. Meat sink models that have been developed 
differ considerably due to the assumptions concerning the interaction of 
the coolant and the mainstream boundary layer at injection. The models 
of Libriszi and Cresci (65) and Kutateladze and Leont'ev (65) both 
assume a uniform boundary layer temperature, T, and account for mass 


addition into the boundary layer at injection. Tribus and Klein (63) 
use a representative temperature profile but do not account for mass 
entrainment at injection. Stollery and El-Ehwany (67) neglect mass 
addition into the boundary layer at injection and assume an average 
boundary layer temperature, T. Goldstein and Haji -Sheikh ( 63 ) incor- 
porate both phenomena of a reduced temperature profile in the boundary 
layer as well as mass addition and predict effectiveness values between 
(65) and (63). 
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The heat sink model analysis predicts that the adiabatic cooling 
effectiveness., can be expressed 
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(A-l) 
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where m 0 ~ entrained freestream mass flow 

m c ~ coolant mass flow 

C = specific heat at constant pressure of the frees tream 
Pco 

C = specific hoat at constant pressure of the coolant 
P C 

The essential differences between all the models stems from how 
m o3 /!Tic is handled. 

Librizzi and Cresci assume the point of injection and the initial 
location of the boundary layer development are coincident. Using a 1/7 
power velocity profile they found 

mjm c » 0,329 (x7M5)°° 8 [Re c (y c /iiJ3"°* 2 
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Kutateladze and Leont'ev assume that the boundary layer starts 
upstream of injection and that its mass flow at injection equals that 


*ttote subscripts on n refer to initials of investigators responsible 
for development of the particular expression. 
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of the coolant. Such an analysis yields 
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Stol'iery and El-Ehwany assume that the total mass flow within the 
boundary layer at the point of injection is zero. This analysis yields 
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Tribus and Klein, neglecting any effects of injection on the main- 
stream flow, modeled a cooling slot using a line heat source. They 
found 


n T „ - 5,76 Pr 2/3 {C r /C„ ) q" 0,S 

1 K Pc p» 


(A-7) 


Goldstein and Ha ji -Sheikh accounted for the temperature variation 
through the boundary layer and for the effects of injection on thicken- 
ing the boundary layer to predict an effectiveness for tangential injec- 
tion of 
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Figure A1 shov/s the vaviation in the predicted coding effectiveness 
as a function of the dimensionless distance from the cooling slot for all 
the models previously described. The curves represent M =0.5 and 
Re^ = 1000. Lribrizzi and Cresci along with Ku tatd adze and .Leont'ev 
both assume an average boundary layer temperature while accounting for 
frees tre am mass entrainment. They predict the lowest effectiveness. 
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Seban, Chan, Scesa (69) 

(Experimental) 
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Dimensionless Distance From Slot, ::/S 
Figure Al. Comparison of Film Cooling Models £2~D) 
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Tribus and Klein use a representative temperature profile but do not 
account for entrainment. Their prediction is the highest. Goldstein 
and Haji -Sheikh incorporate both phenomena and predict effectiveness 
values between the other two groups. Obviously a neglecting mass addi- 
tion at injection is too liberal an assumption, while assigning an aver- 
age temperature to the boundary layer is too stringent. Actually the 
two effects are counterbalancing as shown by the agreement of the pre- 
dicted n between Stoll ery and El-Ehwany {neglecting entrainment and 
assuming T in t^e boundary layer) and Goldstein and Haji -Sheikh. An 
empirical correction from the data of Seban s Chan and Scesa (69) is 
superimposed on the analytical model comparison of Fig. Al, Their 
results are in good agreement with Goldstein and Haji-Sheikh as well as 
Stoll ery and El-Ehwany, 


A. 2. Semi -Empirical Relations 

Hatch and Papell (70) describe a model ^'n which the coolant is 
assumed to remain as an unmixed protective film along an adiabatic wall. 
The effectiveness is described:, using an empirical correlation for the 
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= convective external heat transfer coefficient 
= length of coolant slot 
= flow rate of coolant 

- specific heat at constant pressure of coolant 
= thermal diffusivity of coolant 


Seban and Back (71) describe the slot injection phenomenon as being 
similar to a free jet boundary which becomes altered by the proximity 
of the wall. For thin initial boundary layers at injection* an initial 
mixing region exists near the slot having velocity profiles similar to 
those found in a jet boundary flow. The extent of the mixing within 
this region is strongly dependent on the ratio of coolant to mainstream 
velocity. Far from injection (x'/$ > 40), the boundary layer res sables 
the typical turbulent pov/er law profile. Assuming uniform eddy diffu- 
sivity consistent with experimental results and similarity of the tem- 
perature profiles, an effectiveness relation was derived of the form 
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where Re s - Reynolds number of coolant based on slot height, S 

3* = (x + x 0 ), distance downstream from effective origin or 

boundary layer 

EL - (x * X]), where x-j is distance from slot to where 
effectiveness was unity 


Spalding (72), independently using the same model as Lihrizzi and 
Cresci but accounting for entrainment into the boundary layer, developed 
an effectiveness parameter, g 
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This expression is valid for ell velocity ratios, V^/V^. X has been 
adjusted as shown with properly chosen constants to fit a vast group of 
experimental data ranging from V^/V to = C to 14. 

Ass lining that the boundary layer starts from the point of injection 
and that the thermal and hydrodynamic boundary layer thicknesses are 
equal, Nishiwaki, Hlrata and Tsuchida (73) derived an expression for 
cooling effectiveness of injection through a porous strip along an 
adiabatic plate. 

” ^*0877 B (16 - 11 mo) m3 


(A-17) 


250 


for 0 < m5 < 1 

whore (mS) is the root of the following 

(m6) 3 - 1.5 (mS) 2 + 0.0417 + = 8 (A-18) 

and B = (t) Re®;® (A-19) 

q = weight flow of cool ant per width of plate 


Eriksen, Eckert and Goldstein (74) present a three-dimensional film 
effectiveness model in which cooling from a discrete hole is simulated 
by a moving point source in a semi-infinite medium. An expression for 
the effectiveness downstream of the source can be shown to be 
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where e ~ point source thermal diffusivity 

r 

dg = hole diameter 

z* ® lateral distance from centerline of hole directed 
along X axis 

y Q = horizontal virtual origin of source above X-Z plane 

This analytical model relies on empirical input for both the thermal 
diffusivity and y^. At low blowing ratios (M ~ o.l) the model agrees 
well with their effectiveness data (74) for a single hole. At blowing 
ratios greater than 0.1 and less than 0.5, the model can be useful 
provided reasonable values of y^ are inserted to account for separation 
of the coolant jet from the surface. Provided negligible jet inter- 
action is assumed, the superposition principle can be used with the 


model to predict the additive effect of a row of holes or other rela- 
tively complex injection configurations. 


A, 3. Experimental Correlations 

Consistent with the predicted governing film cooling parameters 
derived from a heat sink analysis, a large portion of film cooling data 
can be represented in the following form where A, b, and c are constants. 



(A-21) 


The next few paragraphs will deal with a few studies that were able to 
describe their data in such a way. All the correlations apply to slot 
type cooling with essentially negligible differences between the coolant 
and mainstream temperatures. 

Meighardt (75), in his pioneer work on aircraft wing de-icing, 
found that the film heating data correlated well with the following 
effectiveness expression for a slot at 3(1° to the surface. 
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Hartnett, Birkebak and Eckert (3). show a similar effectiveness 
dependency on the blowing ratio for a 30° slot and have correlated their 
results for x'/MS 2 . 60 by using 
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Scesa (76) conducted a study of normal injection that supported the 
use of the adiabatic wall temperature approach in calculating heat 
fluxes. Performing both heat transfer and adiabatic wall temperature 


tests, development of the expression 


n ADM 



for M >0.3 


(A-24) 


proved useful in heat transfer calculations by use of the adiabatic wall 
temperatures for downstream locations not near the injection site. 1 

Seban, Chan, and Scesa (69) provide a tangential slot correlation 
that shows a somewhat different dependency of r\^ on (x '/MS) and Re c 

p 

compared to the previous correlations. They found 
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Samuel 

and Joubert 

(78) present the following relation, 

valid for 

M less than 

unity and tangential slot i 

njection 
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for 40 < (^r-) M" 1 * 75 < 400 

Metzger and Fletcher (10), using a transient heat transfer measuring 
technique, calculated local heat transfer rates at a wall with and with- 
out coolant injection. The ratio of the average heat flux to the wall 


1,2 


Surraiarized from review by Metzger (77). 
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with injection to the heat flux without injection, was found to be a 


function of the mass velocity ratio, temperature difference ratio s£, 
x'/3 5 and the geometry. The value of at S’ 1 « 0 corresponds to the 
inverse of the standard cooling effectiveness parameter. 
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They found for angled slots: 


(A- 29) 
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for 0.25 £ M < 0.75, 0 < x7S < 70, % » 20° 



and 
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for 0.25 < M < 0.75, 0 < x'/S < 70, © = 60° 
where a = injection angle. 

The list of existing empirical correlations could go on and on with 

the band of predicted effectiveness values continually widening. Figure 

A2 makes a simple comparison of the previously presented correlations at 

a nominal blowing ratio of M =0.5 and a coolant slot Reynolds number of 
4 

10 . The effectiveness is presented as a function of non dimensional dis~ 
tance from injection. The normal injection of Scesa is understandably 
low compared to the others. For 30° slot injection (Hartnett, Birkebak 
and Eckert) the effectiveness improves from the normal configuration. 

At 20° injection (Metzger and Fletcher) the effectiveness shows even 
more improvement. There is much inconsistency, however, even in the few 
correlations that are shown, Wieghardt's 30° slot data seems more 
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Figure A2. Comparison of Experimental Film Cooling Correlations (2-D) 




effective than the tangential cooling results from Satan, Chan and Scesa. 
Samuel and Joubert's tangential data is shown as less effective then 
Metzger and Fletcher’s 20 c correlation. 

There is much discrepancy among the various empirical correlations 
in predicting the film cooling effectiveness for a given injection geo- 
metry. As shown in Fig. Al, the analytical model predictions also show 
the same inability to universally characterize the cooling effect of 
injection from a specified geometric configuration. This brief review 
of some of the more important cooling correlations was introduced solely 
as a means to familiarize one with the parameters that characterize film 
cooling in general and to emphasize the need to exercise caution in 
applying film cooling data to real problem situations. Mo one correla- 
tion is recommended as the most suitable for a given application. At 
present, with such previous information, one can, at best, make only 
approximate calculations to determine the benefit of a proposed cooling 
scheme. 
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APPENDIX B. EXPERIMENTAL APPARATUS DETAILS 
B.l Tu rbul en ce J3en epajt ij;j^ As_s esnbljf 

To estimate the inlet Tu variation by locating grids of different 
sine at varying upstream locations, the theory of isotropic turbulence 
decay downstream of screens was incorporated (45,46). This theory 
reveals a linear dependency of the inverse square of the intensity with 
the distance downstream from the generating grid, x. 




(B-l) 


where Tu = '7 r “/V 

CQ r CO 

u* = fluctuating component of freestream velocity 

7 m - time average of freestream velocity 

c - constant 

C D “ drag coefficient of screen 

(#-> = value of (ff) at Tu „ ” 0 

" 0 *> 

M s - center- to- center spacing of grid wires 

The following constraints were imposed on Eq. (B-l) in designing 
the turbulence generation system: 

(1) x/M c >.15 for measurements to insure isotropy; 

W ' 

(2) Mg/dy " 2, set as a lower limit for establishing isotropy at 
the shortest distance downstream of the grid; 
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(3) turbulence generation assembly must be as compact as possible* 

(4} the point of intensity measurement downstream from the grid 
assembly must be 1.5 to 2 characteristic model dimensions up- 
stream of the test surface to insure that the presence of the 
model will not affect the value of the intensity at the measur- 
ing location. 

An overall working length for the assembly was set at 0.533m (21 in), 
corresponding to a practical length for a relatively low, estimated tur- 
bulence level. This overall length was then subdivided into 3 sections 
of 0.46m (1,31 in), 0.065m (2,54 in) and 0.381m (13,84 in). Figure B1 
shows an exploded view of the three sections in one arrangement relative 
fco the point at which a hot wire probe can be inserted. The 0.046m 
section serves as a screen holder and allows either of two grids to be 
positioned at eight different locations upstream of the measuring site. 

The three stainless steel rectangular ducts of the turbulence 
section measure 0.114m (4.5 in) in the lateral direction and 0.244m 
(9 >6 in) in height. The main flow channel prior to the turbulence pro- 
duction assembly consists of a 6.096m (20 ft) straight circular 0.254m 
(10 in) OD pipe with a slight contraction in the lateral direction when 
mating to the rectangular turbulence section. 

B.2. Instrumentation 
B.2.1, Mass Flow Measurement 

Orifices, constructed to ASME standards, were the prime means of 
mass flow measurement for the basic flow system. The main by-pass air 
was metered through a 0.0737m (2.9004 in) diameter orifice having a 
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£' - 0.3401 (3* - ratio of orifice diameter to pipe diameter). The com- 
bustor was mixed with supply air from a 0.0598m (2,3533 in) venturi with 
a 8' = Q.57S6 and natural gas through either of too orifices. For low 
flow fuel, a 0.0209m (0.824 in). diameter orifice with 3' = 0.4238 was 
used, while higher temperature conditions required using a 0.0266m 
(1.049 in) diameter orifice with a 0* - 0.2852, 

Downstream of the test section, the flowrate through each leg of 
the cascade was monitored with orifice plates. The upper and lower by- 
pass legs had orifice diameters of 0.076995m (3.0318 in) and 0.076990 
(3.0311 in) and 3' of 0.49980 and 0.49977, respectively. The main test 
channel flow passed through a 0.1268m (4.9933 in) diameter orifice with 
a 3* = 0.62564. 

Film coolant flowrate measurements were handled in three ranges; 
low, intermediate, and high. For cooling at low blowing ratios, 0,00102m 
(0,040 in) diameter venturis with 3' - 0.1333 measured the mass flow to 
both the suction and pressure surface coolant plenums. An intermediate 
set of venturis with 3' = 0.2667 and diameters of 0.00204m (0.080 in) ' 
covered the widest range of coolant flows. For blowing ratios (M) from 
1.25 to 2.0, small Q. 01 27m (0.5 in) diameter turbine flowmeters handled 
from 0.708 to 3.305 liter/sec (1.5 to 7 cubic feet per minute), 

R.2.2. Test Surface Temperature Measurements 
Wall temperature measurements were made using 0.00102m (0.04 in) 
and 0.00127m (0.05 in) diameter, swedge-type copper cons tan tan thermo- 
couples. A small portion of the magnesium oxide insulation was removed 
from the stainless sheathing, exposing the copper and constantan 
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conductors. A stainless tube with an outside diameter slightly smaller 
than the sheath inside diameter was placed around the conductors within 
the outer sheath. This arrangement was he! i -arched, the insert tubing 
providing a filler material for the evacuated space and serving as a 
conductive medium for the thermocouple junction formed during the arcing 
process. The wall thermocouples were furnaced brazed at 1089°K - 1144°K 
(1960°R - 2060°R) into designated hole locations on the test surface and 
were finished flush with the wall. Exposed-ball or bead-type thermo- 
couples were used in the film plenum cavities. The exposed couples were 
positioned with the junctions at the center of the plenum, indicating the 
characteristic mixed film temperature. 

B.2.3. Test Surface Pressure Measurements 
Static pressure measurements were made with flush mounted 5.08 x 
10~\i (0.02 in) copper tubes brazed into the copper test surfaces. 

Plenum pressure ports consisted of 5.08 x 10"\ (0.02 in) holes drilled 
through to the plenum and sealed on the backside with 0.00102m (0,04 in) 
0D stainless extension tubing. 

B.2.4. Inlet Conditions Measurements 
To measure the inlet parameters, three, wedge-shaped, pitot-static, 
total temperature probes were used: two to traverse the inlet conditions 
and the remaining probe to monitor the flow conditions at midstream 
between data points. 

Each probe had a head diameter of 0.0043m (0.189 in) and a total 
length of 0.61m (24 in), 0.559m (22 in) of which incorporated a 0.00635m 
(0.25 in) diameter reinforcement tube. A chromelalumel thermocouple, 
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operational to 811°K (]46Q°R)> as well as the pitot tube were located at 
the apex of the wedge with a static tap on each side of the probe head. 
Due to the unusual length of traverse, 0.244m (9.6 in) vertically and 
0.114m (4.5 in) laterally, the probes were permanently part of a 
0.9127m (0.5 in) diameter support sleeve. L. C. Smith actuator systems 
were used to drive the three probes. 

Figure B2 demonstrates the probe measurement system. A differential 
pressure bet ween both static pressure taps gave an indication of the 
probe alignment. These two pressures were also fed to a United Sensor 
Corporation pressure dividing choke that averaged the nearly balanced 
pressures. The average static condition was referenced to the measured 
total pressure for a dynamic pressure reading, while the total pressure 
was read directly against the atmosphere (gage pressure). Factory 
supplied calibration data showed the total temperature and pressure and 
static pressure were accurate to within 1/255 up to a Mach number of 0.40. 

Prior to the initiation of the actual test program, measurements of 
the cascade inlet turbulence levels were made over a range of inlet 
Reynolds numbers and turbulence-generating grid-locations. 

Figure 133 illustrates the electrical network used in measurement 
of the cascade inlet turbulence intensity, A Thermo -Systems Incorporated 
(TSI) 1220-20 hot film probe was used in conjunction with a TSI 1054 A 
linearized anemometer for measuring the raw turbulent signal. The 
anemometer output was fed to both a Fairchild digital voltmeter and an 
oscilloscope for simultaneous monitoring of the average and instantaneous 
probe output. In addition, the anemometer output was fed to both a 
Hewlitt-Packard 3400 A RMS voltmeter and a DISA RMS voltmeter for 


303 



Figure B2. 

Pressure Measurement Set-Up for Pitot-Static Probe. 
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comparison purposes. Direct monitoring of the RMS signal was also 
obtained using another digital voltmeter. Permanent records of the 
average and RMS turbulent signals were made on a Visicoder, Model 1108, 
high response oscillograph. 


B.3. Measurement Systems 
B.3.1. Raw Data Conversion 

All pressure measurements were made using two methods: strain-gage 
type transducers or a digital automatic multiple pressure recording 
(DAMPR) system. With DAMPR, the pressures to be measured were all 
brought into a pressurized chamber, each separated from the main chamber 
by a pressure sensitive switch. Once the chamber was opened to atmo- 
spheric conditions, the chamber pressure decreased and the time to reach 
equilibrium was recorded with a digital counter. As the chamber pres- 
sure fell and reached the input levels of the multiple unknown pressures, 
the pressure switches opened and recorded a digital opening time with 
respect to the initial chamber depressurization. A comparison of the 
time at which each switch opened with a calibration of "time-to-equi lib- 
rium" for the depressurized chamber gave an accurate reading of each 
unknown inlet pressure. 

Hie raw millivolt signal from each thermocouple was referenced to 
a 338°K (610°R) oven. 

The output from the heat flux sensors was brought up to the milli- 
volt range using operational amplifiers with gains near 1000. As 

mentioned previously, each gage-amplifier set was individually calibrated 

-3 2 

with the gain adjusted to give 8.64 * 10 mv output for each joule/sec-m 


T . A 
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of incident heat flux, Preston, Model 8300, wide band, floating, differ- 
ential amplifiers were used as a buffer system to provide both additional 
gain capability as well as remote zero-setting for the heat flux opera- 
tional amplifiers. 


B.3.2. Data Monitoring 

Excluding pressure measurements under the DAMPR system, there were 
110 basic measurable quantities in this experiment, generating charac- 
teristic output voltages of 0 to 1000 millivolts. The test facility 
was equipped with a VIDAR, Model 5218, integrating digital voltmeter 
and a VIDAR, Model 610, low level scanner to continuously monitor these 
110 channels. The scanner-voltmeter combination simultaneously scanned 
the 110 channels every 5 seconds. 

Following a complete scan, the VIDAR unit transmitted the infor- 
mation from the 110 channels to a Digital Equipment Corporation, Model 
PDP-11, minicomputer. A tape-loaded, data reduction program converted 
the millivolt VIDAR signals to engineering units and proceeded in cal- 
culating basic flow parameters. The output routine from the PDP-11 
was linked to both a Setchell-Carl son. Model 10M915, and a Tektronix, 
Model 4010-1, visual display system. 

After each scan by the VIDAR unit, the PDP-11 sent updated infor- 
mation to the Setchell-Carl son display unit. The following items 
appeared after each scan and were used by the test operator to set the 
desired flow condition at the test section: 

Weight Flows: 

...by-pass leg 
...combustor leg 
...fuel flow 
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...test section, all three channels 

...film coolant flow, suction and pressure vane 

Temperature: 

...total inlet 
...film coolant plenum 

Pressure: 

...total inlet 

Inlet Mach Number 

Inlet Reynolds Number 

Exit Mach Number of Suction Vane 

The output routine was interrupted after a specified number of 
scans and was transferred to the Tektronix Visual display. More detailed 
information concerning the conditions within the test section was dis- 
played on this unit. A permanent record of the Tektronix display was 
made using a Model 4610 hard copy unit. A sample copy is shown in 
Fig. B4. The Tektronix unit was a convenient means for recording test 
section conditions prior to each data point. 

B.3,3. Data Acquisition, Storage and Processing 
Once a data call has been placed, the VIDAR scanned channels and 
the DAMPR signals were simultaneously sent to a central processing 
center where the digital data was recorded on tape and sent to the 
computer for temporary storage in a program data set. 

As with the PDP-11, a data reduction program was permanently filed 
in the library of the center’s IBM Model 360 Computer. Because the 
IBM 360 unit was a tune sharing system (TSS), the data was processed 
from any of the numerous remote IBM 2741 terminals. The test facility 
had both a remote terminal and a Versatech, Model 200A, matrix 
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Figure B4 . Sample of Permanent Record of Test Conditions 
From Textronix Visual Display Unit. 
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printer/pl otter. Any data that was taken could be immediately reduced 
using the TSS system. 

The printed output could be displayed on the 2741 terminal or on 
the Vers a tec. 


B.4. General Test Procedure 

The recording of a typical data point began by setting the overall 
flow conditions within the test section. With a pitot-static, total 
temperature probe positioned at midchannel , just upstream of the test 
surfaces, the total temperature, total pressure, exit Mach number at 
the suction vane, and test section flow-splits were adjusted for proper 
reduced-condition simultation of a high pressure, high temperature 
envi ronment. 

The film coolant flow rate was set and a permanent record of the 
test conditions was taken with the hard copy unit. The set-point probe 
was then fully retracted from the stream, leaving a clear channel. 

The heat flux levels at the test vanes were then recorded, on a 
Visicoder, Model 1103, high response oscillograph, in the absence of 
any upstream probes. This served simply as back-up data in the event 
that some malfunction occurred in the VIDAR system. 

The next step in the data taking procedure was the actual test call. 
However, a brief description should be made of the sequence of events 
that constitute a data call. The call its.elf initiated a sixteen-cycle 
process wherein the vertical and horizontal pitot-static, total tempera- 
ture probes mapped the inlet conditions. Cycle one began with the 
vertical probe fully extended to the top of the test channel, at which 
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po’nt the VIDAR unit scanned the first 68 input channels, sending the 
signals to the central processing center. The vertical probe then 
stepped down the channel, stopping automatically every 0.0152m (0.6 in) 
for a 68 channel scan* At the ninth cycle, the horizontal probe actuated 
0*O127m (0.5 in) into the stream. The two probes traversed together the 
last 8 cycles until, after the 1 6th cycle, the vertical probe was 
entirely retracted and the horizontal prrte reached the opposite wall. 

The 16th cycle was a special scan in that all 110 channels were sampled. 
The computer program then averaged the sixteen, cycle signals from the 
first 68 channels and used til is information with the 16th cycle signals 
from channels 69-110 to perform the required data reduction. Figure B5 
is a diagram of the 16-cycle, probe-actuation process. The same pro- 
cedure can be accomplished without the probes in the stream by manually 
signaling VIDAR that the probes were in the proper location to execute 
a scan. This, of course, requires simulating the probe actuation with 
16 manual calls. 

Continuing, then, to describe the data taking procedure, sixteen 
manual cycles were recorded, followed immediately by 16 automatic cycles 
using the probes as previously described. The reason for the dual nature 
of the cyclic data was due to the influence of the inlet probes on the 
measured heat transfer to the test surfaces. The manual cycles were 
recorded with no probe influence. The heat flux data obtained from 
this manual procedure was temporarily stored. Immediately following, 
the automatic cyclic data was sent for reduction with the temporarily 
stored heat flux signals from the previous manual process substituted 
one-to-one into the automatically actuated cyclic data. In essence, 
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two data points were required for processing a true test condition. The 
printed output data indicated from which reading the measured heat flux 
was chosen. 

The horizontal set-point probe was returned to midstream following 
the probe cycling. Any changes in the flow conditions could be made and 
the procedure was ready to be repeated. 


f 
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APPENDIX C. HEAT FLUX SENSOR PERFORMANCE 

Generally) the miniature, Gardon heat flux transducers used in the 
subject investigation functioned satisfactorily; providing film cooling 
effectiveness information while operating within a severe thermal envi- 
ronment, The nature of the Gardon gage operation {as previously de- 
scribed in IV. 1.4. 5.) allowed heat transfer measurements to be made 
under conditions of large temperature differences through the boundary 
layer (i.e. gage at cool, wall temperature with hot mainstream) , This 
is one important requirement for modeling high temperature and pressure 
gas turbine vane heat transfer. The dimensions of the subject heat flux 
transducers were small in comparison to the test surface characteristic 
length, providing, essentially, point measurements while offering little 
disturbance to the external flow and boundary layer development. 

During the course of the subject investigation some difficulty was 
encountered with the operation of a few gages, especially on the pressure 
surface of the vane channel. It became exceedingly more difficult, as 
time passed, for these gages to reproduce the heat flux reference con- 
ditions necessary to qualify the film cooling data. 

Recall for a moment the data qualification criterion. Heat flux 
measurements under film cooling conditions were referenced to the uncooled 
conditions (i.e. q^/q^) to define a cooling effectiveness. If changed 
more than 3 % during the time the cooling data were taken, no cooling 
effectiveness values could be calculated and the data was discarded. 
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During a period following the dry-wall heat transfer measurements, qjj at 
each gage location on the pressure surface was very difficult to repro- 
duce. 

Figure Cl illustrates the variation in the reference heat flux level 
(as characterized by STq for the pressure surface as a function of 
time during a run and as a function of day to day operation. The verti- 
cal bands represent the variation in STq ^ (i.e. qjj reference) during a 
given day. Four different days during which film cooling heat transfer 
was measured are illustrated by the bands labeled a, b, c, and d. Gen- 
erally, the majority of the reference qjj levels (bands at each of the 
various Re^ represent output from a single gage; 6 gages are represented 
in Fig, Cl) varied much more than 3 % during the time interval in which 
film cooling data were recorded. Consequently, few qpe/q|J values satis- 
fied the data qualification criterion. 

In contrast to the pressure surface, a similar display of the varia- 
tion in the reference conditions is shown in Fig. C2 for the suction 
surface for the same time Interval as Fig. Cl. Host of the bands of 
STq x variation for the suction surface are within 3 % of a mean value, 
thus satisfying the qualification condition. 

By examining both Figs. Cl and C2, a trend of an increasing average 
reference condition can be detected (i.e. mean value of bands increases 
as one proceeds from “a" to "d"). This phenomenon is more pronounced 
for the pressure surface than the suction surface. The reason for this 
occurrence was found to be a sensitivity change in the basic sensor. A 
recalibration of the sensors on the pressure surface verified that some 
gage sensitivities had changed by a factor of 2 over an extended 
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exposure period to the hot combustor gages. A spectrographic analysis of 
a residue which formed on the test surface during a run time of extended 
duration (cleaning of this residue from the test surface using a non- 
corrosive solution was a regular practice following each test) offered 
a plausible explanation for this sensitivity change. 

Large amounts of copper as well as traces of clorides and nitrates 
were detected in the residual specimen. The presence of nitrates sug- 
gests that nitric acids may have been forming in the combustion products, 
creating a corrosive environment at the test vane surfaces. Since the 
heat flux gage sensitivity is inversely proportional to the foil thick- 
ness (48), the gradual decay cf the thin foil (2.54 * 10" 6 m, 0.001 inch) 
resulted in a much larger AT across the foil for a constant incident 
heat flux. 


Fortunately, film cooling effectiveness data was the primary infor- 
mation to be obtained and was of a relative nature. Therefore, if the 
sensing device was linear and remained linear with time, even though the 
sensitivity changed (as was verified by the recalibration of the pressure 
surface gages), the cooling effectiveness would be unaffected. This is 


true, however, only if the time interval between measurements of 


and 


qg for a particular set of coolant conditions is much smaller than the 


time for the gage sensitivity to change significantly. Consequently, all 


data that were taken in which a film cooling data point was immediately 


followed by a reference reading were determined to be accurate due to 


the short time interval between data points relative to the time for a 
sensitivity change. For the cases in which many film cooling conditions 
were set between reference checks, only that data was acceptable in which 
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the reference points, from beginning to end, did not deviate much from 
the initial reference. Small variations were acceptable (less than ~ 3" 
change), and an average reference was defined using the initial and final 
reference values. 


APPENDIX D, CALIBRATION OF FILM COOLING PLENUM 


Although the film coolant temperature was monitored prior to the 
vane plenum chambers, uncertainty in the exit temperature of the coolant 
increases as the wall temperature increases and the blowing ratios ap- 
proach xero. An attempt was made to model those conditions which were 
expected to influence the degree of heat addition experienced by the 
coolant as it passed through the injection holes. 

A simple experiment to approximate fluid heat addition by flow 
through tubes with hot walls consisted of passing heated coolant air 
into a water-cooled vane film plenum. The water maintained a relatively 
constant, cool, wall temperature for various flowrates of heated plenum 
air. A stainless steel, plenum supply-line was electrically heated 
using a heater tape capable of reaching 478°K (400°F) at a maximum volt- 
age level of 110 volts AC. A Variac was used to control the voltage to 
the heater tape such that coolant temperature, measured just at the 
entrance to the plenum, could be varied from ambient conditions to 367°K 
(200°F) . A fine, open-bead, corper-constantan thermocouple was attached 
to the outside of the vane wall with the couple positioned at the exit 
of a coolant hole. This exit temperature thermocouple was small in com- 
parison to the coolant hole diameter and was positioned at the centerline 
of the hole. The vane wall temperature averaged 302°K (83.9°F) and 
varied only 0 . 2 % from this mean over all coolant flowrates and temper- 


atures . 
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A prediction of the heat addition thrqqgffe ; the cool ant holes was 
made using a simple heat balance on the flow through a tube of constant 
diameter and wall temperature. Entrance effects were considered con- 
cerning the local heat transfer coefficient using a relation suggested 
in (79) 

if r\ q ha *^*r 0.18 drt 0.2 

h x = (0.036) X (Re. )°* 8 (Pr) 0 * 4 (X) {X (D-l) 

X a 0 Qg x 

where k = thermal conductivity of coolant 

Tg = bulk coolant temperature 


The differential equation for the bulk fluid temperature passing through 
a tube with a wall temperatures T^> and heat transfer coefficient as in 
Eq, (D-l) is 

dI R 

m C C p c dx"~ = " V^O ^ D “ 2 ^ 

where m c = coolant mass flow rate 

C n - specific heat at constant pressure for the coolant 
P C 

T B 0.18 

Assuming that (*==•) in Eq. (D-l) does not differ significantly from 
l W 

unity over the temperature ranges of interest, Eq. (D-2) can be solved 
with the following result 


where 


- 0.8 

T w " t b _ ." C 1 X 

T w - T C,1 ■ 

Tg . = coolant bulk temperature at inlet (x - 0) 
x - dimensionless distance (x/dg) 


{ D-3) 


and 


.Vi? 


■Li 


22 


0.2059k, 


and 

C 1 " * 0.2 „ 0.8 


m c u c 

where 

k c = thermal conductivity of coolant 


= dynamic viscosity of coolant 


Equation (D-4) assumes that both the Prandtl number and specific heat 
do not change significantly over the temperature range of interest. 

Figure D1 shows the variation in the coolant temperature rise ratio 

* 0.8 , . 

with the parameter C^x for x" =4.4 (value for x « L' for this test). 

-4 

The data is shown for five coolant flowrates ranging from 5.45 x 10 
to 4.32 x lo ~ 3 kg/sec (0.0012 to 0.0095 Ibm/sec). The prediction using 
Eq. (D~l) for h v can be adjusted to fit the data provided the 0.036 

A 

constant in Eq. (D-l) is increased by a factor of 8. Such an adjustment 
is not unreasonable considering the short length-to-diameter ratio of 
the hole ( L ' / d Q « 4.4). Figure D1 is the calibration curve used to 
correct the inlet film coolant temperature for the present study. Over 
the range of coolant flowrates investigated, the change in the exit 
coolant temperature represented, at most, a 7 % change in 0£ from a mean 
value of 1 .03. 


Figure Dl. 


Calibration of Film Coolant Plenum 
For Heat-Pick Through Coolant Holes. 
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